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Preface 



This lecture notes covers the principal elements of classical electromagnetic 
theory embodying Maxwell's equations with applications mainly to situations 
where electric charge can be treated as a continuous fluid. The intention is 
to introduce students to the background of classical field theory and the 
applications of the electromagnetic theory to solid state physics, classical 
optics, radiation theory and telecommunication. 

The goal of this course is to provide a compact logical exposition of the 
fundamentals of the electromagnetic theory and the applications to various 
areas of physics and engineering. The treatment is quantitative throughout 
and an attempt has been made to imbue students with a sound understanding 
of the Maxwell's equations and with the ability to apply them to modern 
problems in physics. 

The organization of the lectures is fairly standard and includes vector 
analysis, electrostatic, magnetostatic, mathematical techniques in the solu- 
tion of the Maxwell's equations and the Laplace's equation, time varying 
fields and applications of the solution of the Maxwell's equations. The ma- 
terial on vector analysis gives greater emphasis to the relationship between 
fields and their sources. 
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1 The Classical Theory of the Electromag 
netic Field 



Classical theory of the electromagnetic field or Classical Electrodynamics, 
formulated by Maxwell more than 100 years ago, is now a well established 
theory. In this context 'classical' means 'non-quantum', but we would 
like to point out that the basic equations of electromagnetism, the Maxwell's 
equations, hold equally in quantum and classical field theory. 

Electromagnetic interactions are ONE of FOUR fundamental types: 



Type of interaction 


Relative Strength 


Strong interaction (nuclear) 


1 


Electromagnetic 


10-2 


Weak interaction (e.g. P decay) 


10-12 


Gravitation 





At our level of discussion there is no relation between these four types of 
interaction (i.e. they cannot be considered as different manifestations of a 
single FORCE). 

Classical EM Theory is particularly interesting because of: 

• Its historical role in the development of Physics in the late 19th and 

early 20th centuries. 

• Numerous applications on 'our' scale of existence: 
Chemical bonding is due to EM forces. 

Communications via EM Waves (radio, TV, telephony, computers). 
Generation and transmission of electrical power. 

Knowledge of the rest of the universe from reception of EM waves. 

1.1 Elementary Aspects of Electromagnetism 

• EM forces are due to ELECTRIC CHARGE which is NOT in turn 

explicable in terms of anything else. 
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Charges are of two kinds called positive and negative. In the static 
limit like charges repel and unlike attract. 

Charges are quantized in units of e :^ 1.6 x 10~^^ [Coulombs]. 

In the static limit the inverse square (Coulomb) law of force holds: 



1 g2gl 



that the charge qi acts on the charge q2 with the force Fq^ . The parameter 6:0 
determines the property of the medium and is called the electric permittivity. 

The Coulomb law holds only for charges whose the spatial di- 
mensions are small compared with the distance separating them. 

The Coulomb law can be tested to great accuracy indirectly by showing 
that no charge rests on the inside of a statically charged hollow conductor^. 

If the exponent in the Coulomb's law 

1 



E 



is not n — 2 but n — 2(1 ± e), the potential inside the hollow conductor 
would be large. Since the potential inside the hollow conductor found by the 
experiments was less in magnitude than a small detectable potential, then 

€ ~ 10~^, the level of sensitivity of the detector. 

• Electric charge is conserved (algebraically) 

q — constant 

whole universe 

In electromagnetic theory none of these things are explained in terms 
of anything else. 



^S.J. Plimpton and W.E. Lawton, Phys. Rev. 50, 1066 (1936). 
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• In the N ON- STATIC case (moving charges) the force is no longer given 
by Coulomb's Law. In general is given by the Lorentz equation: 

Fq^ =q2{E + vxB) , 
where E and B are the electric and magnetic fields due to qi. 

The fields E and B will depend on the frame of reference of the observer 
(F must follow the required relativistic transformation law). 
However, we do not think that q depends on the frame of reference. 

i.e. q does not depend on its velocity with respect to an observer. 

This is because in ordinary matter electrons move much faster than ions, 
their speeds depend on temperature, and electric fields are not observed to 
arise from changes in temperature. 

Why there must be a 5 and how E and B are computed for arbitrary 

motion of charges is the substance of electromagnetic theory. 

THE BASIC IDEA OF ELECTROMAGNETISM IS: 

ELECTROMAGNETIC FORCE ON 

CHARGE 1 ^ FIELD ^ CHARGE 2 
Important conclusion 

Fields are generated by charges - NOT by other fields. 



1.2 Macroscopic Charges and Currents 

We know that electric charge is quantized. The electron is a point charge on 
the smallest scale measurable. We may then speak, on a subatomic scale, of 

a microscopic theory of electromagnetism. On a subatomic scale there must 
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be very strong and rapidly varying electric and magnetic fields on spatial 
scales ~ 10~^ m and temporal scales ~ 10"-^^ s. 

When we measure the fields around a macroscopic circuit, clearly we are 
not looking at these fields. We are measuring fields on distance scales ^ 10~^ 
m and time scales ^ 10~^^ s. The microscopic fields sum to (almost) zero. In 
the macroscopic context it is convenient and justifiable to regard the charge 
as a continuous fluid. 



Charge density 

When we encounter a large number of point charges in a finite volume, it is 
convenient to describe the source in terms of a charge density, defined as 

where Ilq is the algebraic sum of the charge in the volume AV, 

The limit is not to zero but to a AV ^ atomic scale size, which is still 

very small on the laboratory scale. 

If the charge density is represented by a continuous function the total 

charge Q in a volume V is given by 



Q = [ pdV , 
Jv 



Current density 



For many purposes it is necessary to introduce the idea of current density. 
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Let Ia = dq/dt is a current through the area A, Then the current density 
is defined by 

r 1- / . 1. Sq ^ pAv6t ^ - 

J = hm —n = iim - — -n = lim — - — : — nJ — pv , 
A 5tA 5tA ^ ' 

where v — vh. 

The hmit is taken in the same sense as for p. 

Total current through a surface area 

If the current density J is known at every point of an arbitrary surface, 
we can obtain the total current through the area. 




The current through 5 A is: 

61 ^ Jcos{e)5A^ J'fidA . 

Then the current through the total area A is the sum of the contributions 
from all elements of the area: 

Ia^ J J'fidA = J J^dA, 

A A 

where dA is a vector representing the element dA of the surface A, In vec- 
tor analysis it is common to represent a surface by a vector whose length 
corresponds to the magnitude of the surface area and whose the direction is 
specified by the unit vector h normal to the surface. 

In summary, when EM Theory is formulated in terms of p and J as the 

sources, we speak of a MACROSCOPIC THEORY. 
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2 Mathematical Description of Vector Fields 



The study of electromagnetic theory requires considerable knowledge of vec- 
tor analysis. In this lecture, we will introduce vector operations we will need 
for our study of electromagnetic theory. As we shall see, based on this lec- 
ture, it is possible to considerably simplify the formulation of electromagnetic 
theory. 



2.1 Gradient of a Scalar Function 



Let us suppose that $ represents a scalar field and that $ is a single valued, 
continuous, and differentiable function of position. 
The gradient of the scalar function $ is defined as: 

grad $ = V$ = ^— n , 

as 

where n is a unit vector in the direction the rate ^ has its maximum value. In 

OS 

other words, gradient tells us in which direction the change in $ is maximal. 
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For some other direction dX: 

= ' dX = — n ' dX = — cose dX 

as as 

In rectangular (cartesian) coordinates: 

(V^)x = V$ • i = — — n • i = — — cos 6> 

as a5 

= lim = lim — = — 

6{s/ cos 9) Sx dx 

Hence 

„^ d^^ (9$^ 5$,^ 
ox ay az 

in cartesian coordinates. 

The gradient is analogous to multiplication of a vector by a scalar. The 
result, of course, is a vector. We do not usually take a gradient of a vector, 
the result would be a tensor. 

Example 

Consider a scalar function $ = xyz. Gradient of $ = xyz is 

dx ^, dy ^ dz - 

A\ A\ ^ 

— yzi + xzj + xyk . 

2.2 Divergence Function 

The divergence is the scalar which results from operation of V upon a vec- 
tor F in a fashion analogous to the dot product of two vectors. In Cartesian 
coordinates: 

, - _ - dFx dFy OF, 
dtvF = W'F=--^ + —^ + —^ 

ax ay az 

Where there is a positive divergence, there is a source of a vector field. 
When V • F = 0 everywhere, the field F is called solenoidal. 
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Example 

Consider a vector F — xi, 

- dx dx , dx . - 

d^.F=-(^•^) + -0.^) + -(fc.^) = l. 

— » A. 

Consider now a vector F = yi. In this case 

Thus, divergence of a given field is different from zero only if the field am- 
plitude changes in the direction of the field. So the divergence is related to 
how the field changes as you move in the direction of the field. 

2.3 Gauss' Divergence Theorem 

Gauss' law, or sometimes called as the divergence theorem, is stated as 

J ^ ^FdV = ^ F'hdS 

V s 
volume integral closed surface integral 




Remember that h is the unit outward normal over S, 
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Definition of Flux 



The FLUX of F through a surface 5, not necessary a closed 5, is defined as: 

s 

Associated with this is the pictorial representation of fields by lines of 
force^ the direction of F given by the tangent to a line and the strength of F 
given by the fine density per unit area. 

Note that from the divergence theorem: 



F = lim 


[dV !i 








— lim 


dV-^0 








[dv\ 



i.e. the divergence of a vector field is the emanating flux per unit volume. 




A non-zero V • F is then implies a source (if positive) or a sink (if neg- 
ative). And if V • F = 0 there is no source or sink - the field lines have no 
beginnings or ends. 

2.4 The Continuity Equation for Electric Current 

Suppose we have some 'substance' of density p in a volume V enclosed by a 
surface as shown in Fig. 1. 
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Let V — macroscopic velocity of 'substance' . 

Let q — rate of production of 'substance' per unit volume . 

Let I — rate of annihilation of 'substance' per unit volume . 




Figure 1: 



Then, the rate of increase of total 'substance' in the volume V 
—Rate of production - rate of annihilation - rate of transport out through 



^JpdV = j qdV - J IdV - ^ pv^ndS 



We can express this as 

d 

V V V s 

Using the Gauss' divergence theorem, we can write 

^ pv ' h dS ^ j ^ ' (P^) • 
s V 

Thus 

dp 



dV = JqdV- JldV- jV'{pv)dV 



dt 

V V V V 
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Since this relation holds for arbitrary V: 



For a 'conserved' quantity such as electric charge: g = / = 0, and then 

pv = J , 

Hence 

which is well known as the continuity equation. 

When stationary currents are involved, then dp/dt = 0, In this case 
V • J = 0, that is for stationary currents the current density is solenoidal. 



2.5 Curl (Rotation) Function 

Curl (or rotation) is a vector which results from the operation of V operator 
upon a vector in a fashion analogous the cross product of two vectors 

, ^ ^ ^ . OF ^. OF ^ OF 
curl F =\/ X F ^ I X — — h J X — — \- k x 



dx 



dy 



or 



VxF = 



OF. OF. 



y 



dy 

i 

A. 

dx 



J 
d_ 

dy 



dz 

k 

A. 

dz 

F, 



dF^ OF, 



dz dx 



j + 



dz 



dFy OF 



X 



dx dy 



k 



Curl is nonzero when the field increases (or decreases) in a different di- 
rection that the field pointed. If the field is pointed in the same direction as 
that in which is increased, the curl is zero. So the curl is related to how the 

field changes as you move across the field. 



When V X F = 0 everywhere, the field F is called irrotational. 
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2.6 Stokes' Theorem 



Stokes' theorem may be stated in the form 

^F'dl = Jv xF'hdS 

I s 



complete loop 



surface not closed 




However 



curl F 



V X F ' fi = lim 



VdS 



This gives an intuitive meaning to any component of V x F in terms of 
the line integral around a small element of surface. V x F is a measure of 

the vorticity of the field. 



2.7 Successive Application of V 

We can introduce scalar and vector products in which the operator V ap- 
pears more than once. For example, since the gradient of an arbitrary scalar 
function is a vector, we can take the divergence of the gradient 



V • = V • 



dV^. dV^. dVA d^V d^V d^V 

I + —J + —k 1 = + + 



dx 



dy 



dz 



dx^ dy- 



dz- 
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The same result is obtained if we take V • V as a new operator with 

properties 

V = V • V = \ \ . 

dx'^ dy'^ dz^ 

The operator is called the Laplacian and is a scalar. 

The Laplacian may also be applied to a vector, with the result 

^2f^ d^F d^F d^F 
dx^ dy^ dz^ 

It is also possible to form the curl of the gradient, which is identically zero 

V X Vy = 0 . 

The divergence of the curl of a vector is also identically zero 

V' V X F = 0 . 
The proof of the above properties is left to the students. 

2.8 Electromagnetic Field Equations and Electric 
Potential 

We will now illustrate some properties of the successive application of V, 
which will allow us to introduce the concept of vector and scalar potentials 
to the electromagnetic field theory. 

2.8.1 Maxwell's Equations 

Consider coupled differential equations for vector fields E and 5, the Maxwell's 
equations: 

V-E = f , (1) 
V • 5 = 0 , (2) 
V.E = --, (3) 

^ ^? - 1 

& ot 
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By means of divergence theorem and Stokes's theorem, they can be easily 

transferred into an integral form 



9_ 

So 



J E'fidS = 
s 

J B ' ndS = 0 , 
s 

j>E^dl = -^jB-ndS, 
I s 

j> B • fidS — fJ^o J J' ^dS + 
I s s 



According to the Helmholtz Theorem, an arbitrary vector F can always be 
written as 

1 r V ' F 1 r V X F 

F = V / dV + --Vx / dV 

Air Jv r An Jv r 

= Fi + Ft. 

Thus, specification of V • F and V x F is necessary and sufficient to deter- 
mine F. 

Hence, we need four equations of this type to determine E and B, 

In the following lectures, we will discuss how Maxwell arrived at these equa- 
tions. We will also discuss the effect of Einstein's special theory of relativity 
on how we think about electromagnetic fields and Maxwell's equations. Much 
of our discussion will be about 'how to solve Maxwell's equations'. 



2.8.2 Electric Potential 

The solutions of the electromagnetic field equations are not in general sim- 
ple and straightforward. Often their solution is aided by the use of potentials. 

A potential is a quantity from which a vector field can be derived by some 
process of differentiation. 
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Examples 

(1) Let V X F = 0. Since V x V = 0 then the equation V x F = 0 
will have an integral of the form: 

and $ is called a scalar potential. Thus, the field may be derived from the 

gradient of the scalar potential function 

In the electrostatic field V x E = 0. Hence, E = — V$. ^ Then the 
work done per unit charge q moving from point A to point B is: 




AB ' 



A 



The minus sign is inserted in the definition of $ to agree with the definition of $ as 
a potential ENERGY. The negative sign can also be understood physically from the fact 
that E is in the direction that a positive charge moves, hence in the direction of decreasing 
potential. 



23 



Note that if V x E = 0 then: 

i E • dl = Jv xE- MS = 0 . 

Hence 



1 (AB) 




2(BA) 



The work done is independent of the path chosen. Thus a field F with 
V X F = 0 is a conservative field of force. 

(2) Let V ' F = 0. Since V • V x = 0 the equation V • F = 0 will 
have an integral of the form: 

F = V X A 
and A is called a vector potential. 

Note: 

It turns out that V • 5 = 0 always in electromagnetism so there will 
always be a vector field A such that 5 = V x ^1 and such an A is referred 
to as THE vector potential in electromagnetism (though there may be other 
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electromagnetic field functions with zero divergence). 

Note also that just writing V xA = B does not completely specify A even 
if B is known everywhere. One needs to specify V • ^ as well to completely 
determine A, Equivalently we can say that defining V x A = B still leaves 
us free to define V • A, 

On the other hand V x E ^ 0 in general in electromagnetism so it is not 

in general possible to write E = — V$. 

Exercise in class: Vector analysis 

For a given scalar function $ and vectors A, 5, (7, indicate successive steps 
you would follow in the calculation of the following expressions: 

(i) VxV'A; A^VB] AV'B; AxV^B] A^V^ ] 

WAxB'C ] V^'AxB] X V$ ; V-AxVE-C ] 

V X A'BxC , 



(ii) Which of the expressions in (i) are vectors? 

(iii) Which of the four statements is correct: 



1. 


V^A = 


(V • V)/ , 


2. 


V'l = 


V • (Vl) , 


3. 


V'l = 


V • (V • 1) , 


4. 


v'l = 


V X (V • 1) 
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Weekend exercises 

(a) For fields of the form r^r^ (r / 0), find for which values of n the 
divergence is zero. 

(b) For fields of the cylindrical form p'^(j)^ (p 7^ 0), find for which values 
of n the curl is zero. 

(c) If the potential $ satisfies the equation (Laplace equation) = 0, 
show that V$ is both solenoidal and irrotational. 
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3 The Experimental Basis of the Develop- 
ment of Electromagnetic Theory 

In this lecture, we will present the basic properties of the electrostatic field 
in vacuum. The Coulomb law for the force between two point charges is the 
experimental basis for the development of electromagnetic theory. 

3.1 Coulomb's Law — Force between Static Charges 

In 1785, Coulomb investigated the nature of the force between charged bod- 
ies, and the results of his experiments can be formulated mathematically in 
what is known as Coulomb's law 

^2 = T- ^ • 

47r6:o 

We can write Coulomb's law as 

Hence 

^1 = ~A 2 ^ • 

The electric field is an example of a vector field. In principle, we can al- 
ways calculate an electric field using Coulomb's law. However, there is an 
alternative way we can find the electric field. In particular, the field may 
be represented by means of the flux concept. The total flux of E from a 
point charge qi may be readily calculated by integrating E-dS over a surface 
enclosing qi. 



3.2 Derived Result — Gauss' Law 

Consider a macroscopic charge qi closed by a surface S. We will show that 
the flux of the electric field produced by the charge qi is proportional to the 

charge gi, and is independent of the shape of the surface closing the charge. 
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Consider first tlie flux througli a spherical surface: 

SI SI ^ ^ 

dS — ^ Airrl = — 



Next, consider the flux through an arbitrary surface 

S2 52 ^ 

where 5^2 is an arbitrary surface enclosing qi. 

But from the inverse square law and some geometry: 

d^S2 = -r^r-hdS 

qi dS cos 9 qi „ 

where dQ is the solid angle subtended by dS at qi. 

The element of solid angle d^l is independent of where we cut the bundle 
of electric lines of force. 
Thus 

d'ifs2 = d'ifsi and f d^s2 = <f d'^si = — . 

J J Sq 



28 



Furthermore for some arbitrary number of charges 9i, 92? 93? • * * ? 9n within 

a surface 5": 




jd^ ^ j{Ei + E2 + Es + -')'hdS 
q1.q2.q3. 1 

= — + — + — + ••• = —2^g , 

^0 ^0 ^0 ^0 

where Z) 9 is the algebraic sum of all charges within the surface S. 

If q is outside 5", the surface integral vanishes since the total solid angle 
subtended at q by the surface is zero. Thus 

+ (i^2 = 0 . 



Proof: 



q f ' hi dSi 
f ' hi dSi dSi± 



Vi rf 



= -dn . 



Similarly 



d^2 = -T^dn = -d^i . 
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Since 



f • n2dS2 , 
5 — = +ail , 

we find that d^i^i + (i^^2 = 0 and integrating over all S: 



^ = rf* = 0 . 

The physical interpretation of this result is that field lines originating from 
an external charge and entering the surface S must also leave this surface. 

In summary, the Gauss' law says that the total electric flux through a closed 
surface S is: 



* = — , where q = J2 charges INSIDE S . 

Using the definition of the flux, we often write the Gauss' law as 

<f E'dS = — , 

Js Sq 

The power of the Gauss' law lies in the fact that we are free to apply 
it to any closed surface whose shape can be chosen arbitrary such that the 
evaluation of the surface integral becomes a simple straightforward task. The 
Gauss' law is particularly useful in simplifying the calculation of the electric 
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field produced by certain symmetrical charge distributions. We illustrate this 
in the following examples. 

Example of an application of the Gauss' Law 

An infinitely long line is positively and uniformly charged with a constant 
linear charge density pi. Use (a) Coulomb law, (b) Gauss' Law to find the 
electric field about the line. 

(a) If we have to calculate the field due to a static macroscopic distribu- 
tion of charge using the Coulomb law, we divide the macroscopic charge into 
infinitesimal (point) charges dq which produce an electric field dE, The field 

dE is given by the Coulomb field 

dE = -L-%r . 

47r6:o 

Then the total electric field is found by vector addition 

i = j dE , 



Since we are adding vectors, a caution must be employed. We use the fol- 
lowing procedure, which is general and can be employed to any system: 

L Write the expression for the electric field dE produced by the infinites- 
imal (point) charge dq, 

2. Resolve this vector into components dE^^dEy^ and dE^^ 

3. Calculate each component of E by integration, e.g. E^ = J dE^. 

4. Find the resultant E from its components 

E = eJ + Eyj + E^k . 

Return now to our example of the charged infinitely long line. 
Take a small element dl of the line containing a point charge dq. Electric 
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field produced by the point charge at A distance r from dl is given by the 
Coulomb field 



We see from the figure that 



h 



r — 



I — h cot 9 



Hence 



dl — %r7A^ ? ^9 — Pi^^ ? 

sin 9 



and then 



dE = 



Pi 



Pi 



h \ sin^ 9 / 



, ( cos + sin 9j j , 

Aire oh ^ / 

where we have decomposed the unit vector f into two {x^ y) components 



f = COS 6i + sin Oj . 
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Integrating over 6 from ^ = 0 to ^ = tt, as the line is infinite, (at x = 

— oo, ^ = 0 and at a; = oo, 9 = it), we obtain 



E 



Pi 



(— sin 0 + sin tt) i + (cos 0 — cos tt) j 



_ Pi 

Thus, the electric field produced by the charged line depends inversely on 

the distance from the line and points in the direction perpendicular to the 
line. 



(b) Let us now calculate the field by the direct application of the Gauss's 
law. 




The electric field near the uniformly charged line must be radially directed 
because of the symmetry of the problem. The field must have cylindrical 
symmetry because the problem is unchanged by rotating the line about its 
axis. The field must also be independent of position along the line because 
the distance to either end is infinite. This is the ideal situation for the 
application of Gauss's Law. We can apply a cylinder surface of radius h and 
length L centered about the line of charge, see the Figure above. 

According to the Gauss's law 

Js So 

where q — piL\s the charge closed by the cylinder surface. 
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The flux through the cyUnder surface spUts into three fluxes 



E^dS= / E^dSA+ / E'dSB+ / E'dSc ^ 
s J A Jb Jc 

Since E _L dSA^ E _L dSc^ E \\ dSs-, and the magnitude of E is constant 
along the surface 5, the flux through the cyhnder reduces to 



[ E'dS= [ EdSs = 27rhLE , 

Js JB 



i.e. a nonzero flux exists only through the surface 5, and since the cylinder 
is symmetrically positioned about the line of charge, the magnitude of E is 
constant over the surface B. 

Then, according to the Gauss's law 

2ME = ^ , 

which gives 

Pi 



E = 



271 e^h 



Note how simple are the calculations of the electric fleld using the Gauss's 
law. However, we were able to solve this problem because we knew the di- 
rection of the fleld at any point around the line. 



3.3 Biot-Savart Law 



In 1819, H.C. Oersted showed the magnetic effects of electric current. 
In 1827, Ampere showed that quantitatively the magnetic forces in macro- 
scopic circuits can be accounted for by what has come to be known as the 
Biot-Savart Law. 

Experiment 
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Two long parallel wires. Force between the wires: 



where fio = 47r x 10~^ [H/m] in SI units, is the permeability of the vacuum. 

(This defines the SI unit of current, the Ampere.) 

If /1II/2 then the force F is attractive. 

If Ii anti\\ I2 then F is repulsive. 

If one wire is rotated through 90° then F = 0. 

All such observations are explained by the law between current elements 
being: 

_ /^o l2d£2 X {hdii X f) 
We can write the force as 

dF2 = I2di2 X dB , 

where 

,5 ^ohdlixf 

^ ~A 2 ' 
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which is known as the Biot-Savart law for magnetic field produced by the 
current element Iidii, 

The Biot-Savart law allows to compute magnetic field produced by an 
arbitrary current distribution Idi 

- fiol fd£ X f 

B — 



47r 

The method requires integration of small current elements. 

We can simplify the calculations of B by using the following procedure. 
If we replace r/r^ by —V{\/r)^ the integrand becomes 

-dix V(l/r) . 
However, using a vector identity that 

V X 1 — 1 = -V X rf/- dix V (-\ = -dix V (- 
r r \rJ \r 



since V x rf£ = 0, we can write 



5 = V X 



fiQl fdi 



Since we can write 



/- ■ w 

Ji r 



B = V X A, 
as V • 5 = 0 always, we can first calculate A: 



A = 



fiol f d£ 



Att 



I- (7) 

Ji r 



and then using (6), we will find J3. 

The integral for A is easier to calculate than the original expression (5) 
for B, Since the curl operation is readily performed, we may use (7) as an 
intermediate step for finding 5 in a simpler way. 

The vector A is called a vector potential, and will see later in the course 
many useful applications of A in electromagnetic theory. 
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3.4 Current Element and Charge Element 



A charge dq moving with a velocity v is equivalent to an element of cur- 
rent Idi: 

dq — idt = 1 — , 

V 

where dt — the time for all the charge in di to pass out of the volume. 



I 




dl 



Hence 



dqv — Idi . 

The Biot-Savart Law can then be written as 

dF2 = ^dq2V2 X ^ , 

dF2 — dq2 V2 X dBi , 

with 

dBi = —dqi — , 

which shows that magnetic field is produced by a moving electric charge. 
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3.5 The Lorentz Force 



The Lorentz force is obtained putting Coulomb's Law and the Biot-Savart 
Law together: Fe = qE and Fm — qv x B: 

Fem = FE + FM=q{E + vxB) , 

Thus, a motion of electric charges is modified by both the electric and mag- 
netic forces. If the charge is stationary, the force depends only on if it 
moves, there is an additional force proportional to v. 



3.6 Amperes Circuit Law 



Amperes circuit law is a useful relation between currents and magnetic fields. 
This law allows us to calculate magnetic field produced by some currents in 

a very effective way. 




The Amperes law says that for an arbitrary closed path around a current 
carrying conductor, the component of magnetic field tangent to the path is 
proportional to the net current passing through the surface bounded by the 

path 

^B'de = fiol 
closed loop / = total current through loop 
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Proof: 



Consider a long wire of radius a carrying current /. Let P is a point on 
the integration path, see Figure 2. 
The magnetic field at P is: 













1 









Figure 2: The source circuit and the integration path to prove the Ampere 
law. 



s 

Moving P by d£ is equivalent to moving the current circuit by —d£. 
The solid angle subtended by —di^dsai P is: 

{—di X ds) • r —di • ds x r - ds x — f 

/^Zi /^Zi 

(The element of area normal to f is —di x s- r) 



39 



Thus due to the path element di^ the change in soHd angle subtended 
at P by the circuit is: 



X (— f) 



s 



Hence 

^ Att f ^ Att - 
dn = d£' — (tdB = —B • di , 

where integration is around the circuit s giving the magnetic field B at some 
point P as shown in the diagram. 

Now integrating round the closed path: 

jB^dl^ ^/ jd^ . 

If P moves round a closed path (returning to its original position but not 
circulating through the current loop: 

^dil = 0 . 
But if P circulates through the loop: 



and then 

IB' di = ^ I A7r = fiol . 

We conclude that the line integral of the magnetic field round a closed 
loop path is equal to fiQ I where I is the current passing through the loop. 

The Amperes law can be applied in highly symmetric situations to find 
the magnetic field more easily than by computing with the Biot-Savart law. 
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In either case, the result is the same. In case that lack the proper symmetry, 
Amperes law is not easily applied. 

Example of an application of the Amperes Law 

An infinitely long wire caries a constant current /. Use (a) Biot-Savart 
law, (b) Amperes law to find the magnetic field about the wire. 




(a) The element of magnetic field dB due to an element / d£ at distance 
r = rf is found from the Biot-Savart formula: 

dB = ^IdIx ^ 



Att 

— > ^ 

where we note that all dB are in the same (j) direction normal to the direction 
of the current. So we see from this symmetry that the field lines are circles 
concentric with the current. Furthermore, along any such circular path the 
field is constant in magnitude. 

Let us calculate the magnitude of the magnetic field. Since 

^ ^ h 

di X f — dlsinOcj) ^ r — , I — h cot 0 , 

sin 6/ 



we have 



dl — or7,d6 

sin 0 
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and then 

- HqI sin^ e ( h \ , 2 /io/ . .,^7 

an — — — — TT- sin Odd (b — — - — r sm BoBq) . 

Att \ sin^ 9 J Airh 

Integrating the above equation over the length of the wire, we obtain 

B — — / sin 9d96 — (cos 0 — cos tt) 6 
Airh Jo Airh ^ ^ 

- ^^^ 6 



(b) Let us now calculate the field using the Amperes law. 

Since the field lines are circles concentric with the current, and along 
any such circular path the field is constant in magnitude, this is the ideal 
situation for the application of Amperes law: 

(j>B' di= fiQl 27rhB = /xq/ ^ 5 = ^ . 
J 27rh 

Note how simple are the calculations of the magnetic field using the Am- 
peres law. 



3.7 Faraday's Law of Electromagnetic Induction 

Faraday discovered electromagnetic induction by changing magnetic field. If 
we consider a closed stationary circuit located in a varying magnetic field, 
the induced electromotive force around this circuit is equal to the negative 
time rate of change of the magnetic flux through the circuit 

d^ 



where $ is the total magnetic flux through the circuit. 

From the deflnition of the flux 

$ = J d~S , 
s 
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and from that emf S = work done per unit charge, we have 

S 

q 



- = ^{E + vxB) ■ di. 



e 



In a stationary circuit v = 0 (and anyway v x B A. d£ since v || d£). Thus 
vx B-d£ = 0. 



Hence 



€= 4e - di 



and finally 



E ■ dl= -4- [ B - dS , 
dtJ 

I s 



or 



i s 
This is the Faraday's law written in the integral form. 

Using the Stokes's theorem, we can rewrite the Faraday's law in the dif- 
ferential form 

^ ^ dB 
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The Faraday's law tells us that time- varying magnetic fields give rise to 
electric fields. This shows that the fields are related to each other, and we 
then must speak of electromagnetic fields, rather than separate electric 
and magnetic fields. 



Questions: 



1. Prove that the total electric flux through a closed surface S is propor- 
tional to the total charge inside the surface. 



2. Prove the Amperes circuit law. 



3. Derive the integral form of the Faraday's law and then transform it 
into the differential form 

^ ^ dB 



Weekend exercises 



(a) Find the pressure (force per unit area) between two infinite and op- 
posite charged parallel planes of charge density ±cr. 

— » ^ 

(b) Describe the magnetic field associated with E = xj. 
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4 Differential Equations for the EM Field and 
Maxwell's Theory 



We know now that electromagnetic forces are carried by electromagnetic 
fields that propagate at speed c 3 x 10^ ms"-*^. Because of the finite prop- 
agation speed we are forced to assign energy & momentum to the fields i.e. 
we must think of them as real physical entities as against mere mathematical 
conveniences (as is the case for static fields). An electromagnetic system then 
qualifies as static only if all the charges have been at rest longer than the 
time taken to traverse the system at speed c. 

In the 1830's Michael Faraday carried out experiments to measure a finite 
electromagnetic propagation speed. He was unsuccessful due to lack of time 
resolution in his apparatus. Faraday would have had no reason to think that 
the electromagnetic speed was the same as the speed of light (then known). 
In the ISGO's, James Clerk Maxwell, seeking to advance Faraday's ideas about 
electromagnetic fields, by a brilliant process of intuition worked out how to 
generalize certain differential equations deduced from static experiments. He 
produced a set of field equations known by his name today. Maxwell also 
had a theory i.e. a set of qualitative ideas underpinning his equations. The 
theory (unlike the equations) has not stood the test of time. 



4.1 Differential Equations for the EM Field 

Let us take as the source of the electromagnetic field a continuous charge 
and current distribution (p and J). Then we will produce macroscopic field 
equations 

q = J pdV , 1 = j J'fidS . 

We find differential vector equations from the integral forms of observational 
results discussed above. 
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4.1.1 Divergence of E 



( E-dS = — = f ^dV 



Applying Gauss' Theorem to Coulomb's Law: 



j V'EdV = f —dV . 



^0 

Since this must hold for arbitrary V no matter how small: 



£0 



4.1.2 Curl of E 



Prom Faraday's flux cutting rule we had: 



J E-di= - I ^-ndS 



dt 

Applying Stokes's Theorem: 

V xE-ndS = - I ^-hdS 

J dt 

For arbitrary S and n implies 

dB 



V X E = - 



dt 
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4.1.3 Divergence of B 



From Biot-Savart law, we had: 

Air 

Applying the rule 

V- (Ix 5) = 5 ♦ V X 1- 1* V X 5 , 

we obtain 

V . 5 = ^ 1 4 • V X rf/- dh V X 4 I 
However, V x di 



Hence 

V-5 = 0 . 



= 0, since d£ is a constant vector, and 

/% 
r 

V X — =0 for any n , 



4.1.4 Curl of B 



From Ampere's Circuit Law: 



J B ' d£ = fiQl = /io J J-ndS 



Applying Stokes's Theorem, the relation 



J V X B ' h dS — fiQ J J'hdS 
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holds for arbitrary n. Thus 



V X 5 = fioJ . 



However, Maxwell realized that unlike the previous three differential 

equations, this one could not be generally true. 

To see this, take its divergence and remember that V • V x F = 0 for any 

vector function F: 



V • V X 5 = 0 = /xoV • J 



Thus 



V-J = 0 !!! 



We have already seen that conservation of electric charge requires 



dp 

di 



Thus V • J = 0 implies that 



at ~ ' 



i.e. we could never charge or discharge a capacitor. 



Maxwell guessed at the right value of V x 5 as follows: 
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Since 



^ dp ^ 
V-J+Tr-=0 and y9 = 6:0 V • E 
at 

(from the first Maxwell equation for V • E) 

Hence 

V.J + -^=0 or V.(^ + .„^)=0 

If then we (after Maxwell) write: 

^ dE 



instead of 



V X 5 = fioJ 5 

we obtain 

^ dE 

V-Vx5 = 0 ^ V'(J + £o^)=0, 

which is in accord with conservation and motion of charge. 
Note from above that the term that Maxwell added 

dE 

has the dimensions of current density. Maxwell called it the displacement 
current density. 

Maxwell had a theory underpinning his equations in which the displacement 

current was a real physical current - due to 'polarization of the electromag- 
netic ether'. This theory has not survived. Nevertheless the above term is 
still referred to as the 'displacement current'. 

We write the fourth Maxwell equation then as: 

^ r dE 

at 
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The differential form of Maxwell's equations is easier to interpret phys- 
ically and is also useful in deriving the boundary conditions that the field 
vectors must satisfy. 

The Maxwell's equations are self-consistent and no experimental evidence 
for requiring any further modifications has been found. 

Exercise in class: Fields within a capacitor 

A plane parallel capacitor is being charged with a current /. Show that 
the displacement current between the plates of the capacitor is equal to the 

conduction current / in the external charging circuit. Remember that the 

— > 

displacement current density is, by definition, so the displacement cur- 
rent through a surface S is: 



Assume that the external wires are perfect conductors so that E is zero in 
them. Assume the space between the plates is a perfect insulator so no con- 
duction current flows within the capacitor. 

Can you see any curious consequence in this case if the displacement 
current is assumed to be a real physical current (flow of charges)? 



Id 




dE 



s 
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4.2 Maxwell's Equations and Prediction of Electro- 
magnetic Waves 

The basic equations to study electromagnetic theory are Maxwell's Equa- 
tions, which govern the behavior of the time-varying electromagnetic field 

L V • ^ = 
11. VB- 

IIL V X E 

IV. V X 5 

The first equation is just Gauss's law, the second tells us about the non- 
existence of magnetic charges, the third equation is Faraday's law, and the 
final equation is Ampere's law. 

Maxwell's immediate triumph was to predict the existence of electromag- 
netic waves and their propagation speed. The calculated speed came (within 
experimental error) to be equal to the measured speed of light. This predic- 
tion obviously led to the conclusion that light was electromagnetic in nature. 
Thus arose a synthesis of electromagnetism and optics. 



^ P_ 
= 0, 



dB 



dt ' 



dE 
'dt 



4.2.1 The Wave Equation for EM Waves in Vacuum 

In a vacuum there are no sources i.e. p = 0 and J = 0. 

Hence, the Maxwell's equations III and IV reduce to the following: 

VxJ = --. (8) 

dE 

V X 5 = = eoA*o^ • (9) 

We eliminate E ov B between equations (8) and (9) to obtain differential 

equations for E ox B alone, using where required V • = 0 and V • -B = 0. 
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Method: Think of V x and ^ as hnear (differential) operators. By analogy 

with methods of solving linear algebraic equations, applying V x into (8) and 
d/dt into (9), we obtain 

Vx(Vxi^) = 



Hence 



V X (y xEj = -eofio-^ , 

d^E 



VV-E-V^E = - £0^0 
Because V • -E = 0 in the vacuum, we finally obtain 



„o^ I d'^E 



where 

1 



The parameter c has the dimensions of velocity and is numerically equal 
to 3 X 10^ ms-^ 

Equation (10) is the standard form of a three-dimensional vector wave 
equation. The field B satisfies the same equation. 

4.2.2 Plane Wave Solution to the Wave Equation 

The wave equation in a vacuum is 

for X = E,B. 
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Look for plane wave solutions propagating in the z direction. 

In this case, djdx = d/dy = 0, and then 

_o 



dx^ dip' dz^ dz^ 
The differential equations for E and B both have the same form: 

d^x _ 1 d^x 

dz^ (? dt^ 

Such an equation has solutions of the form X = f{z±ct) where / is an 
arbitrary function. 

This solution represents a signal propagating with speed c as can be seen 
from the following discussion. 

Let Xq — f(zo — cto) i.e. X Oit t — to and z = Zq. 

Now examine X at time At later and distance Az further along in z. Since 

a harmonic wave does not change in vacuum, we have 

Xi = f{zo + Az-c{to + At)) 

— fi^Q — cto) — Xq when Az = cAt , 

i.e. the signal propagates a distance Az — cAt in time At i.e. it propagates 
with speed c. 

Proof of solution: 



Let / represent f{z — ct), 

df df d{z-ct) df 



dz d{z — ct) dz d{z — ct) 

where ' means differentiation wrt z — ct. 
Similarly 

d^f _ df _ df d{z - ct) 
dz'^ dz d{z — ct) dz 
dl ^ df d{z -ct) _ 
at d{z-ct) dt 
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Similarly 



Consequently 



Thus 



1 1 

C" = > c = 



2 



which with the numerical values of the parameters 

£o ^ 8.85 X 10"^^ Fm"^ /^o = 47r x 10"^ Hm"^ 

gives 

c :^ 3 X 10^ ms"^ . 
What determines the form of the function /? 

Boundary conditions do. 
Source (charges & current), propagating EM in a vacuum, detectors (other 

charges) 



4.2.3 Harmonic Waves 

In vacuum, we chose a plane wave representation for the propagating EM 
wave. For a harmonic wave 

c = /A=p ^ = 27r/, ^ = y' 
where / = frequency (Hz), A = wavelength (m), 

uj — radian frequency (radians s~"^), and k — propagation constant (m~"^). 



f{z-ct) = f{z--t) = Mujt-kz) 

54 



A plane wave is represented by 

Note that some textbooks on electromagnetism, for engineers in particular, 
use j instead of i for the imaginary number. We will use i throughout this 
lecture notes. 



4.3 The Transverse Nature of Plane Waves in Vacuum 

We now investigate the relations between the amplitudes and phases of the 
electric and magnetic fields in a plane harmonic wave. While it is true that 
the magnetic field satisfies the same wave equation as the electric field, it is 
not independent of the latter, since one must satisfy the Maxwell equations III 
and IV. 

Since V • 5 = 0 always in electromagnetism, and 

_ ^ 55. dBy , dB. 



dx dy dz 
dB^ 

= 0 + 0 + —^ 

oz 

for a plane wave propagating in the z direction, we have 



dz 



= 0 



However, for a plane wave 



— — = -ikB^ 
oz 



Hence, the rhs must be zero, which means that either A; = 0 (zero frequency) 
or Bz — ^ (transverse wave). 

For a propagating wave A; 7^ 0, so the wave is transverse in B, 

In a vacuum V • ^ = 0 and then by the same argument we conclude that 
the plane wave is also transverse in E. In other cases in electromagnetism 
(e.g. for a plasma in a magnetic field) a plane wave may not be purely trans- 
verse in E, 
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In addition: E J- B for a plane EM wave in a vacuum, 



Proof: 



Consider a harmonic wave propagating parallel to the z axis. 
In this case the field components are ~ Q^i^t-^^) 



We will expand Maxwell equation III 



V X E = - 



dB_ 
~dt 



in Cartesian coordinates remembering that d/dt = iu; 



i j k 

d d d 
dx dy dz 

BEE 



As ^ 

= -iLo{BJ + Byj + B^k) 



(11) 



For electromagnetic plane waves propagating along the z axis in a vac- 
uum: 

d 



E-f — B-f — 0, 



Hence, Eq. (11) reduces to 



9 d ^ 



dx dy 



dz 



= —ik . 



I J k 
0 0 -ik 

Ex Ey 0 



= -iuj{BJ + Byj) . 



Comparing the left and right-hand sides: 

X component 



ikEy — —iojB 



X 



B^ — E, 



—ikE^ — —icuB 



y 



y component 

Consider a scalar product 

E^B = {iEx+jEy)^{iBx+jBy) 

—k k 
= (iEx+jEy)'(i — Ey+j-Ex) 

k k 

E E H — E E =0 

^ ^x ^y ^ 5 



k 



y 



Bq. E^ 



to 
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which means that E -L B, 

Note that 

E _u; _ 

In electromagnetic theory when E and B are related, their ratio is always a 
velocity characteristic of the problem in hand. 

This is as far as Maxwell took the subject. It was for others like Heinrich 
Hertz 1884 to show how to solve Maxwell's equations with source terms p, J 
included (i.e. the generation of electromagnetic waves). We will consider 
this later. 

You should be aware that we have not derived Maxwell's equations from 
the static limits like Coulomb's Law and the Biot-Savart Law. The solutions 
to Maxwell's equations include the static limits as special cases but many 
more. Maxwell's equations have the status of postulates suggested by exper- 
imental results. 

In summary, we have the following important results for related electric 
and magnetic fields propagating in vacuum: 

1. The electric and magnetic fields propagate in a form of plane waves, 
so-called electromagnetic (EM) waves. 

2. The plane EM wave is transverse in E and B, 

3. The electric and magnetic fields are perpendicular to each other. 

4. The ratio E/B is constant and equal to the velocity of the wave, that 
is equal to the speed of light in vacuum. 
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Questions: 

(1) Starting from the Maxwell's equations derive the continuity equation, 
i.e. show that conservation of charge is built into the Maxwell's equations. 

(2) Using the Maxwell's equations show that B satisfies the same wave 
equation as E, 

(3) Show, using the proof of solution of the wave equation, that f{z-\-ct) 
represents a signal propagating in the negative z direction with speed c. 



Weekend exercises 

(a) Demonstrate that the Coulomb field E — -^^^t for a stationary 

point charge, follows from the Maxwell's equations. 
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5 EM Theory and Einstein's Special Theory 
of Relativity 



Special relativity (1905) grew out of Einstein's meditation on electromag- 
netic theory and the properties of space and time. Historically, the insights 
of Einstein's theory follow after electromagnetism. Logically however, special 
relativity contains more general statements about nature than electromag- 
netism. Electromagnetic field theory is just one of a possible set of field 
theories that are compatible with the Einstein theory of space and time. It 
is evident that relativistic effects are important if we were to calculate the 
field of a charge moving with a speed comparable to that of light. What is 
not so obvious is that special relativity offers insights in to aspects of elec- 
tromagnetic theory even in the case of the low speed charges we consider in 
this course. 



Two such aspects are: 



(1) The unity of the electromagnetic field i.e. the field is a single entity 
with 6 components (represented by two vectors E and 5, each with three 

components). 

(2) Understanding the nature of causal relationships in electromagnetic 
theory, e.g. 

^ ^ dB 

Does this mean that a time changing B causes a spatially changing E? 



5.1 Lorentz Transformation Equations for Space and 
Time 

The principle of relativity 

1. The laws of physics are the same in all inertial reference frames. 
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2. The speed of light in vacuum is independent of the uniform motion of 

the observer or source. 

The constancy of the velocity of light, independent of the motion of the 
source, gives rise to the relations between space and time coordinates in dif- 
ferent inertial reference frames known as Lorentz transformations. 

Consider a stationary reference frame S and a inertial frame moving 
with a velocity u parallel to the x axis. 




Let x^y^z are coordinates in the S frame and x^^y^ are the coordinates 
in the frame which is moving with a constant velocity u parallel to the x 
axis with respect to S, 

The time and space coordinates in are related to those in S by the ID 
Lorentz transformations: 



x^ = 7(x — ut) 

y' = y 
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where 




is the Lorentz factor. 

The above transformation corresponds to a situation of parallel to the x 
axis. Later in the course, we will consider the general case of the velocity u 
of the frame in an arbitrary direction. 



5.2 Force Transformation Equations 

We will demonstrate how one might infer the law of Biot-Savart from appli- 
cation of special relativity to Coulomb law. 

Consider a particle which is moving with velocity v = v^i + Vyj + v^k in 
the S frame and is acted on by a force with components Fx^ Fy and F^. Then 
in the frame: 

JP^ — TP V- TP - TP 

2 . 1 



y 



2 X 1 



(1 - F. 

pf _ V QlJ_JP — 



1-^ ' 7(1-^) • 

Suppose Fj^^Fy^Fz represents the velocity independent Coulomb force. 
Then in the frame (source of the field now moving) the force is no longer 
velocity independent. In electromagnetic theory we say that there is now a 
magnetic force and we define a magnetic field B that determines the mag- 
netic force. 

We now present a detailed calculation that illustrates how the form of 

the Maxwell's equations is determined by nature obeying Einstein's special 
theory of relativity. 
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5.2.1 The Force between Two Charges Moving with Constant Ve- 
locities 



Invariance of electric charge 

In ordinary matter, electrons move with much greater speeds than protons 
yet there is no associated electric field. This implies that electric charge is 
independent of velocity unless electromagnetic laws modified in some compli- 
cated way (see discussion by King in Physical Review Letters 5, 562 (I960)). 

Consider charges qi and q2 moving with velocities u and v in an inertial 
frame S, No loss of generality occurs if v is taken parallel to the x axis. 




Now consider a frame moving with velocity v along x axis, that q2 is 
stationary in S\ Assume that at time t = 0, the frames S and overlap. 

From the principle of relativity, in the frame Coulombs law holds. The 
force on qi seen in 5' is therefore 

We shall transform this expression to find the force observed in the S 
frame in which the source of the field (^2) is moving. 

We shall see that what we normally call the MAGNETIC FIELD arises 
as a natural consequence of relativistic invariance with no extra assumptions. 
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The X component of the force is 



and similarly for y and z. 

The Lorentz transformations at t = 0 are 



X = 



z 



^x 

y 

z 



vx 



= -r 



To do the transformations, we also need transformation from r ' to r. Clearly, 
there is an "axial symmetry" here and we should involve the angle 6, 



/2 /2 I /2 I /2 2 2 I 2 I 2 

r = x +y +z = + y + z 



2 / 2 , 2/^ + \ 2 



= r 



7 



x^+[l- 



+ z') 



v 



= 7' - sin" ^ ) = I 1 - ^ sin" 6 ) , 



v 



where sin 6^ 
Hence 



{y^ + z^)/r. 



1 

\ 2 



r' = 7r 1 sin 9 



Substituting into F^, we obtain 
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where 

1 92 



g — 3 

U r^Z^6 M _ £_ gj]^^ 



Similarly 



7 7 



The force transformations are 



a) a; component 



— F - F - I 



(l-^) ' c2(l- 

Thus 



b) y component 



y 



7(1- 



vu 



X 



c2 



= 7(i-5)f; = 7(i- 



c) z component 



Fz = qigz (1 - 



VUx 



Substituting for Fy^^ in equation 



= qiQX + qigy^ + qiQZ— 
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Note: Here is the germ of the magnetic field. The last two terms are typical 
second order relativistic terms ^ jc? , In a nonrelativistic calculation we 



would have — F!^, 



We now combine results a), b) and c) into a single vector equation for 
the force F. First, note that 



Next, we can write the x component as 

vu 



= qigx ( 1 - — 1 + qigx- 



X 



C2 



VUy VUz 



= 1 - 



and with the y and z components 



Fy = qigy-- 

= qigU- 
these three components combine into 



C2 



y 



', z 

C2 



F = gipll--^ lr + gip-(wf) 

= qigf+ [v{u-r) — f{u • v)] 

- ^ Qig ^ ,^ ^ 

r = qigr -| — —u x [v x r) 

& 

We can write this equation in the form of the Lorentz force 

F = qx(E^uxB\ 

where 

1 92 ^ 



E = gf = 
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and 



-* V X qf V X E 



Note that as ^ 0, 7 ^ 1, and then 

il/ — > — T 

Moreover, the ration of magnitudes of magnetic to electric term in the force 
equation is uvjc? ^ i.e. magnetic forces are second order relativistic effects. 



5.3 Electric and Magnetic Field Lines of a Moving 
Charge 

5.3.1 Electric Field Lines 

Let j3 — v/c. Then, we can write the electric field as 

E = -^ ■"^-^'> .r. 

"reo r8(l-^2sin^«)' 

For a given 9^ the electric field E still varies as but the field lines are 

crowded in the direction perpendicular to v. 

In the forward direction 0 — Q and then 



1 q{l-f^') 



^ ^ ~7~Z < ? 



where E^ is the static electric field (at ^ = 0). 

In the perpendicular direction, ^ = 7r/2, and then 



The electric field lines radiate from the present position of the charge. 
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static field 




5.3.2 Magnetic Field Lines 

From the relation between electric and magnetic fields, we find 

V X E 



B = 



C2 



1 q{l — (f') V X f 



In spherical polar coordinates 

B = Brf + Bee + B^^ = B^^ 



since B ^ v x and f^O^cj) are unit vectors. 

The magnetic field lines form concentric rings about and there is sym- 
metry about the plane 0 — 7^/2, 

In the non-relativistic case of <C 1, /5 ^ 0, and then 

g 1 q V X r 

which is the Biot-Savart law. 

Applied to a continuous line current /: 

^ 1 I dl X r 
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V 



The constant l/(6:oC^) is normally written /xq - the magnetic permeability of 
free space. 



5.4 Field equations for an elementary point charge mov 
ing with uniform velocity 

(1) Equation for the total electric flux 
Consider the total electric flux 

= J E^hdS , 



We will use the axial symmetry and break sphere up into rings lying between 6 
and e + de. 



§ 
I 
1 
I 




r \ 



V 
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Since E h, the flux d'^E is 



= E-hdS = EdS . 



However 



dS = 27r(r sin 9) rd9 = sin BdO . 



Hence 



d^^E = 



1 



q{l-f3^) sin Ode 



3 
2 



27vr^ sin OdO 



2eo 



1 - /32 sin^ O) 



Thus 



/ 



TV 



sin Ode 



E 



^=0 (l - 13^ sin^ e) 



To calculate the integral, put cos ^ = a;, so that sin 6d9 = —dx, and then 



/ = 



-/ 



dx 



If dx 



dx 



where a = y^l — Performing the integration, we obtain 



/= - 



1 



X 



a 



Thus, including the Umits of the integral 



-1 



/ = 



-/ 



dx 



( (1-^^2 + ^2^2)3/2 



1 -1 



+ 



1 



1 
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and then 



9(1 - P^) 2 

2eo 1 - P- 



So 



Hence, the electric field produced by a moving charge satisfies the Gauss's law. 
(2) Magnetic flux through a closed surface 

Consider the magnetic flux through a closed surface 

= j B-hdS . 



If we choose a sphere centered on q to calculate the integral, we find that B 
is perpendicular to n everywhere, and then = 0. 
Similarly, V • 5 = 0 always. It is easy to see that 



\(e xv-v-V X E) = Q 



since V x v = 0 as i; is constant, and v is perpendicular to V x 
(3) Spatial E derivative related to temporal B derivative. 



We shall show that for a point charge under uniform velocity 



V X E = - 



dB 

'at 



In spherical polar coordinates 



V X £ = 



r sin^ 



d {sin OE^) _ dEe 



+ - 



e_ 

r 



1 dEr d{rE^) 



said dcj) dr 
d (rEe) dEr 



dr 



06 
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Since 



Eq ^ — 0 ^ 



then 



VxE = — —^0 
r ou 



with 



p 1 3(1 - 



K 

(i - /?2 sin^ e) 



3 

2 



Next 



^ = r-^W(l-/52sin2^)"^(-2/?2sin^cos< 

Si^y^^ sin 26 



2 [I -13'' sm^e^'^ 



Hence 



V X E = , ^ r </> • 

To calculate dB/dt^ we use the theorem of partial derivatives, 

then from the maximum change of y 

dy = ^dx+^dt = 0 , 
ax at 

we obtain 

dy dy dx 
dt dx dt 
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Thus 



dB _ dB 
dt dx 

Alternatively, to see this "physically", remember that the field pattern 
moves with constant velocity v. Let a stationary observer measure the change 
in the field 5 in a time interval dt. This change is the same as he would 
observed at a fixed time by moving a distance dx — —vdt^ i.e. 

dB in dt = dB in dx — —vdt . 



Hence 



and then 



Now 



and 



Since 



_ as 

dx vdt ' 



dB dB 

— = —V — 

dt dx 



- 1 q(l-f3^)vsme ^ 

B = : (p 



2 



dB _ idB^dr_ ^ 
dx dx dr dx 



sin 0 = = — 



we can write 5^ as 

Ka Ka 

Ba. — 
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where 



Next 



Hence 



and then 



2\l/2 



dx dx 

1 / 9 9 9\-l/2 X 

= - [x^ + z^] 2x = - = cos6 

2 v / r 

SKa r 



(r2 - /J2a2)5/2 



55,^ 1 3g(l - t;r2 sin 9 cos ^ 



5rc 



1 3g(l-y^^) t;sin2^ 
4^eo 2cV3(l-/?2sin2^)^ 



= —V (p = 



dt dx - 47reo 2c^r^ {l - (3^ sin' of 

1 3q{l- 13^)13^ sin 2^ 



47r£o 2r3 (l - ^2 sin2 ^) 



Comparing with V x we see that 

dB 



V X E = - 



dt 
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5.4.1 Electromagnetic Induction 



It has been known since about 1831 when Faraday first waved a magnet 
near an electric circuit and played with transformers that when the magnetic 
flux through a circuit changes, an electromotive force (emf) S appears in it. 
Faraday gave the rule 



The nature of this "phenomenon" is, however often misinterpreted. 

Consider a point charge moving near a closed circuit as shown in the 
figure. 




Because of the 6 dependence of the electric field on side (1) is larger 
than that on side (2). Thus, there is a net driving force round the circuit. 
Calculate the resulting electromotive force in the circuit, which is equal 

to the work done on a charge in the circuit 

r F ' dl r ^ - - 

£ = Wq = J = j{E + uxB)'dl , 

where u is the velocity of the charge g. 

Now since uxB \s perpendicular to both u and 5, it is also perpendicular 

to rfZ, and then 
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Hence 



J E-dl = JvxE-ndS= JvxE-dS 

-I " 



However, it is obvious in this example that the changing magnetic flux is not 
the CAUSE of the emf. The changing magnetic field and the electric field 
have a common CAUSE through the charge q. 

We can conclude that: Electric and magnetic fields do not produce 
each other - they are both due to electric charges. 

It is often thought in the textbooks, however that e.g. in a transformer 
the changing produces 8, It happens because the flux cutting rule is an 
extremely powerful one for calculating the integrated electric field of electric 
currents. 

— * — * 

The Faraday's rule £ = —d^M/dt^ which arises from V x ^ = —dB/dt 

should not be thought of as a casual relationship. What it means is that if 

a charge moving with a constant velocity produces a time varying magnetic 

field then that charge also produces a spatially varying electric field. 

(4) Relation between spatial variation of B and temporal varia- 
tion of E. 



Since 



B = 



V X E 



V x(vxE) = (^E'V)v-{v'V)E + v{y 'E)-E{V ^v) 
and V is constant, we obtain 



and then 



V X 5 = X (vx E 
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Now 



as 



and then 



dE_dE 
dx dt 



Vx — V Vy — Vz — yj and — = —v 



For points in free space 



^ ^ ^ dE 

at 



dt dx 



In summary 

Maxwell's equations for a point charge moving with uniform velocity are 

V X E 

V X 5 

These equations arise from the necessity for the correct relativistic trans- 
formations between frames in uniform relative motion. If the postulates of 
relativity are correct and Coulomb's law gives the field of a stationary charge, 
these equations follow, and the force on a charge is 
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Questions 

(1) Show that magnetic and electric fields of a charge moving with a 
constant velocity v are related by 




(2) Show that the magnetic field lines produced by a charge moving with 
a constant velocity v form concentric rings about v, 

(3) Show that the electric field produced by a moving charge satisfies the 
Gauss's law. 

(4) Explain the statement: Electric and magnetic fields do not produce 
each other - they are both due to electric charges. 
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6 Energy and Momentum in the Electromag- 
netic Field 

Since the electromagnetic field propagates at finite speed it must be assigned 
energy and momentum or else we must say that energy and momentum are 
not conserved in the finite time of propagation. 



The field takes a time At = r /do travel from the source charge to the 



6.1 The Energy Conservation Theorem — Poyntings' 
Theorem 

Energy may be transported through space by means of electromagnetic waves. 
We expect energy flow in the direction of propagation of the wave: E x B. 



We will show that the power flow across an element of area dS is given 
by c^SqE X B • dS, To derive this relation consider a volume V bounded by 




detector (another charge). 
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a closed surface 5, and containing a charged material. 
Consider the expression 

V'{ExB) . (12) 
If we employ the vector identity 

and use the Maxwell equations III and IV: 

^ ^ dB ^ ^ r 9E 

V X E = -— and V x 5 = fioJ + ^o/^o^ , 

we can write Eq. (12) as 

V • (i X 5) = -5 . ^ - /io^ • /- SofioE • ^ , 



or 



fiQ dt \2 2 fiQ 

On Ihs we put l//io = ^o^^ ^iid integrate the equation over some closed 
surface S enclosing a volume V, Then, we obtain 

Soc^V '{Ex B)dV ^ - i E' JdV - ^ / ( -^o^^ + ^— ) dV . (13) 

J Ot J \2, 2 /in / 

V V V ^ 

Now we apply Gauss' theorem to the Ihs of the above equation, and find 

SqC^ {E X B) ' hdS Energy flux 

s 

— — j) E * JdV Rate of doing work by field on the current 

V 



V 



This interpretation of the terms in this equation is suggested by some 
special cases: 
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6, 1.1 Ohmic Heating 



Prom the circuit theory, we know that in a resistive medium Ohm's Law is: 
V = IR and R — TZi/A^ where TZ is the resistivity of the medium. 




Thus, / = V/R = {VA)/{m). Since E = V/i {E = -V V), we get 

I/A = ^E . 

Defining the conductivity a — Ohm's Law becomes: 

I/A = J = aE 

The rate of conversion of electromagnetic field energy in to heat is then: 

y2 ^2 f 
A 



R 7^4 



where V = £A is the volume of the resistive medium. 

Hence 

^=aE^=EJ = E'J, 
Thus, E ' J is the rate of heating per unit volume in this case. 
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6,1,2 Electrostatic Field Energy Density 



It is well known that the work required to charge a capacitor of a capacitance 
C to a voltage V is ly = ^CV'^. According to the field theory of electro- 
magnetism this work done corresponds to conversion from some other form 
of energy into electrostatic field energy 

S = W = ICV^ and C = , 
2 d ' 

and the electric field in the capacitor is given by V = Ed. Thus: 

where V = Ad is the volume of the capacitor. 

Hence 

e 1 ^2 

6.1.3 Magnet ost at ic Field Energy Density 



It is well known that the work required to energize an inductor of L [henrys 
to a current I is S = \LiP and that the magnetic field within a long solenoid 
of self-inductance L — fion^Ai is 5 = fioul. According to the field theory 
of electromagnetism this work done corresponds to transformation of some 
other form of energy to magnetic field energy. 

Thus, for a solenoid of length £ and area of cross section A we have energy: 

e = w = Ili^ 

2 

2 fxo 

or 

£ _ IB^ 
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All our observations are in accord with supposing that the expressions 
for energy density found in electrostatic and magnetostatic situations are 
applicable in general. 

The expression 

N^Soc^E X B) 

is referred to as the Poynting vector. 

It represents the energy flux in the electromagnetic field, i.e. the energy 
flow per unit area (measured normal to the flow) per unit time. 
The units of N are thus watts/square meter. 

The energy flow equation can be converted into the form of a differential 
continuity equation or energy conservation law. From Eq. (13), we have 

at 

where 

is the energy density of the EM field. 



The physical meaning of the differential continuity equation is that the 
time rate of change of electromagnetic energy within a certain volume, plus 
the energy flowing out through the boundary surface of the volume is equal 
to the negative of the total work done by the fields on the source inside the 
volume. Thus. J • ^ is a conversion of electromagnetic energy into heat 
energy. 



Example 

An Earth satellite transmits a power P = 5 W. What are the tj^ical E 
and B in the radiation field at a distance of 1000 km from the transmitter? 
The mean flux (Poynting vector magnitude) at a distance r meters is: 

N = = = 4 X IQ-^" [Wm-'] . 

47rr2 47rl0i2 ^ ^ 
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Since 



E ±B and jB = 



E 



N = SqcE' 



Hence 



4 X 10"^^ = 8.85 X 10"'^ X (3 X lO'^) E 



-12 



Ec^h2x 10"^ Vm 



-1 



Thus, an antenna consisting of a straight piece of wire 10 cm long would have 
a voltage V :^ :^ 1.23 x 10~^ x 0.1 or :^ 12 /xV induced in it if it were 
placed parallel to E, 
Then 



B = E/c= 1.2 x 10"7(3 x 10^) ^ 4 X 10"'^ [T] . 

It is a very weak field. For, example magnetic field produced by a hair 

dryer is of order 10~^ T. Magnetic field produced by a large laboratory elec- 
tromagnet is about 5 T. 
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6.1.4 No Fluxes from Static Fields 

Consider a source of electrostatic field E and magnetostatic field 5. If it is 

arranged so that £^ _L 5 it is sometimes asked whether we should expect to 
see an energy flux of SqC^E x 5. 




"e 


B 







source_of 
static"! 



The question is easily answered in the negative by noting that the elec- 
trostatic field is E^ 0 and the magnetostatic field is 0, B. Each of these has 
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E X B = 0, The E^ B obtained by mentally combining the electrostatic E 
from one source and the magnetostatic B does not constitute an electromag- 
netic field. 

Electromagnetic fields do not interact with each other 

The fields interfere (as the term is used in physical optics). Beyond the 
interference region the two plane waves are the same as they would have been 
had there been no interference. If there are charges in the region of inter- 
ference then this is a different matter. The charges may absorb and re-emit 

radiation and so permanently modify the radiation fields. 



6.2 Phase Relationships in Harmonic Waves 

Only the in-phase components of E and B contribute to net energy flow 
averaged over a whole cycle of the radiation. 

li E — Eq cos(u;t)i and B — Bq cos(u;t)j then: 

N = SqC^E X B = SqC^EqBq cos^ {(jjt) k . 



N — —SqC? EqBq — SqC^ Ej.^sBj.^s 5 



Since, cos^(cc;t) = I, we obtain 

^ 1 
2 

where Ej^^ns = Eo/\/2 and B^ms — 

If however E — Eq cos{u;t)i and B — Bq sin(ct;t)j then: 

N — SqC^E X B — SqC^EoBo cos(cc7t) sin(ct;t) k . 



Since, cos(a;t) sin(cc;t) = 0, we have TV = 0. 



6.3 Momentum Flux 

To obtain an expression for the momentum carried by the electromagnetic 
field, we may employ the relativistic energy-momentum relationship 

c2 2 2 I 2 4 
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Since for electromagnetic radiation mo = 0, we obtain 

S 

p= - . 
c 

Thus the momentum flux of the electromagnetic field is: 

M = = SqcE X B . 

c 

6.4 Electromagnetic Energy Flow 

We now consider two simple examples illustrating applications of the Poynt- 
ing vector to circuit theory to show how the field theory provides an alter- 
native way of viewing some circuit problems. 



6.4.1 Energy Flow into a Resistive Wire 

Consider a wire (resistor) of length carrying a current /. Let 1/ is a po- 
tential difference applied along a wire. 

Circuit theory calculation 




According to circuit theory the power dissipated in the wire is P = VI^ 
where / is the current flowing through the wire. Thus, according to the cir- 
cuit theory, energy flows along the wire. 



Field theory calculation 



Lets look at the same problem from the point of view of the field theory. 

According to electromagnetic field theory the flow of energy should be de- 
scribed by the Poynting vector N = Sqc^E x B. 
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N 

A 

t r 



Calculate the electric and magnetic field produced by the current. The 
electric field propagates along the wire, and is given by 

V 

E = -VV = jz. 

From the Amperes line integral theorem, we find the magnetic field around 
the wire 

27va 

Hence 

N = eoc'E xB = Soc'^^ f 

I zTra 

since z x (j) — f. 

Since Sofioc^ = 1, we finally obtain 

N — — — r . 

Thus, the field theory predicts that energy flows into the wire from the air 
not along the wire. The energy is in the fields, the wire provides boundary 
conditions and guides the fields. 

If the energy is in the fields, it means that the electromagnetic energy 
goes out of a battery into the air, and then goes into the wire from the air. 
This is exactly the case we will show in the following example. 
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The total rate at which field energy flows in to the surface is given by 



^AT • dS 

s 



N-rdS= iNdS 



/ 



= iV 0 dSside = 



VI 

27da 



27r£a = VI , 



[N ± dS on the ends of the cylinder), which is in agreement with the result 
of the circuit theory. This result demonstrates quantitatively that the power 
which heats the resistor enters through the sides not through the wires. 

6.4.2 Energy Flow out of Battery 

In the above example, we have shown that the energy enters the resistor from 
the air. Then, a question arises: If the energy enters the resistor from the 
air, how does the energy get out to the air from a source of energy (battery)? 

Consider a battery, which provides energy to the resistor. 




Inside the battery J and E are in opposite directions. The magnetic field 

circulates around the battery, so we see that the Poynting vector points 
out into the air, not along the wire. 
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6A.3 Propagation of Energy along a Wire 

Above we showed that field energy flows in to a wire so that it can be dis- 
sipated as heat. It is well known that electromagnetic energy can be trans- 
mitted along a wire with very little loss. This is a different situation. If we 




consider the case of a perfect conductor we find that there will be a current 
wave along the wire with surface charges induced producing a radial electric 
field. Then E x B is parallel to the wire and the field description of energy 
transmission is that it is transmitted in the space around the wire. In the 
space around the wire E and B are in phase and N is always in the same 
direction. Within a perfect conductor we will show that E and B are 7r/2 
out of phase so the mean N averaged over a cycle is 0. So there is no net 
energy transmission within the perfectly conducting wire. 



Exercise in class: Energizing of a capacitor 

Consider the energizing of a plane parallel capacitor with circular plates. 

Show that circuit and field calculations agree as to the rate of energizing of 
the capacitor, i.e. Pc = Pf where: 

Pc = VI — rate of doing work (by current / and voltage V between the 
plates) in charging the capacitor according to circuit theory. 
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Pf = J Sqc^E X B ' dS — rate of energy flow into the surface of the ca- 

s 

pacitor according to field theory. 

From which direction does the energy flow in to the capacitor according 
to field theory? 

We will keep the calculation simple assuming the plates of the capacitor 
to be uniformly charged. Under what conditions is this assumption likely to 
be true? 



Questions: 

(1) Using the Maxwell's equations derive the continuity equation for the 
Poynting vector. 

(2) Show, using the field theory calculation, that the power dissipated 
along a resistive wire is P = VI ^ the same predicted by the circuit theory. 
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7 General Solution of the MaxwelPs Equa 
tions 



The electromagnetic fields of charges in uniform motion are effectively bound 
to the charges. The fields of accelerated charges on the other hand can prop- 
agate as electromagnetic (EM) waves at the speed c and can have a life of 
their own (until absorbed by some other charges). 

Consider a simple model showing why only accelerating charges can emit 
electromagnetic waves. 

Charges in uniform motion produce a constant magnetic field, i.e. when 
V — const then B — const. EM wave is produced when B 7^ const. From the 
Faraday induction law 

8{t) — ——f- / const 5 

where ^ = S(t) • 5 is the magnetic flux. 

B 7^ const when v / const (acceleration). 



7.1 Difficulty of the Direct Solution of Maxwell's Equa- 
tions with Time Varying Sources 

The existence of an electromagnetic field implies the presence of currents and 

charges. If the currents and charges are known, we can find the fields solving 
the Maxwell's equations. 

Consider the Maxwell's equations 

1. V • ^ = 
II. V • 5 = 

III. V X ^ 

IV. V X 5 
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In general case, E{f^t) and B{f^t) depend on (r, t), the charge and current 
densities also depend on (rjt). 

First, we will try to separate the Maxwell's equations into an equation for E 
alone or B alone. 

Assuming in the usual way that space and time operators commute, we act 
with on III and Vx on IV, and obtain 

1 d 

— — V X E- V X (V X 5) = -/xqVxJ. 
at 

— * 

Eliminating E by subtraction 

1 d'^ 

V X (V X 5) + = /ioV X J 

and using the vector identity for double x product, and II 

V X (V X S) = - + V(V • B) = -V^5 , 

we obtain 

V'^-^^^=-/^oVx/. (14) 
Similarly, elimination of B gives 

Equations (14) and (15) are in the form of coupled wave equations known 

as inhomogeneous Helmholtz equations. We see that the current density J 
enters into these equations in a relatively complicated way, and for this rea- 
son these equations and are not readily soluble in general. 

In the absence of currents and charges, J = 0, p = 0, and then the above 
equations describe a free EM field, and can be solved separately. 
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Solution of the wave equations in the absence of currents and charges is given 
in a form of plane transverse waves 

U^Y^Uk e-<"^'-^"1 , (16) 
k 

where U = {E^ 5), k = cjfc/c, and Uk is the amplitude of the E or B wave 
propagating in the k direction. 

The general solution of the wave equations, in the presence of space and 
time varying currents and charges, is more readily attained via the electro- 
magnetic potentials. 



7.2 Scalar and Vector Potentials 



Generally, we do not find fields E and B directly by integration of Eqs. (14) 
and (15). We rather first compute scalar and vector potentials from which 
the fields may be found. 

We will illustrate here the advantage of working with the scalar and vec- 
tor potentials. 

Scalar potential is a quantity from which a field can be derived by a pro- 
cess of differentiation, e.g. in electrostatics 

E = - V$ , 
where $ is the electrostatic potential. 

Introduce the vector potential A defined such that the Maxwell's equa- 
tion II remains unchanged. Since V • 5 = 0, we can always write 

5 = V X A . 

Substitute this relation to the Maxwell's equation III, and since V x V$ = 0, 

where $ is an arbitrary scalar function, we find 

E^-^A-V^ . (17) 
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In the static limit of dA/dt — 0, the scalar function $ reduces to the 
familiar electrostatic potential. 

The above equation shows that the electric field depends on the specific 
choice of the potentials. 

We can define new potentials without changing E and B 

= $ - I^A . (18) 

Proof: 



at 

d 



+ va) - V ($- 



dt 
d ^ 

— A- V$ = E. 
dt 




5' = V X 1' = V X 1+ V X (VA) 
= V X A^B . 

as required. 

The transformation (18) is called a gauge transformation^ and the invari- 
ance of the fields under such transformations is called gauge invariance. 

The definition B = V x A does not completely define A despite the fact 
that B is completely defined. The vector potential A is arbitrary to the ex- 
tent that the gradient of some scalar function can be added. Thus, infinite 
set of possible potentials $ corresponds to an infinite set of possible vector 
potentials. However, the Maxwell's equations should be independent of the 
specific choice of potentials. 
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7.2.1 Lorentz Gauge 



Recall the Helmholtz Theorem which says that any vector field can be 
written as a sum two terms 

- 1„/*V-F,_ 1„ f V xF 

F — — 



V / dV + —V X / dV 

Jv r Air Jv r 



= Fi + Ft, 

— * — * 

where Fi is called the longitudinal part of the field and has V x = 0, 

while Ft is called the transverse part and has V • Ft = 0. 

We see that V • F and V x F together determine F but neither do alone. 
Thus, if we define V • A, we complete the definition of A. This is called 
"choosing the gauge of the potential". The above is an excellent illustration 
of the power of the Helmholtz theorem. This theorem enables us to recog- 
nize basic common properties of vector fields independent of their individual 
physical properties. 

Derive differential equations for A and 

From the Maxwell's equation I and (17), we find 

- <9 - 
dt 

The electric field E will satisfy the Maxwell's equation I when 

From the Maxwell's equation IV and 5 = V x ^, we have 

Vx(Vxi)=;.„/+i|(-|i-V$) , 

which can be written as 

-VM-+V(V.l) = ;.o/-i^A--iv|*. 
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Hence 



We have already shown that 

-|(V-A)-V2$ = p/£o. (19) 

The freedom of choosing A and $ means that we can choose a set of 
potentials to satisfy the condition 

ot 

This is called the Lorentz gauge and defines V • A. This equation is some- 
times called the Lorentz equation. 

Under the Lorentz gauge, the Maxwell equations reduce to two uncoupled 
wave equations 

- 1 5^ - 

and 

7.2.2 Coulomb Gauge 

Another useful gauge of the potentials is the Coulomb gauge or transverse 
gauge 

The origin of the name "Coulomb gauge" is in equation (19) that under the 
condition V • A = 0 reduces to the Poisson equation 
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that determines the Coulomb potential due to the charge density /?. 

Before we proceed further, we will show that the solution of the Poisson 
equation is of the form 



It can be proved in the following way. 
From the Coulomb's law 

Attsq J 

and using the relation 

-v-= " 



2 ' 



we can write the electric field as 



E = -y^V / -dV = -V$ , 



Attso J r 

where 



— f^dV , (20) 
veo J r 



Aire 

Since the electric field satisfies the Maxwell's equation I, we find 
as required. 

Now, we can find the wave equation for A under the Coulomb gauge. 

Under the Coulomb gauge, the vector potential satisfies the inhomoge- 
neous wave equation 
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The term involving the scalar potential is called "longitudinal" as it has 
vanishing Vx. This suggests that it may cancel the longitudinal part of the 
current density J. According to the Helmholtz Theorem, the current density 
can be written as 

J = -—V/ dV + --Vx / dV 

Att Jv r Att Jv r 

= Ji + Jt^ 
Using the continuity equation 

and the solution of the Poisson equation, we find the longitudinal part of the 
current density 



- — V / dV = -—V / — ^ 

47r Jv r Att Jv r 

I d _ f p d 



An 



at Jv r ot 



Then 



Hence, the inhomogeneous term in the wave equation for A can be ex- 
pressed entirely in terms of the transverse current. 



1^' 

This is the origin of the name "transverse gauge". 



7.3 Solution of the Inhomogeneous Wave Equations 

We have shown that the Maxwell equations can be reduced to two indepen- 
dent wave equations for the potentials A and In fact, we have four scalar 
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equations for {A^, Ay, A^,^). Each of these equations has the same form. 
We will illustrate the solution on one of the equations, $: 

A general solution may be found by considering two limiting cases: 
(a) Electrostatic limit d/dt = 0 

In this limit the wave equation for $ reduces to the Poisson equation 
whose the solution is 

Hr) = -^f(^dV'. 



(b) Source free limit p = 0 

In this case, the wave equation reduces to the homogeneous equation 

V $ = 0 . 

This equation has a spherically symmetric solution of the form 



/(/ - r/c) 



where /(t — r/c) is an arbitrary function of the retarded time t — r/c. The 
retardation r/c is equal to the time needed for the electromagnetic wave to 
pass the distance from the source to a given point in space. 

Proof: 



If there are no boundary surfaces, the potential can depend only on r, 

and must in fact be spherically symmetric. Thus, in spherical coordinates 
only the radial part of the Laplacian will contribute to the wave equation 



1 d 



r 



5$ 
dr 
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Hence, we have 



df df dt, 



dr dtr^ dr 



-V, 

c 



where tj. —t — r /c and /' = df /dtj.. 



Thus 



= — r 




c c \ c 



1 



rc 



where /" = d^f/dti. 

Since d/dt = d/dtr, we have 

1 

as required. 



1 



We can construct a general solution of the wave equation by noting that 
it must represent a spherical wave outside the source and reduce to the 
appropriate static limit. This solution is 

p{t - r/c) 



dV , 



where r is the distance coordinate from the source (from the charge pdV) at 
the time when the potential wave left it. This exhibits the causal behavior 
associated with the wave disturbance. The argument of p shows that an 
effect observed at the point r at time t is caused by the action of the source 
a distant r away at an earlier or retarded time — t — r / c. The time r/c is 
the time of propagation of the disturbance from the source to the point r. 
Thus, the Maxwell's equations satisfy the causality principle. 
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7.4 Alternative Solution: Green Functions Method 



The wave equations all have the basic structure 



V^$(r, t) 



c 



1 a2$(r,t) 



47r/(r,t) , 



where f{r^t) is a known (source distribution) function. 

To solve this equation, we will introduce the Green function of the equa- 
tion and solve it as an inhomogeneous Helmholtz equation. 

Suppose that $(r, t) and /(r, t) have the Fourier integrals 



When we insert it into the wave equation, we find that the Fourier trans- 



where k = uj/c is the wave number. 

The advantage of working in Fourier components is to remove the deriva- 
tive over time. 

Green function 

For a unit point source the potential satisfies the Poisson equation 



The function 1/r = G{r) is called a Green function of the above differen- 
tial equation. 

In analogy, we can define the Green function of the wave equation 




form $(r, (x;) satisfies the inhomogeneous Helmholtz wave equation 





A7r6{r) . 
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The Fourier transform gives 

(V^ + k^) Gk = -47r5(r)e^^*« . 

where Gk is the Fourier transform of the Green function G{r^t — to), which 
we are trying to find. 

If there are no boundary surfaces, the Green function depends only on r, 

and then the Laplacian operator in spherical coordinates depends only on r 
giving 

1 d? 



r dr 



2 



Everywhere except r = 0, the function rG^e ^"^^^ satisfies the homoge- 
neous equation 

^ {rGke-'^'') + e{rGke-'^'') = 0 , 
whose the solution is 

rGke-'""^' = Ae^^" + ^e"^^^ . 
Thus, the general solution for the Green function is 



Gk — 



Using the inverse Fourier transform, we find 



I POO p,-^""" 

^ ' ^ 27r J-oc r 



where r — t — t^. 

The integral 



1 



27r ^-oo 

is the delta function 5(r =F r/c). Thus 



J — oo 



G{r,r) = ^S{r =F r/c) 
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The Green function is a casual response function, and has the same property 
as the scalar potential of a point source. 



In summary 

The general (retarded) solutions of the Maxwell's equations are 

- - - dA 

B = V xA , E = -V$ - — , (21) 

at 

with 

In practice, we calculate the scalar and vector potentials from Eqs. (22) 
and (23), and then find the electric and magnetic fields from Eqs. (21). 
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Questions 

(1) Find a differential equation for the vector potential A in the Lorentz 
gauge. 

(2) Explain, why the Coulomb gauge is often called "Transverse gauge". 

(3) Prove that the homogeneous wave equation 

1 d^^ ^ 
V $ = 0 

has the solution of the form 

, f(t - r/c) 
r 

where f{t — r/c) is an arbitrary function of the retarded time t — r/c. 
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8 Solution of Laplace's Equation and Bound 
ary Value Problem 



In the previous lecture, we have shown that electric and magnetic fields 
are readily calculated with a help of the scalar and vector potentials. In 
this lecture, we will illustrate applications of the scalar potential to physical 
problems involving bounded fields. 

There is a class of problems in electromagnetism in which a field can be 
derived without involvement of the complete set of the Maxwell's equations 
as the gradient of a scalar potential which satisfies Laplace's equation 

= 0 . 

The condition for this to happen is that 

V*F = 0 and VxF = 0. (24) 

Since, V x F = 0 is the condition for F = V^, where $ is an arbitrary scalar 

function, then V • F = 0 means that 

V • (V$) = = 0 . 

Thus, the scalar potential $ contains all the necessary information to com- 
pletely specify the field of the properties (24). 

Examples 

1. Electrostatic problems involving Laplace's equation 

Since in general V • F = p/e^^ and VxF = —^B we see that the require- 
ment for Laplace's equation to be relevant is that /? = 0 and d/dt = 0, i.e. a 
source-free region and static conditions. Of course there must be a source of 
charges somewhere or there would be no field anywhere. The typical situa- 
tion where solution of Laplace's equation is relevant in electrostatic is where 
we have source-free non-conducting regions between statically charged con- 
ductors. 
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2. Magnet ost at ic problems involving Laplace's equation 

Since in general V • 5 = 0 and V x B = fJ^oJ + we see that the 

requirement for Laplace's equation to be relevant is that J — 0 and d/dt = 0, 
i.e. a source- free region and static conditions. Again, there must be currents 
somewhere or there would be no fields anywhere. The typical situation is 
to be calculating the magnetic field in the non-conducting region between 
constant currents. 



8.1 Uniqueness of the Solution of Laplace's Equation 

As we will see below, general solutions of Laplace's equation are in terms 
of some constants which are usually found from boundary conditions for a 
given problem. A question arises: What boundary conditions are ap- 
propriate for the Laplace equation to ensure that a unique and 
well-behaved (physically reasonable) solution v^ill exist inside the 
bounded region? Our experience leads to believe that specification of the 
potential on a closed surface defines a unique potential problem. This is 
called Dirichlet theorem or Dirichlet boundary conditions. 

8.1.1 Dirichlet theorem 

Consider a volume V completely bounded by a closed surface S, Within S 

there is a potential $ satisfying = 0. The Dirichlet theorem says that 
the value of $ is uniquely determined by the variation of $ over S, 

Proof: 

Suppose, to the contrary, that there exist two solutions $i and $2 satis- 
fying the same boundary condition, i.e. V^$i = 0 and V^$2 = 0 within 5, 

but $1 = $2 on S. 
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Let U = $1 — $2 is the difference between the solutions. Since $i and $2 
are known to be sokitions of the Laplace equation, then from the linearity 
of the operator V'^U = 0, i.e. U is also a solution of the Laplace equation. 



We will prove that U = 0 inside the volume. To show this, we introduce 
a vector F = UVU, Then using the property 

V • F = V • UVU = C/V • (VU) + VU'VU 

and the Gauss' Divergence Theorem, we get 

/ V • FdV = / C/V • VUdV + / VC/ • VUdV 
Jv Jv Jv 



= UV^UdV + {VUf dV 



= / UVU'dS , 
Js 

Now RHS= 0 because C/ = 0 over 
Also, the integral 



Jv 



because $1 and $2 both satisfy the Laplace equation throughout V , 

Hence 



Since the integral from a positive function is always positive, VU must be 
zero for the integral to be zero. Thus VU = 0 and consequently, inside U 
is constant. Since U = 0 on S", so that inside ]/, = $2- 



8.2 Solutions of Laplace's Equation 



There are different methods of solving the Laplace equation 
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• Method of Images 

• Green functions method 

• Variational method 

• Method of lattices 

• Numerical Monte- Carlo simulations method 

• Method of separation of variables 

• Solution in spherical coordinates 

We will illustrate last two methods which can be applied to a large class 
of problems in electromagnetism. The other methods can be applied to spe- 
cific problems. For these methods it is necessary that the boundaries over 
which the potential is specified coincide with the constant bounding surfaces. 



8.2.1 Method of Separation of Variables 



In cartesian coordinates the Laplace equation for the scalar potential can be 
written as 

dx^ dy^ dz^ 

Since re, z are independent variables, the solution of the Laplace equation 

is of the form 

^x,y,z)=X{x)Y{y)Z{z) . 

Substituting this into the Laplace equation and dividing both sides by XYZ^ 
we obtain 

1 (PX 1 (PY 1 dPZ _ 
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This equation can be separated into three independent equations. To show 
this, we write this equation as 



Both sides of the above equation depend on different (independent) variables, 
thus are equal to a constant, say —a^: 



1 d^X 

1 d'^Y 1 d'^Z 



= -q;2 



= -a' 



Y dy"^ Z dz^ 
The second equation can be written as 

= a. — . 

Y dy"^ Z dz"^ 

Again, both sides depend on different variables, thus are equal to a constant, 
say —P^: 



1 d^Y 



= -0- 



Y dy^ 

Zdz^ ^ 

Hence, after the separation of the variables, we get three independent ordi- 
nary differential equations 

1 2 



1 rf^y 

Y~d^ 



+ = 0, 

1 d'^Z 



The solutions of these equations depend on whether a'^ and are positive 
or negative. If we choose o:^ and to be positive, the solutions of the 
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differential equations are 

X{x) = ^ (AfeC^"^ + i^fcc-^'-) 



P 

The solutions can also be written in the form 

-^(-^0 — [^k sin(ax) + 5fc cos(ari;)] 

k 

y{y) = E sin(/3y) + A cos(^y)] 

Z{z) = Yl si^h (^^a^ + P^z^ + Fp cosh j^y^c^^ + ^2^ 

The above solutions are in a general form, where the constants 

Bk^Cij Di^ Ep and Fp can be found from specific boundary condi- 
tions. 

Consider two examples: 

1. We have a solution with known boundary conditions, find the problem. 

2. We have a problem with specific boundary conditions, find the solution. 

Example 1. 

Consider the following two-dimensional solution of the Laplace equation 

$(2;, z) = X{x)Z{z) = Vo sin(ax) sinh(az) 

with the lower boundary ^rnin = 0 and the upper boundary ^max = Vq. 
In what circumstance is the above the solution? 

Consider z) in some limits. X = 0 for x = 0 or ax — tt, i.e. x — tt /a. 
The solution thus satisfies the boundary conditions along the "vertical 

lines for a — T^jh. 
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Since sinh(az) = 0 for z = 0, the boundary condition along the lower bound- 
ary is satisfied. 

For the solution to satisfy the upper boundary condition, the shape of 
the upper boundary must be such that 

Vo sm{ax) sinh(az) = Vb 
for all points z on the line, i.e. 




Since sin i^^j ^ 0 at the edges, it is greatest (= 1) at the center x = 6/2. 
Hence, sinh (^) must be equal to one at x = 6/2. This happens when 

arcsinh(l) ^ 0.885 

0.8856 ^ 
= 0.2826 

TT 

We usually have reverse problems to the above that we have a set of 
electrodes which constitute equipotential lines or surfaces, and need to find 
the appropriate solution of the Laplace equation. This is illustrated in the 
next example, where we will try to find potential inside a rectangular box 
whose three sides have potential equal to zero, and the remaining side has a 
potential Vq. 



TTZ 

T 

z — 
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Example 2. 



x=b 




This two-dimensional problem has a general solution 

z) = ^2 [^n sin(Q;2;) + cos{ax)] 

n 

X [En sinh(Qf2:) + Fn cosh(az)] . 

The boundary condition $ = 0 at = 0 can be satisfied by Bn = 0. 
The boundary condition $ = 0 at z = 0 can be satisfied by = 0. 
To have $ = 0 at or = 6 we must have ab = nvr. 
Hence the solution reduces to 

oo / .N 'nirz 



z) = ^ sin f — 7— j sinh 

n=l 



where Kn — A^En^ 

To find Kn we apply the remaining boundary condition $ = Vo at z = a 



00 . 



a) = Vo = X] sin ( — t~ ) sinh 

n=l \ 0 J 



b 



This is a Fourier series in x and in the usual way we use the orthogonality 
properties of sine functions to calculate Kn- 



L 



• / A\ ' / A\ JO for m 7^ n 



sm{m(l)) sm{n(f)) d(j) — 



TT for m = n 
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f r A\ t A\ f 0 for m 7^ n 

/ coslmo) cosino) do = < „ 

Jo \ V 7- / 7- I TT tor m = n 

I 



27r 

sin(m0) cos(n<^) d(j) — 0 for all m and n 



^0 

Multiplying both sides by sm{m7rx/b) and integrating over x — 
we get 

/Vosin(^)da. = f:i^„sinh(^) 

7X — -L 

/ m7rx\ f mTX\ 

Jo y-r) • 

All integrals on the RHS are equal to zero except for m = n 



riTTX 



— cos 



^ b 1 f n7Ta\ 2 /^'7r2:\ . 

iv„ sinh — — / sin — — ax 



b J lo nir \ b J Jo \ b J 

which we can write as 



n ^ M ' u ^^^^ 

[1 — cos(n7r)J = Kn sinn 



riTT 



1 r f2n7rx\ 
X / - 1 — cos — - — 
Jo 2 1 \ b J. 



2n7TX\ 



dx . 



The cos (^^y^) integrates to zero over the range 0^6, giving 

Vob ' u ^^^^^ b 
1 — cos(n7r) = Kn sinh 

riTT 

and finally 

2Vo 1 — cos(n7r) 



b 



n 

riTva 



^TT sinh^ ^ J 

If n is an even number then cos(n7r) = 1 and Kn = 0. 

If n is an odd number then cos(ri7r) = — 1 and 1 — cos (nvr) = 2. 
Hence 

Kn = — 7 >- for odd n 
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and then 



^ 4Vo sinh(^^) , ^nirx 

^[x, z) = T sm 

oddn sinhl- 



nira 



8.2.2 Solution of the Laplace Equation in Spherical Coordinates 



In this lecture, we continue the discussion of boundary-value problems and 
will illustrate solution of the Laplace equation for general problems of spher- 
ical symmetry. 

In spherical coordinates (r,^, <^), the Laplace equation can be written as 



= — 



1 d 



+ 



1 



5$ 
dr 
d 



sin 6 06 



sin 6 



5$ 



+ 



1 



sin 6 d(jP 



= 0 . 



Multiplying by r^, the Laplace equation can be written as a sum of two 

separate parts 

d ( .d^\ 1 d ( . m\ 1 ^ 



dr \ dr 



sin 6 dO 



sin 9 d(j)'^ 



The first part depends only on r, whereas the second part depends only 

on 6, (j). Thus, the solution is of the separable form 

$ = R{r)Y{e, (j)) . 

Hence, substituting $ = R{r)Y(6,(j)) and dividing by R{r)Y [6 , (j)) , we 
obtain 



1 d 



dR 



r 

Rdr \ dr 



1 

Y 



1 



sin 6 dO 



Of. JY 
sinp— — 

m 



+ 



1 d'^Y 



sin^ e 
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Both sides of the above equation depend on different variables, thus must 
be equal to the same constant, say —a\ 

r^— \ +aR = 0 



dr \ dr 
1 a / . dY\ 1 d'^Y 



sin QdQ\ 09 ) sin' 9 dcj^^ 

Thus, the Laplace equation splits into two independent differential equations. 
We will call them (A) and (B). 

(A) Consider the equation for Y , 

Multiplying both sides by sin^ 6^, we get 

d f dY\ d^Y 

This equation contains two separate parts, one dependent only on 9 and the 

other dependent only on (j). Therefore, the solution will be of the form 

Y{d,cj>) = X{d)^{cl)) . 

Hence, we get 

I , ^d ( , dX\ . . Irf^* 
— sm6^— - sme^— - - asm 6/ = --r^^ • 
X d9 \ d9 J ^ d(j)^ 

As before, both sides must be equal to a constant, say m^: 

1 ^d ( dX\ 2 
— sin 9—- sin 9—— — a sin 9 — m 
X d9 \ d9 J 

Id^^ _ 2 

(Al) First, we will solve the equation for which we can write as 
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It is the familiar differential equation for a harmonic motion. The solution 

of this equation is of simple exponent form: 



^(</)) = Aexp{im(j)) , 

where A is a constant. 

To determine the constant m note that in rotation, 0 and 0 + 27r corre- 
spond to the same position in space: ^(^) = + 27r), which is satisfied 
when 



exp{im(j)) = exp[im(0 + 27r)] 



This leads to 



exp(i27rm) 1 , 



that is satisfied when m = 0, ±1, ±2, .... 

Hence, the constant is not an arbitrary number, is an integer, 
(A2) Now, we will find X{e), 



Using the solution for the differential equation for X{6) can be writ- 
ten as 



1 / . ^dX 
smu 

sin ede\ de 



-\a + - 



X = 0 



Introducing a new variable z — cos 9, we can rewrite this equation as 



/ ^\ d^X dX I m 
(l - ^ ) - 2z^ I a + 



dz^ 



or 



dz 



(1 - 



dX 
dz 



dz 



- a + 



\ — Z' 



X = 0 



X = 0 . 



The above equation is known as the generalized Legendry differential equa- 
tion, and its solutions are the associated Legendry polynomials. For m = 0, 
the equation is called the ordinary Legendry differential equation whose so- 
lution is given by the Legendry polynomials. 
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Lets look into the solution procedure of the above equation. This will 

allow us to find a and X{9), 

We assume that the whole range of z (cos^), including the north and south 
poles (z = ±1), is in the region of interest. The desired solution should be 
single valued, finite, and continuous on the interval — 1 < z < 1 in order to 
represent a physical potential. 

The differential equation for X has poles at z = ±1. In order to find the 
solution of this equation, we first check what solution could be continuous 
near the poles. 

Lets check a possible solution near z = L Substituting re = 1 — z, then 

dx — —dz and in terms of x the equation takes a form 



_d_ 
dx 



x{2 — x) 



dX 
dx 



x{2 — x) 



X = 0 . 



We look for a solution in the trial form of power series in x 



oo 



X{X) ^x'Y^^n 
n=0 



Substituting this into the differential equation for X, we get 



2s^aQX^ + (5 + l){2sai — sa^ + 2ai)x^ + . . . 



- a + 



x(2 — x) 



(ao + aix + . . ,)x^ — 0 . 



Near x ^ Q^we can replace x{2 — x) by 2x^ and obtain 



2s^aQ \x^ + • • • 



= 0 



This equation is satisfied for all x only if the coefficients at x^^ x^^^^ . . . are 
zero. From this, we find that 



1 



m 



We take only s — +i 



as for 5 = — I 



m 



the solution for X{x) at x = 0 



would go to infinity. We require the solution to be finite at any point x. 
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Thus, the solution that is continuous near 2: = 0 is of the form 



00 



n=0 

or in terms of z 



oc 



n=0 



Using the same procedure, we can show that near the pole ^ = — 1, the 
continuous solution is 



00 



n=0 

Hence, we will try to find the solution in the form 

00 

n 



n=0 



Substituting this equation into the differential equation for X{z) and collect- 
ing all terms at the same powers of z^, we obtain 



J2 {{n + l)(n + 2)fe,+2 - n{n - 1)6, 

n 

- 2{\m\ + l)nbn - {a + \m\ + m^)6„| z"^ = 0 
Hence, we get a recurrence relation for the coefiicients 



(n + |m|)(n + \m\ + 1) + a 

[n + + 2) 



We have two separate solutions for even and odd n. For bo 7^ 0, we put 
bi = 0, and the solution is given in terms of even n. For 60 = 0, we put 

bi ^ 0, and the solution is given in terms of odd n. 

We cannot accept both the even and odd solutions at the same time, because 

in this case the solution X{z) would not be a single valued function. 
For example, for bo ^ 0, we have a = —\m\ — m^, but for bi 7^ 0, we have 
a = —2 — 3\m\ — m^. If we would accept both of the solutions at the same 
time, the potential would have two different values. 
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We check now whether the series is converting when n ^ oo which would 
ensure that the potential is finite. 

Since 6^+2 > ^n? the series diverges for 2: = ±1. Therefore, in order to get 
the potential finite everywhere in the space, we have to terminate the series 
at some n = no. In other words, we assume that 6no+i = . . . = 0. 

The series terminating at n = no indicates that 



(no + |ni|)(no + |m| + 1) + Of = 0 



Introducing 



I = no + 



m 



we see that / > |m|, and 



a = -l{l + 1 
Thus, the solution for X{z) is 



Xlm{z) 



1 - Z 



2\ 2 



h\m\ 



l-\m\ 



z 



n 



n 



where the sum is over even n when / — |m| is an even number, and over odd 



n when I — 



m 



is an odd number. 



First few solutions 



Xii{z) 



60 - boP^iz) 
biz = hiP^{z) 
boVT^ = boPl{z) 



where P^{z) = 1, P^{z) = z, Pl{z) = Vl - z^, 
polynomials of the order /. 

Useful examples [in terms oi 6 (z = cos 6)] : 



. are the associate Legendry 



P°(cos^) 
Pi°(cos^) 

Pi (cose) 



= 1 



cos^ 

sin^ 
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P^{cose) = ^ [3 cos(2^) + 1] 
P2^(cos^) = ^sin(2^) 
P|(cos^) = |[l-cos(2^)] 



Orthogonality of the Legendry polynomials is an important property: 



-1 

unless m = n and I — for which 

[P,"^(cose)]' dicosO) = ji^^^j; 
With the above notation, the solution for Y{9^(j)) is of the form 



I m 



(B) What left is to find the radial part R{r), 

With a = —/(/ + 1), the differential equation for R is of the form 

|(,.^f)-,. + :). = o. 

Dividing by r and introducing a new function U{r) — rR{r)^ we obtain 

d^U l{l + 1) 



C/ = 0 . 



Lets first check the asymptotic solution for r ^ 1. In this case we can ignore 

the second term in the differential equation, and find that the asymptotic 
equation has a solution U{r ^ 1) = Cr^ where (7 is a constant. 
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Following this asymptotic behavior, we will try the general solution of a 

form 

U{r) = . 

Substituting this into the differential equation, we obtain 



[s{s - 1) - l{l + 1) 



= 0 . 



This equation is satisfied for any r when 



5 = (/ + 1) 



or 



s — —I , 



Thus, the general solution is of the form 



C/(r) = Cir^+' + C2r-' , 



or 



Thus, solution of the Laplace equation in spherical polar coordinates is of 
the form: 




I m 

xAimPr{cose)e 



The solution can be written as 




I m 



X 



aim cos{m(l)) + bim sin(m<ji))] P^'^(cos 6)} . 



Example 1: Potential inside an uniformly charged sphere 

Consider an example of boundary- value problem with azimuthal symmetry: 
A conducting sphere of a radius a in an uniform electric field. 
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The conducting sphere is an equipotential volume (else there would be 
electric fields driving currents till it became so). Take the zero potential on 
the sphere. 

The boundary conditions to be satisfied are: 

1. The potential on the surface of the sphere V{a^ 6^ (j)) — 0. 

2. The potential at infinity is the uniform field potential (no effect of the 
sphere at infinity), so $ = —Er cos 9 at infinity. 

Since the applied potential is independent of the angle phi^ the induced 
potential will also be independent of 0. Thus, we can set m = 0 in the 
general solution and get 

The boundary condition at infinity is satisfied for all constants Cu = 0 except 
for / = 1 (remember Pi = cos6'). 

$(r, 6) = CiirPi(cos^) + ^^i ' 

As r ^ oo, the potential ^{r^9) CnVcosO — —Er cos 9. 
Therefore Cn — —E, 
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The other boundary condition to be satisfied is $(a, ^) = 0 on the surface 
of the sphere. 

$(a,6') = 0 = -EaPi{cos6) + 2^ — - . 

We can determine the coefficients C21 using the orthogonaUty properties of 
the Legendry polynomials. Multiplying the above equation by P]^[cos9) and 
integrating over cos 9^ we obtain 



The first term vanishes by orthogonality of the Legendry polynomials. 
All terms in the summation vanish except that for I — k. Thus 



0 




0 




Since the integral is nonzero, then C2k = 0 for A: 7^ 1. 



For k 



1: 




2 (/ + m)! 



2Z + 1 {l-m)\ 



which for / = 1 and m = 0 is equal to 2/3. Thus 



0 = -aE- + 



C21 2 



a? 3 



from which, we find C\ 



21 — 



Eo?, Hence 



9) = —Er cos 9 + Ea 
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The first term is just tiie potential of a uniform field E. The second term 
is the potential due to the induced surface charges or, equivalently, is the 
potential of the induced dipole moment p — AireoEa^ 

pcosO 



More useful exercises: Before you approach the tutorial problems on this 
subject (tutorial set number 6), try to solve the following problems 

1. Consider a two-dimensional region with boundaries at a; = 0,6 and 

z = 0, a, as shown in the figure 3. The boundary conditions are 



dz 



at z = 0 , 

0 at re = 0, 6 , 
Vo at z = a . 



Find the potential at any point inside the two-dimensional region. 




Figure 3: 

2. Consider a conducting sphere of radius The surface of the sphere is 

kept at a potential 

^R,e,(j)) = Vosin^sin<^ . 
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Using the above as a boundary condition, find the potential at any 
point inside the sphere. 



Questions 

(1) Explain, why in the two-dimensional case, we take from the general 
three-dimensional solution of the Laplace equations 

$ = X{x)Z{z) , 

or 

$ = Y{y)Z(z) , 



but not 



$ = X{x)Y{y) . 
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9 Electromagnetic Antennas: Hertzian Dipole 



In this lecture we will illustrate an application of the vector potential A in 
calculations of the electric and magnetic fields produced by a source system 
containing time varying charges and currents. 

Consider the retarded solutions of the Maxwell's equations 

^ ^ _ fjX 

B^V xA , E = -V$ - — , 

at 

with 

AtvSq J r 



Atveq^ J r 
The above solution holds for the Lorentz gauge in which 

at 

Assume that the charges and currents vary sinusoidally in time 

J{r,t) = /(r)e^"' . 

In this case 

V- A= — , 



C2 



which gives 

$ = V- A . 

iuj 

Thus, the scalar potential can be eliminated from the field equations 
leaving only the dependence on A. Hence, we can express both E and B in 
terms of the vector potential A alone. We have 

E=-V{V-A)-^ , B = VxA. 
tu> at 
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This result may seem rather strange at first, since normally we should expect 
to need both the scalar and vector potentials in order to completely determine 
the field. The explanation and in fact an another way of saying the same 
thing is that time varying charge must satisfy the continuity equation 

V ' J = -— = -tujp , 

so that 

P= - 



Then 



$(r,t) = / ^ ^ dV 



1 fV 'J 



-I 



dV 



1 V- / - dV^-—V'A 



Thus, specification of J alone is sufficient to completely determine all sources, 
and hence a solution for A in terms of J contains all the necessary informa- 
tion to completely specify the time- varying field. 



9.1 Field of an Element of Alternating Current 

Consider a linear element AZ carrying an alternating current / = /q exp{icjt). 
The current element may be viewed as two charges Q and —Q oscillating 
back and forth. 

Assume that AZ is much smaller than the wavelength A = li^cjuj. In this 
case we can ignore the phase variation along AZ. 

An understanding of the properties of such an antenna is of great inter- 
est since, in principle, all radiating structures can be considered as a sum of 
small radiating elements. Moreover, many practical antennas working at low 
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z 




frequencies are very short compared with the wavelength. 

If Al <^ A, then there are three spatial regions of interest: 

• The near field (static) zone A/ <C r <C A 

• The intermediate field (induction) zone A/ <C r ~ A 

• The far field (radiation) zone A/ <c; A <C r 

We will see that the fields have different properties in the different zones. 
In the near zone the fields have the character of static fields, with a strong 
dependence on the properties of the source. In the far zone, the fields are 
transverse to the radius vector and fall of as r~-^, typical of radiation fields. 

The retarded current element 

J{t - r/c)dV = T{t - r/c)dl = dlloe'^^''^'^ , 

where k = uj/c. 

Thus, the vector potential is 

^ = ^^T^- ' 

where A/ = J dl. 
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In the near zone, where r <^ A (or kr <C 1), the exponent exp{—ikr) can 
be replaced by unity. In the far field zone fcr ^ 1, the exponential oscillates 
rapidly, and in this region it is sufficient to approximate exp{—ikr) ^ 1 — ikr. 
In the intermediate zone, all powers of kr must be retained. 

In cartesian coordinates the vector potential has only one component, 
say along the current element. 

Referred to spherical coordinate system, A has components 



Aj. = Az cos 6 = 



loAl cos^e^^^*-^^^ 



Ae = -A,sme = —^ , 26 

= 0. 

In order to find the fields E and 5, we have to calculate V • A and V x ^4, 
that in spherical polar coordinates are given by 

^ 1 d jr^Ar) 1 d{Aesme) 



+ 



dr rsin^ dO 

1 dA^ 
rsiu.6 d<j) 



V xA = 



rsm.6 



a (A^ sing) _ dAe_ 



6 

+ - 

r 



r 



1 dAr d{rA^) 



sm9 d(j) 
'dirAo) dA, 



dr 



dr 



09 



Since A^ — Q and there is no 0 dependence of Aj. and Aq^ i.e. OAr^o/dcj) 
the above equations reduce to 



= 0. 



V • A = ^ — + 

<^2 ^)r 



1 d {Ae sin 9) 



ae 
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V X A 



d (rAe) OA 



dr 



d9 



r 



Hence, the magnetic field of the current element is 



5 = Vxl=5^(^, 



where 



Bs — - 
r 



d{rAo) OA, 



dr 



(27) 



and Br — Be = 0. 

Thus, the magnetic field is perpendicular to the radius vector at all dis- 
tances. 

Calculate the magnitude 5^. Substituting Eq. (26) into Eq. (27), we 
obtain 



Brh — 



-JqA/ 

IqAI 



d (sin 0e^M-M^ ^i{u;t-kr) q ^^^q 



dr 



ifcsin 



r de 

^{{ijjt—kr) 



Hence 



B^ — 



IqAI 

AtTSqC^ 



ik 



1 



sin ee^^'^*-^^) 



The magnetic field is composed of two terms: the near zone term ~ 

and the far zone term ~ 1/r. In the limit of a; ^ 0, the near zone term 



reduces to the Biot-Savart formula. The far zone term is only present for 
an oscillating field {uj ^ 0) and therefore it is radiation field arising from 
accelerated (oscillating) charge. 

Calculate electric field of the current element: 

^ c2 dA 
E = -V V ^A)-— , 
10} at 
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where 



and 



OA 

~dt 



V-A = — 



ld{r^Ar) 1 d{^meAe) 



dr 



+ 



r sin 9 dO 



IqM \ 1 5(rcos^e^('^*-^^)) 



AttSqC^ 1 
1 d 



r sin 6^ 06 



sin^ 6^ 



dr 

^{ujt—kr) 



2 sin e cos (9e'("*-'='') 



sin 9 



IqAI 



1 ifc 

-IT + — 



r 



cos ^e^^"'-^"^ . 



Since in spherical coordinates 



„ ^ d -Id 6 d 

V = f h 6* ^ , 

dr r 89 r sin 9 d(j) ' 



we obtain 



To Alcove 



IqAI cos ^ 

47r£:oC^ 

2ik 



ik 



k'^ ik 



r 



J,(ujt—kr) 



V(V- A 



-I ^ 47r£oc^ 



1 ik 

^3 2 



sin^ e*('^t-fe'-) , 



[v(v-l)" 



= 0 . 
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Hence the radial part of the electric field, Er, is 



E, 



V(V • A) 



IqAI cos 6 



dA 



dt 



2 2ik 

,^3 fjp>*2i 



J r 



r 



^i^Lot—kr) 



loAl cos ^ e^'^'*^*-^^) 



IqAI cos 6 

AtTSqC 



2 ik iuj 
+ ^ + 



J>{ijjt—kr) 



Since k — oj/c^ the 1/r terms cancel and then Ej. simpUfies to 



/o A/ cos 9 



ikr 



Similarly, we find the 9 component of E as 



Ee — 



V(V • A) 



dA 



dt 



IqAI sin 9 
/o A/ sine 



1 ik 



J e 

j,{ijjt- 



. Jo A/ sine ^^{u;t-kr) 



1 



1 



ik 



ikr 



r- 



r 



i{ujt—kr) 



The radial part of the electric field contributes only to the near and in- 
termediate zones, whereas the angular part contributes to all of the zones. 

The part is the Coulomb type contribution. It is similar in nature 
to a static field surrounding a small linear-current element and an electric 
dipole. 

Proof: 



The Coulomb or static field is for (j ^ 0, In this limit the contribu- 
tion is 

/oAZsine 1 



Eo — 



Aireoc ikr^ 



J.{ujt—kr) 
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loAlsinO 1 



loAlsinO 1 



1 — ikr H — {—ikr'Y + 



where we have taken only the real (physical) part of the field. 

Since /q = Aq/ At and Al/At = c, we get 

Agsini9 1 

^ ~A 2 ' 

as required. 

Electric and magnetic fields in near and far field zones 

Consider first the near field zone (r <C A). In this limit the magnetic and 
electric fields are 

Bnear = SmO^^^-'-'^^ 

Enear = -i^^ (cOS + sin e^^"*"^^) 

ATreockr"^ ^ ^ 

Since the magnetic field is real and the electric field is imaginary, the Pointing 
vector involving the near-zone field components is a pure imaginary quantity. 
It does not represent any flow of energy. This imaginary quantity represents 
energy that oscillates back and forth between the source and the region of 
space surrounding the source. 

Consider now the far zone or radiation components of the magnetic and 
electric fields. 

AttSqc r 

Brad = B^ci> , B^ = ^ e^(-*-'='-) . 

Note that: 
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1. The electric and magnetic fields oscillate in phase. 

2. The electric and magnetic fields are orthogonal to each other. 

3. The ratio -5^=0, the value for plane waves in free space. 

4. The electric and magnetic fields are transverse to the radius vector at 
all distances. 

5. The Poynting vector N = c^SoErad x Brad is a real quantity and is in 
the direction of the radius vector, indicating that the energy of the field 
propagates away from the current element. 

These properties show that in the far zone the field is in a form of plane 

waves. 

9.2 Power Radiated from the Current Element 

The power flux at any point is given by the Poynting vector 

N = c^SqE X B . 
Then, the total power radiated across a sphere of radius r is 

W = [ N^dS , 
Js 

where 

dS = sin 6 dOd^ . 

Only those partial products m E x B which vary as will have net radi- 
ated power. The other partial products are small as they fall off more rapidly 
than Thus, the only part of the fields entering into the expression for 

the radiated power is the far field zone part (radiation component) consisting 
of the terms varying as 1/r. 
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The volume Idl is energized and the energy flows but there is no net 
energy loss over a cycle in the equilibrium situation. 

The radiation components of E and B are in phase and average over 
time is 



where (^^)o amplitudes of and B^, 

Hence, the time averaged Poynting vector is 



1 



N = - 



^ sin^ 9 



2^ ^^167r%oC^ 



Thus, the total power is 



sm 6 — — 



-^2 



sin^ ^ 



w = 



r / r^iV sin e ded<j> 
Jo Jo 



27r p 
^0 



0 SSqC 



A 



d(f) / si 

^0 



sin^ ec/e 



Integrating, we get 



W = 27r 



^0 



'Af 


^4 


_ TTli 


'Af 


A 


3 




A 



We can write the total power radiated in terms of the power absorbed in an 
equivalent resistance, called the radiation resistance: 



2 3£oc 



A/ 

T 



-^0 = 9-^-^0 • 



where 



2tt 

3£oC 



T 
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is the radiation resistance. 



Since l/{eoc) — y^/xo/^o = 377 or 1207r, we obtain 



R = SOtt^ 



Al 

T 



For example, if Al/X ^ 0.1, then R = O.Stt^ ^ 8 ohms. 

This example shows that for a current element which is 10% of the wave- 



length long, the resistance is very small. 



Thus, if A//A 1, the radiation losses are negligible, that the radiated 
power is very small. 

This also explains why ordinary circuit theory often works so well al- 
though it ignores the loss of energy by radiation from AC currents. 

Appreciable power would be radiated only if the current amplitude /q 

were very large. A large current, on the other hand, would lead to large 
amounts of power dissipation in the conductor, and hence a very low effi- 
ciency. 

We can conclude that current carrying systems that have linear dimen- 
sions small compared with the wavelength radiate negligible power. An ef- 
ficient antenna should have dimensions comparable to or greater than the 
wavelength. 



9.3 Gain of the Dipole Antenna 



A further property of the dipole antenna that is worthy of consideration is 
the directional property of power radiated in different directions. 

The gain or directivity function of a transmitting antenna is the ratio of 
the Poynting flux to the flux due to an isotropic radiator emitting the same 
total power W: 

9t{0, (j)) = — - — , 
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where 



ISO 



w 

Airr 



2 



is the energy fliix uniform in all directions. 
For the infinitesimal dipole: 



w = 



^0 



SCqC 

36:oc 




-^ 2 



sin^ 9 



1 2 



Hence 




Figure 4: The polar radiation pattern of the dipole antenna. 

The directivity function gri^^ (f)) defines a three-dimensional surface called 
the polar radiation pattern of the antenna. The function varies as sin^ 6^ and 
hence the radiation is most intense in the 9 = 7r/2 direction and zero in the 
directions ^ = 0,77. The maximum gain then is 1.5 for directions defined by 
9 — 7r/2, in the equatorial plane of the dipole. The gain function is indepen- 
dent of 0. 



We can conclude, that the directivity function griP^^) is a measure of 
how effective the antenna is in concentrating the radiated power in a given 
direction. 
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Questions: 

(1) Show that in spherical polar coordinates, the magnetic field of a short 
current element lAl = AZ/q exp(i(x;t) has only an azimutlial component of 
the form 

sin e e^('^*-^-) 0 . 

(2) Show that in the far field zone of a radiating short current element, 
the electric and magnetic fields oscillate in phase and are orthogonal to each 
other. 



B = 



IqAI 

Air Sod 



ik 1 



(3) Given the expressions for the EM field of a Hertzian dipole, show 

that the total radiated power from the dipole is 



w = ^ 



Al 

T 



(4) Show that the time averaged Poynting vector of the EM field emitted 
by a short current element is maximal in the equatorial plane of the element. 
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10 Electromagnetic Theory of Polarizable 
Materials 



If an uncharged dielectric (insulator) is placed in an electric field, the field 
within the dielectric will be found to be modified by charges induced by the 
external field. The charge is induced by rearrangement of bound charges 
within the molecules of the dielectric. In dielectrics these charges are a set 
of molecular dipoles. The dipole may exist permanently or may be induced 
by the external field. Dielectrics with permanent dipole moments are usually 
electrically neutral due to random orientation, in the absence of electric fields, 
of the dipole moments. An example is H2O. We will study electromagnetic 
theory of polarizable materials in terms of dipoles, which we will treat as 
building blocks of dielectric materials. 

10.1 Potential and Electric Field of a Single Dipole 

Mathematically, it is convenient to deal with the dipole not as just a pair of 
individual plus and minus charges, but as a separate object on its own rights. 

Suppose that opposite charges ±g are separated by a distance d. We 
will find the potential $ at a distance r and angle 6 under the assumption 
that r ^ called the dipole potential. 



+q d -q 




We define dipole moment as the product of the charge times the separa- 
tion 

p = qd . 
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Since dcosO at r ^ we get for the potential at point A 



Airsori 47r6:or2 
9 r2 — ri q d cos ^ 1 p - r 



(28) 



We can find electric field of a dipole using E = — V$. Since the potential 
of the dipole depends on r and 6^ it is convenient to work in the spherical 
coordinates in which the electric field is given by 



E^Erf + EeO + E^(j) , 



where 



Ej. — 



(9$ _ 1 2pcos9 
dr ^ttSq r 



3 



1 (9$ 1 p sin 6^ 

hjQ ^ ^ 

r dO 47r£:o 
r sin U o(p 



Hence 



E = , ^ , f2cosgf + sing e] . 
47r£:o^ ^ ^ 

Figure 5 shows a sketch of the electric field lines of an electric dipole moment. 




Figure 5: The electric field lines of a dipole moment. 
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10.2 Polarization Vector 



If there are N dipoles per unit volume, the dipole moment per unit volume is: 

i=l 

and is called the polarization. 




The electric potential set up by an arbitrary volume distribution of electric 
dipoles can be calculated by using the potential produced by a single dipole, 
Eq. (28), and the above definition of P. The electric potential at an arbitrary 
distance r from a volume element dV containing such dipoles is 

1 P-r 
rf$ = - -dV, 

where r is the unit vector from dV toward X, and we have assumed that r 

is much larger than the extent of the volume element dV . 

Let r be the unit vector from X toward dV , (We want to integrate over V 

with the position of X fixed). In this case, we change f — f, and obtain 

rf$ = -^P- l-^r 1 dV 
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The result can be transformed into one that has an interesting physical 

interpretation. 



Noting that 



f „1 



we have 



p.i_4Up-vi. 



Applying a vector identity V • {^A) — $V - A + A - V^, we thus obtain 



IT J \ T" \ 



Hence 



Transforming the first term using Gauss' divergence theorem, the potential 
becomes 



S V 



On comparing of Eq. (29) with the general form of the potential, Eq. (20), we 
deduce that the first term on the rhs of the above equation is the potential of 
a surface charge density P • h. The second term is the potential of a volume 
charge density — V • P. 

10.3 Maxwell's Equation for V • ^ in a Dielectric 



In general, the electric field in a dielectric can be found from the Maxwell's 

equation I: 

So 
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+a -a +a -o^ 

Figure 6: Surface (left picture) and volume charges (right picture). Surface 
charges exist because there are no neighboring charges at the end surfaces of 
the material to cancel them out. Volume charges exist because the number 
of dipoles per unit volume changes, that there is an incomplete cancellation 
of charge density from adjacent dipoles. 

In the dielectric it is convenient to (mentally) separate the polarization 

charges from whatever other charges might be there also. The other charges 
are usually referred to as the free charges or conducting charges to distinguish 

them from the bound charges in the dielectric. 
We can write 

60 V 'E^ pf+pp , 

where p/ is the "free" charge density and pp is the polarization charge density 

throughout the volume. If we express pp in terms of P [pp — —V * P), we 
have 

Hence 

V-(so^ + P)=P/. (30) 

Now it is common practice to drop the subscript /, but one must remember 
that p now stands for the charge density not counting the polarization charges. 
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It is common practice to define 

D = eoE + P , (31) 

where D is called the dielectric displacement 

The reason for the name "dielectric displacement" can be easily under- 
stood if we refer to the Maxwell's theory. 

Take time derivative of both sides of Eq. (31): 

dD _ dE dP 

~dt ~^''~dt^~dt ' 

We know from the Maxwell's theory that the first term on the rhs of the 

above equation represents displacement current density, and the second term 

— * 

is the polarization current density. Therefore, dD /dt can be called a gener- 
alization of the displacement current density, and then D can be regarded as 
the dielectric displacement. 

Introducing the dielectric displacement, Eq. (30) becomes 

V-^ = p . 

We can read this equation that the source for the field D is the free charge 
density p. 

We turn now to a consideration of the macroscopic effects of the polariz- 

ability of dielectric materials. We will consider only ideal dielectrics. 

Ideal dielectrics can be divided into following categories: 

1. Homogeneous - properties independent of the position. 

2. Isotropic - properties independent of direction. 

3. Linear - polarization proportional to E. 

4. Stationary - properties independent of time. 



143 



Case of simple isotropic and linear dielectrics 

Ordinary dielectrics (glass, teflon, plastics etc.) are linear in polarization 

for fields not strong enough to cause dielectric breakdown i.e. P (x E. 
For these materials, we can write 

where x is called the dielectric susceptibility. Then, we can write 

D = €oE + P = £o(l + x)E = SoSrE = sE , 

where Sr is the relative permittivity or dielectric constant^ and e is the per- 
mittivity. Hence 

V • ^ = V • (sE) = p 

or 

e 

if e is independent of position i.e. if £: is a permittivity of a homogeneous 
dielectric. 

The same result would be obtained by replacing Sq in Coulomb's law by e. 

r — r 

Atts r^ 

The ratio s/sq then represents the relative shielding of qi from q2 by the 
polarization charges induced in the medium. 

The theory of the molecular structure of a material will yield an estimate 
of X ^iid e. The dipole moment p of a molecule will be proportional to the 

local electric field E so we can define a molecular susceptibility Xm such that 
V = XmSoE. 
Then 
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where N is the number of molecules per unit volume. 
Thus 

e = SoSr = 6:0(1 + Nxm) 
and s is a quantity directly measurable from the measurement of capacitance 

~ ~d ~ ~d~ ~ ' 
where Co is the capacitance without the dielectric. 

We can summarize, that filling a capacitor with dielectric multiplies its ca- 
pacitance by Sr^ 



Exercise in class: Capacitor filled with a homogeneous dielectric 

A plane parallel capacitor has charges +cr and —a per unit area on its plates. 
The capacitor is filled with a homogeneous and linear dielectric of dielectric 
constant 6:^ = 1 + 

Show that: 

1. The electric field within the dielectric is: 



€rSo 



2. The polarization charge per unit area on the surface of the dielectric 

adjacent to the surface of the negatively charged plate is: 



1 + X 



3. The capacitance of the capacitor is C = SrCo where Co is the capaci- 
tance of the same capacitor without the dielectric (i.e. a vacuum or air 

between the plates). 



145 



10.4 Dense Dielectrics: The Clausius-Mossotti Rela- 
tion 



The Lorentz theory of polarizabiUty of dense dielectric materials distinguishes 
between the mean electric field E and the local electric field Ei^^ seen by a 
typical dipole. The tj^ical dipole is considered to be at the centre of a small 
sphere that has been excavated from the dielectric. Ei^c is thus the mean 
field E minus Epiug where Epi^g is the field of the spherical volume excavated. 
We will show that: 

Eioc = E + 



3eo 



If there are N molecules per unit volume each of polarizabiUty a the 
poarization vector: 

P = NaEioc = Na\E+ ^ 



and by definition: 

Let E be the mean field throughout the dielectric. 
Let Epiug be the field due to the spherical plug alone. 
Let Eioc be the field in the spherical hole. 

^ ^ ^ ^ P 

E = Eioc + Epiug and from above Epiug = - — 

Thus 

Ejoc — E + 



3^0 



The argument now is that each molecule is at the centre of a small hole 
and the field acting on the molecule is thus Eioc If ^ is the molecular 
polarizabiUty, its induced dipole is thus: 

p = aEioc 
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If there are N molecules per unit volume then: 



P = Np = NaEioc = Na\E-^ 



P_ 



By definition P = x^o^ = i^r ~ ^)£oE. Substituting for P: 



(Sr — l)eoE = Na 



{Sr - l)eo = Na(l + 



3eo 

£^ — 1 



Na 



{Sr + 2) 



— 1 



Na 



which is known as the Clausius-Mossotti relation for a dense dielectric. 



We can solve the Clausius-Mossotti equation for Sr^ and obtain 

Na 



Sr — I = {Sr + 2) 



1- 



Na 
3^0 - 



= 1 + 



2Na 



Note that 



1 + 



2Na 
3eo 



1- 



oo as 



Na 
3^ 



1 , 



which is called the Clausius-Mossotti catastrope. 



Removal of the plug leaves polarization charges, whose field tends to line 

up the dipole parallel to the external field. The system is self-polarizing 
(Clausius-Mossotti catastrope) if ^ ^ 1. 



3^0 



10.5 Time Dependent Fields and the Complex Dielec- 
tric Susceptibility 

The induced polarization charges do not produce any currents inside the 
dielectric. There is no DC current in response to a DC electric field, but if P 
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is changing with time (because E is changing with time) there will be an AC 
current density: 

r das . dP 

dt dt 

Thus, dP /dt plays the role of polarization current density. 

Let us consider what will happen if a dielectric is introduced into a time- 
varying electric field. 

With an alternating electric field (due say to an AC voltage across a ca- 
pacitor filled with a dielectric) the polarization P may lag in phase behind 

— * 

the driving field E. This means there is internal friction and heat dissipa- 
tion as is discussed below. The "capacitor" will exhibit resistive as well as 
capacitive properties. 




Consider the work done in charging a capacitor. We start from the circuit 
theory quantities and will express them in terms of the field quantities. 
The work done in charging the capacitor is 

dt dt ' 

where V is the voltage. 
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Since 



Q = CV = CEd 



eA 
d 



Ed , = AeE = AD 



we obtain 



Q = A{€oE + P) , 



and then 



^ = Ed^MeoE + P) = EdA^(eoE + P) 
at at at 

Since dA = V is the volume of the capacitor, we can write 



dW 
~dt 



^d , dP 



dV , 



or 



dW f d /I r 
= I - i T^eoE^) dV + J 



dt 



dt \2 



V 



V 



where we took into account a possibiUty that the electric field and polariza- 
tion can vary across the capacitor's plates. 

The first term in the above equation is the rate of doing work building 
up E field. The second term is the rate of doing work on the dipoles by E, 

Thus, the supplied energy to the capacitor is used to build up the electric 
field inside the capacitor and to polarize the dielectric. Consider separately 
both terms. 

First term: 



li E = Eq cos{(jt)^ the first term takes the form 



J ~Eo cos (cut) SqujEq sin (cut) dV 



V 



The element of work done per unit volume unit time is 



dW 



= —SqujEq cos (cut) sin (cut) 
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Averaging over a cycle, we get 

27r/u; 

dW T-,2 / / N . / X 7 

-— - = —SqujEq / cos(ct;t) sm[ujt) dt = 0 , 

Work is done building up the field in one part of the cycle but the stored 
energy is given back in another part. 

Second term: 

If P = x^oE = x^o Eq co8{ujt) the same zero net energy conversion aver- 
aged over a cycle will happen with this term. If there is internal friction 
there will be a phase difference between P and E. 
Write 



P = x^oEo cos{u;t + (j?) 



Then 



P = X^oEq cos^ cos{ijt) — x^qEq siiKf) sin(ct;t) 



and 



dP 



(JX^qEq co8(j) sin(ct;t) — (jx^oEq sin^ cos(cc;t) 



dt 



Hence, the work done per unit volume per cycle will be 




27r/u; 



0 



27r/u; 




0 



Since the first integral is zero, we get 



dW 



dV 



27r/u; 

ojXSqEq smcj) J cos^{ujt)dt 

0 
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The integral on the rhs of the above equation is positive and dW/dV must 
be positive corresponding to energy dissipation (or the dielectric would keep 
getting energy from its interior and building up the field with it). 



Thus, sin0 must be negative, so — tt < 0 < 0. 



10.6 The Complex Susceptibility and Permittivity 

In the usual way use complex exponentials to represent amplitude and phase, 
We can write 

E = Eoe^""^ then P = Poe^^"^*"^^ {(f) is positive) . 
We can write the complex polarization in different forms 

P = Poe"^^ e^"^* = (Po cos (J) - iPo sin (;/>)e^'^* , 

P = (^cos0-i^sin<^)Eoe^-S 

P = eo{x'-^x")Eoe'^'-eoXcEoe'^'^eoXcE, 
where Xc = ~ ^x'' is a complex susceptibility. 

With the complex polarization, the dielectric displacement takes the form 

D = eoE + P = eoE + SoXcE = £o(l + Xc)E , 



or 



D = eo{l + x'- ixlE = SoSrE = ScE , 

where Sc is a complex permittivity, and Sj. is a complex relative permittivity 
or dielectric constant 

£c = £0(1 + XO - ^^OX" • 
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Example - equivalent circuit of the lossy capacitor 

The imaginary part of the dielectric susceptibility corresponds to net en- 
ergy dissipation i.e. in the circuit theory language the dielectric susceptibility 
adds a resistive component to the capacitor. The material filling the capac- 
itor could also have some ordinary ohmic conductivity (due to the presence 
of 'free' charges as well as 'bound' charges in the material). Let us calculate 
the magnitude of the resistance. 



+Q = SA 



E 



-Q = -SA 



Let / be the charging current in the external circuit, and 

r - - dP 

I, = jj,-dA = A- 

be the polarization current in the dielectric. Let 

f ^ ^ V 
Ic = j Jc' dA^ AaE = Aa— 

be the conduction current in the dielectric due to its finite conductivity a. 
If S is the charge per unit area on the plates then: 

S — (charge supplied by /) 

— (charge removed by Ip) 

— (charge removed by 7^) . 



Then 



Q = SA^ J Idt- J Ipdt- J Ic dt 
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and 



S V 
So d 



Hence 



S — Cq— 

a 

and also remember that P = x^oE — X^o^? we obtain 

Q = = A = J Idt- J Ipdt- J Icdt 
or taking a derivative in time, we get 



^0 



AdV 

d dt 



I-Ip-Ic 



= I-A 



XSq dV .aV 



d dt 



-A 



d 



Putting 



dV 
'dt 



= iuV and X = x' ~ ^ x" 



and solving for /, we get 



/ = 



d 



a 



(1 + x' - ix")^^ + - 



V 



Separating real and imaginary parts and putting Co = e^A/d (the capac- 
itance there would be if the dielectric were lossless): 



V 



Since the capacitor transmits some charges through the internal dielectric, 
in the circuit theory this system is equivalent to a parallel circuit 
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Comparing with the above result for current flow in the lossy capacitor we 
see that the effective capacitance is (7o(l + xO the effective resistance is 



where we remember that Cq is the capacitance in the absence of losses. 



10.7 Added Note - The Loss Tangent 

The properties of a dielectric material are usually specified by giving its di- 
electric constant and its loss tangent tan S, 



If we write: 



(1 +x')-^ {x" + 



a 



V 



The quantity in [ ] brackets can be defined as a complex relative permittiv- 
ity Sj.: 

er = l+x'-^{x" + —) . 

This is a generalization on the previous definition of complex relative 
permittivity to include the effects of ohmic conductivity. Then we can define 
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a generalized permittivity and dielectric constant e = SqEi 
Now if we write: 



— cos 5{1 — i tan 5) — K{1 — i tan 5) 



Then the standard form is 



e — £0^(1 — tan 5) , 

where tan 5 is the 'loss tangent' and is read out on some AC bridges as an 
alternative to reading out the resistive property of a lossy capacitor. 
Since K' cos 5 = (1 + xO and sin 5 = (x'' + ^) it follows that: 



tan 5 — - 



(J 



+ 

In this equation, tan 5 includes the effects of finite conductivity and the effects 
of polarization damping force. 



Questions: 

(1) Show that the electrostatic potential due to a distribution of electric 
dipoles of moment per unit volume P throughout a volume V enclosed by a 
surface S is that of a volume charge density — V • P together with a surface 
charge density P • n. 

(2) Show that the polarization of a dielectric driven by a time varying 
electric field lags in phase the driving field. 

(3) Show that the phase difference between the polarization and a time- 
varying electric field results in a complex permittivity of the dielectric. 
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11 Magnetic Fields in Magnetizable 
Materials 



We have previously discussed how the polarization of dielectrics by an exter- 
nally applied electric field is equivalent to volume and surface distribution 
of charge. Analogously, a magnetic field can act on molecular scale current 
loops to produce macroscopic effects. It was Ampere who first suggested that 
the magnetism of matter was due to the cooperative effects of currents circu- 
lating in atoms (and not, as previously thought, due to a separate magnetic 
charge called poles). 




11.1 Magnetic Polarization Currents 

We define the magnetic moment of a current loop as 

fl = I Ah , 

which is equal to the product of the area of the plane loop and the magnitude 
of the circulating current. The vector direction h of the moment is perpen- 
dicular to the plane of the loop and along the direction set by the right-hand 
rule. 

In a material body we define a macroscopic dipole moment per unit vol- 
ume 

i 

which is called magnetization. 
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Theorem 



If a volume V enclosed by surface S has a magnetic dipole moment per 
unit volume M (which may be a function of position), the macroscopic mag- 
netic fields so produced are equivalent to those of: 

• A volume current density Jy = V x M. 

• A surface current density J5 = M x n. 
Proof: 



The theorem is proved by showing that the vector potential due to the 
dipole distribution in a volume V closed by a surface S can be written in the 
form 

V s 

From the solution of the Maxwell's equations (static) vector potential of 
a current loop is given by 

Consider a current loop of radius a, shown in Fig. 7. 
In polar spherical coordinates 

^ /N /\ 

dl — ad(f) (f) — —a sin (j) d(j) i + a cos (j) d(j) j , 

However 



r — 



[x — a cos (f)) + (y — a sin (j)) + z 



1/2 



which for a <C i? can be written as 



r — (x^ -V + — 2ax cos (f) + — 2ay sin (j)) ^ 



(r^ — 2ax cos <p — 2ay sin 4>) 
ax cos 0 ay sin 0 



1/2 



R 1- 



R^ R^ 
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Figure 7: Current loop of radius a. 



Hence 



T — R 1 H h 



where we have used the Taylor expansion of 1/(1 — a?) = x -\- 
Thus 

/ r^TT / ax cos (j) ay sin (j) 
AttSoC^R Jo \^ R^ ^ 

X i^—a sin (f) dcj) i + a cos (j) dcj) . 

Since 

p27r r27v r27v 

/ sin (j) dcj) — cos (j) d(j) — l sin ^ cos (j) dcj) = 0 ^ 
Jo Jo Jo 

the formula for A simplifies to 

^ = / f-a^ysin^^ i + a^o^cos^^ j) (i0 , 

Attsoc^R Jo ^ ^ 
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which can be written as 



T 

A 



Next, since 



r27r r27r 

■yi I sin (j) d(j) + xj / cos (j) dcj) 
Jo Jo 



p27r /•27r 

/ sin^ (j) dcj) = cos^ (j) d(j) — tt ^ 
Jo Jo 



we obtain 



A = 



Using the relation 



k X R— k X 



R R 



we get 



where we have used the result 



R_ 



If we consider a set of dipole moments and change, for a convenience, the 
direction of R into —R^ we can write that the vector potential dA produced 
by a set of dipole moments contained in a volume dV is 



dl = M xV- dV 

Attsqc^ R 



where M is the magnetic dipole moment per unit volume. 



159 




Then, the vector potential produced by the whole dipole moments 

tained in the volume V is 



1 

R 



Using a vector identity 



V X {(j)A) = V(j)X A-\-(j)V X A 



we have 



_ 1 „ 1 V X M _ M 

M X V- = -V— X M = — V X — 

R R R R 



and then the vector potential is 



A = 



1 



AtTSqC' 



V X M 
R 



dV-—^ /Vx 



M 

'r 



dV 



In order to proceed further, we introduce a Theorem 



X 



, „ M f Mxn 



This is an application of the more general theorem 



- VxFdV= FxhdS 
V Js 
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Proof: 

Let C be a constant vector. Then 

V . (F X C) = (V X F) . (7 - (V X C) • F 

= C-(VxF) 

We will prove the general theorem by using the divergence theorem 

/ V^(FxC) dV = [(FxC)^hdS 
Jv Js 

For an arbitrary constant vector C, and using (32), we get 

C' f V xF dV = [ C'V xF dV 
Jv Jv 



V'(FxC) dV^ / FxC'fidS 

'V Js 



Hence 



/ V xF dV = - C xF^hdS 
V Js 



However 

C X F 'h^C ' F xh 

and then we obtain 

[ VxF dV = -C' f FxhdS 
Jv Js 

Since this is true for arbitrary (7, we finally have 

V xF dV ^- I FxhdS 
v Js 

as required. 



Thus 



- 1 r V X M 1 f Mxn 

^ = -A 2 / — ^ + / — -- dS , 
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or 



A = ^ f — dV + ^ ( — dS 
Attsoc^ Jv R Aitsqc^ Js R 

The effective (Ampere) currents associated with a dipole moment per unit 
volume M are 

(i) A current density Jy = V x M throughout the volume. 

(ii) A surface currents Js — M x h. 




— * 

11.2 The Magnetic Intensity Vector H 

When we were dealing with dielectric materials in the presence of electric 
fields, it was convenient to introduce the displacement vector D in order to 
eliminate the necessity of taking the electric dipole polarization P of the 
material into account explicitly. 

A similar procedure is used for the magnetic materials, and we will illus- 
trate it here for both static and time- varying fields. 

For a static current distribution, the Maxwell's equation IV reduces to 

V X 5 = ^qJ . 
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In a medium where there are magnetic polarization currents as well as con- 
duction currents we can write 



J — Jc ^ Jm 5 

/ = X + V X M . 

Hence 

V X 5 = /io(X + V X M) , 
which can be written as 

V X (5-/ioM) = /io/c , 

or 

V X I — - M I = X 

This shows that the vector B/ fiQ — M has as its source only the conduction 
current J^. Therefore, to eliminate the necessity of dealing directly with the 
magnetization M, we can define a new vector 

H = — -M ^ 

which is called the magnetic (field) intensity vector. 

Then the Maxwell's equation IV for static fields in magnetic materials 
takes the form 

V X ^ = X . 

In dealing with magnetic materials we often know Jc but not M (well not 
directly anyway.) Think e.g. of an inductor filled with some magnetizable 
material like iron. Then H becomes useful. It is a way of avoiding a detailed 
calculation of the polarization currents. The magnetic intensity H is the 
magnetic analogue of the dielectric displacement D in the electric case. We 
may drop the subscript c and write 

V X # = /, 

but we should remember that J is now not the total electric current density 
everywhere. 
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11.3 Linear Isotropic Magnetic Materials 



For most materials (excluding ferromagnetics) the magnetization M is pro- 

— * 

portional to the applied external field. Hence, at any point the vectors jB, 
M, and H will be in the same direction and we may write 

B = ^lQ^lrH = /io(i + Xm)H = ^lH , 

where ji^ is the relative permeability^ Xm is the magnetic susceptibility^ and 
fi = /xo(l + Xm) is called the magnetic permeability. 
Since 

B = tio{H + M) , 
we have defined Xm such that 



f} fJ^oXmB XmB 



Hence, we get 

B 



iUo(l + Xm) ' 



Thus, if we know the material we use, we can find H, 

Example A solenoid filled with magnetizable material 

In the expression B — ^qNI we should include the Ampere surface cur- 
rents as well as the conduction currents in the wire. Since, B — fi^I , where 
/ is the total current per unit length 

B = fioiNI + M)= fioiNI + XmH) = /^(iV/ + ^) , 



which can be written as 



A'' 
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or 

Note that if Xm is positive then B is greater than it would have been in the 

absence of the magnetizable material. Evidently in this case the macroscopic 
Ampere current is in the same sense as the conduction current in the solenoid. 
Prom the definitions 



fl /io(l+Xm) 

1 _ ^ 1 ^ ^ 

1 ~l~ Xm 1 ~l~ Xm /^r 



B = firl^oNI — fiNI = fiH . 

The effect of filling the solenoid with a material of relative permeability fir is 
to multiply 5 by a factor fi^ (assuming the current in the wire remains the 
same) . 

Note: H is independent of the presence or absence of the magnetic ma- 
terial. 



Consider a solenoid filled with a magnetic material. By the definition 

= — -M . 

When we remove the magnetic material, M = 0, and then 
With the material present 

H^ — -M = ^- XmH = flrNI - XmH . 

Hence 

H = {l+Xm)NI-XmH 
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from which, we find 
and then 

H ^ NI 

as before. 

This gives rise to the notion of H as an inducing field and 5 as a resul- 
tant field. This concept is much used in the study of magnetic properties of 
materials. 

The effect of magnetic materials on inductors 

Because the relative permeability multiplies the magnetic field by /x^ 
for the same current, the self inductance of a solenoid is multiplied by /i^. 




Let is a number of turns/meter. Hence, Nl is the number of turns, 
and then the flux is 

= Lol = nBoA = (iV£)(/ioiV/)A . 

Thus 

with relative permeability ^r-, so that B ^ fJ^rBo- 



166 



11.4 The Magnetization of Iron 



Iron is not a linear isotropic material. It is subject to saturation of the 
internally produced magnetic fields because all the internal current loops are 
lined up. Iron does not have a unique value of magnetic susceptibility because 
of strong nonlinearities. In ferromagnetic materials the relation between B 
and H is usually presented graphically in terms of hysteresis. 




Work done in magnetization 

Think of the case of the solenoid of length £ and cross-section area A filled 
with magnetic material. The applied voltage V is: 

dt dt 

where Nl is the number of turns. 
The rate of doing work is: 

P = VI = AiNI^ = V 

dt dt 

where V = Al is the volume of the material. 

In the time dt that it takes to change B by dB^ the work done is dW = 
Pdt = V H dB. Thus, the work done per unit volume is dW — H dB, 

Hence, the work done per unit volume in one cycle of the hysteresis loop 
of a non-linear material such as iron is given by the area of the loop. 
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IRON 



FERRITE 



It is sometimes useful to write 



HdB = 



Hoj Hd{H + M) 

= J H dH + pLo J H dM 

= ^fioH^ +f^o J HdM , 

where the first term on the rhs is the work to establish magnetic field, and 

the second terms is the work by the field H to establish magnetization dM. 



11.5 Time Dependent Magnetic Fields and Energy Loss 

In a time-dependent magnetic field the magnetization M may not stay in 

phase with the driving field H. This corresponds to internal friction and 
heat dissipation. 



Let 



H = Hq COS^Ljt) . 



Then 



M = Mo cos{ujt + (j)) , 
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where Mq and (j) represent the ampUtude and phase of the magnetization 

response to the magnetizing field. 

Thus 

M — Mq cos (j) cos(u;t) — Mq sin 0 sin(u;t) . 

Now the work done in magnetization per cycle of the AC current producing 
the magnetizing field is 

27r/u; 

dM 



/r rfiyi 
HdM^fjLo J H-^d^ • 



Since 

dM 



— —ojMq cos(f) sin(cc;t) — cjMq sincj) cos{u;t) , 

(jLL 

we have W = Wi + W^^ Consider the term W\\ 

W\ — —fiQEoMQij COS (f) J cos(ijt) sm(ijt)dt . 

t=o 

This term is zero - it represents reversible energy conversion to and from H, 

Consider now the term W2: 

271 /uj 



W2 = —fJ^QHQMQCJsincj) J QOS^{ujt)dt , 



0 



This term represents work done against internal friction during magnetizing 
and demagnetizing the material. 

Since / cos^(cc;t) dt is positive, sin^ must be negative so that work is done 
on the material i.e. (j) is negative. 

The complex magnetic susceptibility etc 

Let 
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which for a physical field can be written as 

H^Hq cos{ujt) = Re Hoe"""* 



Then 



M = Mo e'^""'-^^ = Mo e-'^ e' 



iuut 



or 



iujt 



M = (Mo cos ^ - iMo sin 0)e 

— — cos (p — I—— sin (p I iioc 

This result can be written in terms of real and imaginary susceptibility 

M = (x' - ^x")^oe^"* = (x' - ^x")H . 
Using this complex number notation, we find 

B = fio{H + M)=f,oH + f,o{x'-ix")H 

= [/io(i + xO- W]^ 

and then 

5 = (/i — i/x )if = /iif , 

where fi' = /xo(l + xO? /^'' = mo"x^\ and /x is the complex permeability. 



11.6 The Ferromagnet 

We have so far considered magnetic materials, diamagnetics and paramag- 
netics in which the magnetization is a function of the external field, i.e. 
M B. There is a class of materials, called ferromagnetics in which macro- 
scopic magnetization exists even in the absence of the external field. 
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Consider a homogeneous ferromagnetic material. In this case B inside the 
material is due solely to the M 

or in terms of the magnitudes 




Using Ampere's law § B ^ M — /xq/, we obtain 

Bt = /ioM£ i.e. B = /ioM 

Thus 

H ^ M = 0 

in the region where M 7^ 0, i.e. inside the magnet. 
Outside the magnet 

B B 

H= 0= — , 

Mo Mo 
i.e. H is just a scaled replica of B, 
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Note: 

V • 5 = 0 always and the lines of B form closed loops. Then 




Thus H has a source (field lines start and stop) where M varies i.e. at the 
ends of the magnet. 

For a ferromagnet 

V • # = - V • M = , 

so we can think of volume density of 'magnetic charge' giving rise to 

the H field. It must be stressed that this equivalence is purely mathematical, 
and does not prove a physical existence of magnetic charges. 
Moreover, for the ferromagnet 

V X ^ = X = 0 . 

Note the similar mathematical properties of H here to those of the E field 
in electrostatics 

V-^=— and Vx^ = 0. 

Historically, magnetostatics of permanently magnetized materials (ferro- 
magnets) developed via use of the H field and 'magnetic charges' or 'poles'. 
One obtains a law analogous to Coulomb's Law for the force between mag- 
netic poles. 

Exercise in class: Plane magnetized material 

An infinite plane surface divides the universe into a vacuum on one side 

and a magnetic material on the other. Within the magnetic material there 
exists a uniform magnetic moment per unit volume M which is parallel to 
the surface. 
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(a) Show that while the direction of the magnetic induction vector B is 
different on the two sides of the surface, its magnitude is given everywhere by: 



B 



M 



(b) Find the magnitude and direction of the magnetic field B everj^where 
due to an infinite plane parallel slab of material of thickness d which is perma- 
nently uniformly magnetized with dipole moment per unit volume M lying 
parallel to the bounding surfaces. 

(c) Find the magnetic intensity H everywhere. 

11.7 Maxwell's Equations in Dielectric and Magnetic 
Materials 

The Maxwell equation containing a current density term is 



We think this is always true provided J is the total electric current den- 
sity. Applying this in a region where there may be electric and magnetic 
polarization effects we can write 



dE 



J=J^ + Je + Jm = J+ — + VxM. 



Thus 



V X 5 = floJc + fJ^o 



dp 



+ /xoV X M + fioSo 



dE 



dt 



dt 



V X 
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and finally 

^ ^ dD 
at 

where we must remember that J represents the conduction current only. 



In summary: The Maxwell's equations in materials are of the following 
form 

V-D = p, 

V B = V-H = 0 , 

W X E = -n 



Vxif = J + 



dt ' 
dP 
'dt ' 



but in general 



These equations are supplemented by appropriate constitutive relations which 
connect the electric field E and the magnetic induction B with the displace- 
ment field D and the magnetic field H 

D = eE and B = fiH , 

These relations carry information about the material. 



Questions: 

(1) Show that inside a ferromagnet H — Q. 

(2) Derive the Maxwell's equations for the EM fields in electric and mag- 
netic materials. 
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12 Poynting's Theorem Revisited 



We have seen in Section 6.1 how energy of the electromagnetic field may 
be transported through vacuum (empty space) by means of electromagnetic 
waves. We have shown that the direction of propagation of energy is deter- 
mined by the Poynting vector. In this lecture, we will reconsider the Poynting 
theorem taking into account propagation of the electromagnetic field in mag- 
netizable materials. 



12.1 Poynting Vector in Terms of E and H 

In a vacuum 

^ B B . ^ 
i7 = — = — = eoc^B . 

Thus the Poynting vector is 

TV = Soc^E xB = ExH , 

The cross product E x H also turns out to be the correct expression for 
the Poynting vector when magnetizable materials are involved and is the ex- 
pression most commonly quoted for it. 

Consider, as before in Section 6.1, a flow of the energy through a surface S. 

Using the Gauss's theorem and a vector identity 

V • (1 X 5) = 5 ♦ (V X 1) - 1 ♦ (V X 5) , 

we obtain 

J^{E xH)'dS = J^V'{ExH)dV 

= J H ' {V X E) dV - J E ' {V X H) dV 
Now substitute from the Maxwell equations 
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and obtain 

j{ExH)^dS^- j ii — dV - j E~ dV - j E'JdV , 

The Ihs is the rate of flow of field energy out of volume V , 

1. First term on the rhs is the rate of work in establishing the magnetic 
field in V . 

2. Second term is the rate of doing work in establishing the electric field 
in V, 

3. Third term is the rate of doing work on the currents in V , 

Note that the rate of doing work in establishing the magnetic field is 

f ^ dB 

— # 

Thus, to change the field by an amount dB we have to do work per unit 
volume 

dWu = H'dB 

and then the total work is 

Wm = J H'dB = J d{fioH + /xqM) 

The first term on the rhs is work in energizing field, and the second term is 
work in aligning magnetic dipoles. 

In a previous discussion we saw that / H • dM can involve loss of energy 
from the field to the material (c.f. magnetization of iron). A similar result 
holds for the electric field. 
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We go further with Poynting's Theorem 
Put B — fiH^ then we get 

Hence 

-J E'JdV . 

Thus, the rate of flow of field energy out of volume V is equal to the rate of 
changing energy of the EM field plus the arte of doing work on the currents 
in V, 

12.2 Poynting Vector for Complex Sinusoidal Fields 

It is well known that the electromagnetic field (e.g. light) is a real physical 
quantity (observable). However, in the electromagnetic theory it is advanta- 
geous to represent the real electromagnetic field by complex sinusoidal quan- 
tities because of its mathematical simplicity. In addition, what we usually 

measure is the average intensity of the field, {E* -E)^ which is a real quantity. 
We usually write 

E = Eoe^"* and H = FqC^^* , 

where Eq and Hq are complex quantities including both amplitude and phase 

information. We understand that the electric and magnetic fields are given 
by the REAL PARTS of E and H, 



dtJ 2 u 
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The power of the complex exponential scheme lies in the fact that for oper- 
ations such as summation, subtraction, integration etc., we take real parts 
AFTER the operation. For example: 

ReEi + ReE2 = Re (Ei + E2) 
^ dEi ^ dEo ^ d / ->■ \ 



Some care has to be taken in evaluating the Poynting vector. 

The Poynting vector N — E x H ^ but if we write complex exponential ex- 
pressions for E X we must remember that 



N = ReE^ X ReHc ^ Re (e^ x , 

where we use a subscript c to indicate that we are writing a complex expo- 
nential. 

Proof: 

We can write 

Ec = E^e'"^^ = (Er + iE^ (cos ujt + i sin uot) 

— {Ej. cos ut — Ei sin cot^ + i {Ej. sin ujt + Ei cos ut^ . 



Similarly 



He = Hq^"^^ = [Hr + iHi) (cos ujt + i sin ujt) 

— [Hr COS cot — Hi smut) + i [Hr sincut + Hi cos ujt) , 



where Er^ Ei^ Hr and Hi are real vectors. 
However 



ReEc = Er cos (x;t — Ei sin c^t 
Re^c = Hr cos ct;t — Hi sin ct;t 
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Clearly, if we calculate Re(^^c x He) we get extra terms in addition to those 
in the expression for KeEc x ReiJ^, and then 

ReEc X ReHc ^ Re{Ec x He) , 

as required. 

There is however a useful expression for the MEAN POYNTING VEC- 
TOR in terms of the complex exponential Ec and He 

^ =^Re{E,x H:) = ^Bb{e: X H,) , 
where "bar" over means average over whole cycles of the sinusoidal field. 
Proof: 
Since 



N = ReEe X ReHe 
= (^Er COS (jt — Ei sin ojt^ x i^Hr cos ojt — Hi sin uj-t^ 

= Ej. X Hj. cos cot + Ei X Hi sin cot 
— {Er X Hi + Ei X Hr) cosLjtsinojt 



and 



we obtain 



I f'T 1 

— / cos^ ujt dt — — 
T Jo T 

[ cos cc;t sin cc;t dt = 0 
Jo 



7V= ^{Er xHr + E,x Hi) . 



1 

/ sin^ ujt dt ^ - 
Jo 2 



On the other hand, take 

H: = (Hr - iHi) e-'^\ 

Ee = (Er+tEi)e'-\ 
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and then 

{Ec X HI) = {Er xHr + EiX Hi) + i{Ei x Hr - Er X Hi) . 

Hence 

^Re (e.xHI) = ^{ErXHr + EiXHi) ,= 7V 

as required. 

In summary: The average Poynting vector N of complex exponential fields 
satisfies the relation 

N = ReEc X ReHc = ^Re (^^ x H^) . 
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13 Plane Wave Propagation in Dielectric and 
Magnetic Media 



Now, we shall examine in some detail how existing radiation is modified by 

the material it passes through. We will find that the conductivity is the most 
significant parameter. 

We have seen that in a lossy dielectric the properties of the dielectric can 
be described using a complex permittivity and similarly in a lossy magnetic 

material use a complex permeability. 

Thus, for a lossy material the Maxwell equations can be written as 



0 , (V • ^ = 0) , 



V X E = — — = - 



dt 

VxH - J + e 



dB dH 



dt ' 



dE 
'dt ' 



where /x are complex quantities that characterize the material, and J is 
conduction current only, i.e. J = aE, 

Consider a plane wave propagating in the z direction 



E = ^oe^(^*"^^) . 



Then 



W X H = aE + iujsE = (cr + iijje)E 



Since for a plane wave propagating in the z direction the derivatives d/dx 
and d/dy of E and H are zero, we have 



0 



3 
0 



k 

dz 



(cr + iu)e)E 



y 
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Prom the LHS, we see that 

(V X = 0 {a + ioje)E^ = 0 . 

— ^ ^ 

Hence = 0 unless (cr + itue) = 0. Thus, E J- k. 

Since E^ = 0, we have that V • ^ = 0. 

Also 

^^^ = -^^ 

can be used to show that E J- H, 

In addition, we can find from the above equation that = 0. 

Since V • £^ = 0, the quantity e occurs only in the equation for V x 
and it is common to proceed as follows: 

^ dE ^ 
V X i7 = J + e— = aE + iujeE 

at 

W X H = iu; (^-^ +e^E = iu;(^e- E 

Now 

€ — t— = £ — le — t— = 6— lie H — ] = s 

which gives 

V X H = icjsE . 

Physically what has been done is to lump together the conduction current 

and the lossy dielectric current term. To an external observer they are in- 
separable. Only using some theory of the internal structure of the dielectric, 

they can be separated. 

We can summarize our findings of the propagation of the EM wave in a 
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conducting material that the electric and magnetic fields of the propagating 

wave satisfy the following equations 

V-E = 0 (33) 
V-^ = 0 (34) 

V X E = -A*^ = (35) 

V X # = iu)eE (36) 



where we have left the bar off e. 

We now proceed to solve Eqs. (33)-(36). We will see that the solution leads 
to a dispersive equation. We look for plane wave solutions, and will try to 

find how k behaves. 

The procedure is as follows: Taking Vx of (35) and using (36), we obtain 

V X (V X E) = u^pleE 
-\J^E = uj^y^eE 

Hence 



Since 

d^E 



= {jj^iisE 



= -k'E 



we finally obtain a dispersion equation 

7 2 2 

k — ijo lie 

This dispersion equation is not as simple as it looks. We cannot just say that 
phase velocity is 

LJ 1 



= - = 



k ^JJie 
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as e and are complex quantities and then k is complex. 
What does an imaginary value of k mean? 
We can write the complex number k as 



k — a — ij3 ^ 



where a and (5 are real. 

Then 



Clearly from this, the phase velocity is 



a 



and (5 — — Im(fc) is the attenuation coefficient. 

Physically the losses come from (1) Conduction currents, (2) Lossy dielectric, 
and (3) Lossy magnetic material, which are lumped together in e. 

We now find a and 

We can write 



k — uj (fie) 2 = oj 



k = 



CO 



— I 



e'-i[e" + - 

UJ 



a 



fj^ [e + —j +e fi 



it.' - ^/) 



Let 



p = e'fi' -fi" ie" + 



a 
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Then 



k — ij (p — iq)^ 



k = 



k = 



(J 



(J 



1 

'p'^ _^ 2 ^-ia.rctaii{q/p) 



1 
2 



where 6 = ^arctan q/p. 
Hence 



a = Re(fc) = cj (^p^ + 
P = — Im(fc) = cc; (^p^ + g 



J ^ cos^ 



1 

2\ 4 



sin6' 



which in general is quite comphcated. 



13.1 Wave Refraction and Skin Effect 

We will illustrate propagation of an EM wave in dielectrics and conducting 
materials on two examples: 

Example 1: The low loss dielectric and no magnetic effects 

For dielectrics, a is negligible and then we disregard magnetic properties. 

£" + -<£' 

1 

Hence 



k = 



(J 



CJ 



1 

2\ 2 



-i& 



1 

2 
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where 



(3 = arctan 



Then 



A; = - 



c \ So 



Hence 



a — — \ — cos 



c V ^0 \ 2e^ J c 
We can find phase velocity and refractive index 



Vp = — — C\ —f < c tor £ > £0 5 



n = — = J- > 1 , 
t^p V £0 

which are not affected by the losses. 

Thus, inside the dielectric, the EM wave will propagate with a phase ve- 
locity Vp < c and will be refracted. 

Similarly, we can find /? 



^ = -Im(fc) = -y- sinj^^^ 

Smce for small 0^ sin 6^ ^ 0^ we get 



cc; 6:" + a/cj 
c 2y^e^eo 
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Thus, losses (absorption of the wave) are small. 

We can summarize, that the theory predicts that the refractive index for 
a lossless dielectric is given by 

n = y/n^ . 

Table below compares theoretical values of n with that obtained experimen- 
tally. An excellent agreement is observed, except for polar molecules (e.g. 
water). 







n (experimental) 


Air 


1.00029 


1.00029 


Argon 


1.00028 


1.00028 


CO2 gas 


1.00047 


1.00045 


Benzene 


1.49 


1.48 


Ethanol 


5.3 


1.36 


NaCl 


2.47 


1.54 


Water 


9.0 


1.33 



Example 2: Good conductor 
A good conductor, e.g. Cu, Ag. 

(1) There are no dielectric and magnetic losses e^^ = /x'' = 0. 

Consider the Maxwell equation IV: 

OE 

V X H ^ J + e— =(TE + iujeE . (37) 
Definition of a Good Conductor 

Good conductor is when the conduction term on the rhs of the Maxwell 
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equation (37) dominates over the displacement current, i.e. when a ^ uoe. 

For example, for copper at a; = 1 MHz = 27r x 10^ rad 
cr 5.8 X 10^ 



ujE^ 27r X 10^ X 8.85 x IQ-i^ 



= 1.0 X 10 
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Consider general expression for k\ 



k 



— % 



t fl — fl [£ + 



UJ 



1 

2 



k ^ uj{ [e^^i! - 0] - i 



UJ 



1 

2 



k = 



1 

cr 1 2 



Now we might as well drop the dashes on £, understanding that they are 
real quantities, and obtain 



k — (J <efi — ifi— 



a I 2 



Cc? 



k 



a 



1 

2 



k = iJJ^J^l 



1 + 



.6:cc;. 



1 

2 



^-iarctan^ 

e ^'^^ 



Remembering that for a good conductor ^ ^ 1, we get 



k — UJ^J^l 



— e 2 



1 

2 



Since cos(7r/4) = sin(7r/4) = l/v^? we obtain 
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a = Re(fc) = ./^ = /3 = -Im(fc) 



Note a very heavy attenuation of this wave 

E = Eoe-^^ giM-a^) 

Since 

27r 



a — ^ — 



A 



the distance for attenuation "e" fold (i.e. amphtude falls to a factor ^ of its 
original value in a distance ^) is 



1 A 

z — — 



i.e. 




where 6 is called the skin depth in the conductor. It is the distance the wave 
must propagate in order to decay by an amount e"-*^. This effect is sometimes 
called "skin effect" as with an increasing a the current flows in a narrower 
and narrower layer, until in the limit of cr ^ oo a true current exists only on 
the surface of the conductor. 
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Questions: 

(1) Derive the special form of Poynting's Theorem apphcable in certain 
material media 

{E X H) ^ dS ^ - j H dV - j E dV - j E'JdV 

Interpret the above equation in terms of energy storage and energy flow etc. 
State qualitative meaning of each term in the equation. 

(2) Prove that the useful formula for the mean Poynting vector for sinu- 
soidal fields is 



where Ec = Eq exp^icjt) and = Hq exp(— icjt). 

(3) Show that the amplitude of a plane wave propagating in a non- 
conducting material is damped with the rate /5 which arises from the imagi- 
nary parts of the complex permittivity and permeability. 

(4) Show that in a good conductor an EM wave propagates on the surface 

of the conductor. 



Weekend exercises 
Find the skin depth for cooper at 60 MHz. 
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14 Transitions Across Boundaries for Elec- 
tromagnetic Fields 

14.1 Applications in dielectrics 

Across boundaries between different material media there are sharp changes 
in electrical properties cr. On a macroscopic scale the fields may have to 
be regarded as varying discontinuously across such boundaries. The source of 
such discontinuities will be the surface polarization charges and currents P^h 
and M X h discussed previously. 

We will use the Maxwell's divergence equations I and II to in- 
vestigate the transition of normal field components, and the curl 
equations III and IV to investigate the transition of tangential field 
components, 

14.1.1 Normal Component of B 

Consider the Maxwell equation II: V ♦ 5 = 0. 

Using the Gauss' divergence theorem, the Maxwell's equation can be writ- 
ten as 



where S is an arbitrary surface closing some area on the boundary plane. 




s 
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In order to evaluate the integral, we consider a thin cylinder of area 6S 
and thickness 5i including the boundary 

{B2'hi)6S+{Bi'h2)SS+ I B'hsdS = 0, 

sides 

Since B is finite everywhere and we are interested in the transformation of 
the field at the boundary, 6£ 0, the integral 

J B'hdS ^0 as S£ ^ 0 , 

sides 

Hence 

^2 • ni + 5i • ^2 = 0 . 

However 

^1 = -^2 = ^ 5 

and then 

n . (52 - Bi) = 0 or B2± = . 

Thus, the normal component of B is continuous across a boundary 
between two different dielectric materials. 




Figure 8: Polarization currents at the boundary between two different mate- 
rials. 

We can understand this result intuitively by noting that the B fields of 
polarization currents M xh are parallel to the boundary and so do not affect 
the normal component of 5, see Figure 8. 
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14,1.2 Normal Component of H 



Since B = /xiJ, we have 

Thus, when fii ^ fi2^ the normal component of H is not continuous across a 
boundary. This result for the normal component of H is the consequence of 
different magnetizations of the materials. 

14.1.3 Normal Component of E 

Since in nonconducting dielectrics V • 5 = 0, an identical argument to the 
above will show that the normal component of D is continuous across 
a boundary. 

In the case of dielectrics we write D = sE, 

Then eiEi± = S2E2± can be a useful boundary condition. 




The electric field of dielectric surface charge P-nis normal to the bound- 
ary. This causes the discontinuity in the E field. 

14.1.4 Tangential Component of E 

We will use the Maxwell's curl equations III and IV to investigate 
the transition of tangential field components. 
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For the tangential component of we will apply the Faraday induction 
law (Maxwell equation) to a closed path such as in the Figure, in which the 
sides perpendicular to the boundary are made infinitely short compared to 
the parallel sides L. 




Consider the Maxwell equation III: 



V X E = - 



dB 

~dt 



Integrating both sides of this equation over the surface a and applying the 
Stokes's theorem, we obtain 



/ E • dl = - I 



dB 
'dt 



ho da . 



Hence 



E2 • k + Ei-ii) L+ f E- dl = 

^ J ends 



dB 
~dt 



fiQ da , 



(38) 
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where ho is the unit vector normal to the surface a, and ii and £2 are unit 
vectors along the paths L on the side (1) and (2), respectively. 

As rfZ ^ 0 the rhs and the second term on Ihs of Eq. (38) go to zero, since 
E and B are finite everj^here. In this limit, the area enclosed by the path 
approaches zero. 

Since, ^2 = —h = t, we then have 

(^2 - ^1) • t = 0 . 

But ^ • t is the component of E tangential to the surface. Since this is true 
for any t, we obtain Ei\\ = £^211? and we conclude that: 

The tangential component of E is continuous across a boundary 
between two different dielectric materials. 

The continuity of a tangential component can be written in an equivalent 
form as 

n X (^2 - ^1) = 0 . 



A A 

n 




® nxE 



Explanation 

From the Figure above, we see that the cross product 



h X E 



= Esm{7r/2-9) ^EcosO , 



i.e. it is the tangential component of E. 
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14,1.5 Tangential Component of H 



Consider the Maxwell's equation IV: 

dE 



V X H = Jc + e 



dt ' 



where is the conduction current, i.e. not counting polarization currents. 

Stokes's theorem then gives 



/ 



H ' dl = I Jc ♦ nda + £ [ • fida . 

a J a ut 



T T 

Both terms on the rhs go to zero as SI 0 because Jc and dE/dt are finite. 
Hence 

(#2-^1) -^ = 0, 

or 

n X (^2 - ^1) = 0 . 

Thus, the tangential component of H is continuous across a bound- 
ary. 

14.1.6 Tangential component of B 

The tangential component of B is not in general continuous across a bound- 
ary because of the presence of the magnetic polarization surface currents 
M X n, which do not have a finite current density as they flow in an infinitely 
thin surface layer. Thus, if we examine the corresponding Maxwell equation 
for V X 5, the term in the integral involving J may stay finite as 5j 0. 



Since 

H = Soc^B - M , 

we have 

H2-Hi = eoc^ (B2 - Bi) - (M2 - Ml) . 
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Take scalar product of both sides with t, and using the fact that 

(h2 -HA-i = o, 



we obtain 



(^2 - Bi) 



t = 



1 



eoc 



(M2 - Ml) ■ i 



If Ml 7^ M2 the difference between the two surface currents generates a 
discontinuity in the B field (see Figure below). 



Ml 



® M,xni 



A 



Summary: Field components that are continuous across a bound- 



ary 



1. The normal component of D, 

2. The tangential component of E, 

3. The normal component of B, 



4. The tangential component of H, 
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Questions: 

(1) Prove the following boundary conditions at a bounding surface be- 
tween two dielectrics: 

(a) The normal component of B is continuous across the boundary. 

(b) The tangential component of E is continuous across the boundary. 

(c) The tangential component of H is continuous across the boundary. 
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15 Reflection and Transmission of Waves 
Across a Boundary 



At a boundary between two different media, the general boundary condi- 
tions cannot be satisfied by a transmitted wave alone. There has to be a 
reflected wave also. What happens is that the incident radiation is absorbed 
by charges in the boundary and reradiated in all directions. The waves in- 
terfere destructively except in two directions along those of the reflected and 
transmitted waves. The directions, amplitudes and phases of the reflected 
and transmitted waves can be derived from the general boundary conditions 
already obtained. We consider linear isotropic materials characterized by 
electric and magnetic constants and cr. 

15.1 Representation of Plane Waves in Different Di- 
rections 




\ 

\ 

\ 

\ 

\ 

\ 

Let hp is the unit normal to the phasefront (wavefront) of a propagating 
wave, and r is a position vector that is independent of hp. The vector f 
determines position of a wave-front as seen by an observer O. 



199 



Suppose that the wave propagates in the z direction. Then 



E — Eq exp[i (cut — kz)] . 



Since 



z — n ' r ^ 



we obtain 



— * — * 

E — Eq exp[i {ijjt — h ' r k)\ . 



Thus, the observer can distinguish different waves at the point A by differ- 
ent n's. 



15.1.1 Representation of B in terms of E 



For the field of a plane wave, we have 

^ k 

B = —hp X E , 

where hp is the unit vector in the direction of propagation of the field. 
Recall that the relation (39) arises from the Maxwell equation III: 



(39) 



V X £ = - 



dB 
~dt 



To show this, we expand both sides of this equation, and obtain 



t 

0 



J 

0 



Ex E. 



y 



k 

d_ 

dz 

0 



dB ^.dB 



X 



^.dB 



y 



dt 



dt 



where we have used the fact that the wave propagates in the z direction. 

Comparing the coefficients standing at the same unit vectors, we find that 
the X component gives 



dE 



y 



dz 



dB, 
~dt 
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from which, we obtain 



ikE„ — —iuB^ or Ey cj 



Bx k 



The y component gives: 



dE, ^ dBy 

dz dt 



from which, we obtain 



—ikEx = —iuBy or 



Ex cj 



By k 



Thus 



/C V ^ k 

Since E^ B^ hp^ are mutually orthogonal, E x B gives the direction of hp 
(Poynting result) of the propagation direction. Then hp x E gives the direc- 
tion of B, Hence 

^ B k 

H = — = — h^ X E . (40) 

In the next few sections, we will use the continuity conditions for E and H. 
With the relation (40), we will be able to limit the analysis to the electric 
field alone, as knowing properties of E^ we can find from Eq. (40) properties 
of H, 



15.2 Directions of Reflected and Transmitted Waves 

Here, we will show that some familiar elementary results on reflection and 
transmission can be derived from the Maxwell equations. 

The most useful results in this connection concern the continuity of tan- 
gential components of the E and H fields. However, most of the results can 
be derived from the fact that the tangential component of E is continuous 
across a boundary. 
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Consider properties of an EM wave at a boundary between two materials. 
Let the origin of position coordinates f be located on the boundary S, 
The electric field for incident, reflected and transmitted waves are 

Eq exp [i {ujt — hi ^ f ki)] , 
El exp [i [uot — hr ' f ki)] , 
E2 exp [i {cut — fit ' f ^2)] . 

These equations determine the three electric fields relative to the direction 
of observation r. 

Why do we need the presence of the reflected wave in the propagation of 
an incident wave through the boundary? 

The answer to this question is provided by the requirement that the tan- 
gential components of E and H must be continuous through the boundary. 
Suppose that E is polarized in the plane of incidence. Then 

Ei cos 6i — Ej. cos Or = Et cos Of , (41) 

and 

Hi + Hr = Ht, (42) 

Since 

k ^ 
H = — E , 
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the relation (42) takes the form 



Ei + Er 



E: 



(43) 



ki 



t 5 



as for a dielectric /xi = fi2 — jJ^o- 

Now we see that without Er^ Eq. (41) gives 



Ei 



cos 6t 



Et 



COS 6i 



while Eq. (43) gives 



Thus, without we would get two different values for Ei or Et^ which we 
cannot accept as both continuity conditions Eqs. (41) and (42) must be satis- 
fied at the same moment. Hence, we conclude that the continuity conditions 
for E and H will be satisfied only if Er / 0. 

What are the relative directions of hi^hj. and n^? 

From the boundary condition for the tangential component of E: 



This relation must hold over the whole surface S for all r (subject to fi-f^ 0). 
Thus, the exponential phase factors must all be the same. Otherwise if it 
was true for one r it would not be true for other f 's, but we have a freedom 
of choosing r. Hence 




I.e. 



h X IEq exp [i (cut — hi - f ki)] + Ei exp [i (cut — hr - r ki) 




kihi ' f 



kihr • r — k2ht - r 
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from which we have 

hi 'f — fir 'f ^ 

and 

ki 

These relations will help us to prove that 

Incident, reflected and transmitted rays are coplanar 

i.e. hi^ hr^ ht are coplanar. 

To show this, we use the relation 

f — —h X (h X r) 

Proof: 

Since 

h X (h X f) — (h ^ r)h — {n » n)f 
and n • r = 0 as the vector r lies on the plane we obtain 

fix {fix r) — — f 

as required. 

Hence 

• [n X (n X f^] = • [n X (n X r)] . 
Interchange (•) and (x) products 

[hi X h) ^ [h X r) — [hr x h) - [h x r) . 
This must be true for all f in plane S, Thus: 



Hi X n = X n 
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A A 

rij^xn 



This implies that hr is in the "plane of incidence", i.e. the plane containing h 
and hi. 

Similarly 

kihi • r = k2ht • r 

implies that 

ki hi ^ h X (h X = k2 ht ' h X (h X f) . 
Interchange (•) and (x) products, we find that 

ki {hi X h) ' {h X r) — ^2 {ht x h) - {fix r) 
for all r m S, Thus 

kihi xh — k2ht X h , 
which implies that n^, n, and ht are coplanar. 

The coplanar property of the waves is observed in any experiment. This 
is an another example of a remarkable triumph of the Maxwell's electromag- 
netic theory 



205 



15.3 Snell's Law of Refraction and Angle of Reflection 

Since 

hi X h = hr X h 

we have that 

sin 9i — sin Oj. 

Thus 

Oi — 9j. . 

The angle of incidence equals the angle of reflection, another familiar law of 

elementary optics. 

Moreover 

ki hi X h — k2 ht X h 

from which we have 

ki sin 6i — ^2 sin 6t 

and then 

sin Oi ^2 
sm ut ki 

This is the well-known law of refraction in optics, called the SnelFs law. 

In the case where ki and ^2 are purely real (e.g. in dielectrics), the 
refractive index has a simple physical interpretation 

k2 uj/ki Vi 
ki uJlk^ V2 

i.e. the refractive index in equal to the ratio of phase velocities. 
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Questions: 

(1) Prove that in the reflection and refraction at a bounding surface, the 
direction of incident, reflected and refracted waves are coplanar. 

(2) Using the results of Question 1, derive the famiUar laws of elementary 

optics: 

(a) Angle of reflection equals to the angle of incidence. 

(b) Snell's law of refraction. 

(3) Show, using the continuity conditions for E and H that both reflection 
and refraction takes place in the incidence of light on a boundary between 
two dielectrics. 
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16 FresnePs Equations 



The boundary condition on tangential E does not give sufficient information 

to calculate E^ and Et in terms of Ei, For a given Ei there are still two 

— ♦ — * — * 

unknowns in the equation for continuity of tangential E viz Ej. and Et. We 
need a second relation between Ei^ E^ and Et. This can be obtained from 
the continuity of the magnetic field. 

We know that B tangential is not continuous across a boundary because 
of the presence of M x n surface currents in magnetized materials. To allow 
for such possibilities we can use the more general condition that tangential 
component of H is continuous, i.e. 

n X {Hi + Hj}j = h X Ht . 

Since 

ff = — = — rzp X E , 

where hp is a unit ray vector in the direction of propagation, we have an 
equation 

n X I rii X Ei-{ Ur x Er 




= n X -^ht xEt) , (44) 



which together with 

hx (^Ei + Er) =hx Et , (45) 
contain sufficient information to determine Ej. and Et in terms of Ei. 

The solution of these two simultaneous vector equations is greatly facili- 
tated by considering separately, electric field components parallel and normal 
to the plane of incidence. 

Equations (44) and (45) also provide a simple explanation of why we need 
reflected field at the boundary to obtain the correct results for the field am- 
plitudes. 
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16.1 Ei normal to plane of incidence 



In this case, the incident electric field Ei is purely tangential to the boundary. 
Since the materials are isotropic, the induced fields Er and Et will also be 
tangential to the the boundary. 




Thus the condition 



n X (^Ei + Er^ —h X Et 



gives 



Eq -\- El — E2 . 



Note then that 



h ' Eq = h ' El = h ' E2 = 0 . 



The other equation for Eq^ Ei and E2 comes from the H condition through 
the relation 

n X I rii X Ei -\ rir x Er 



ijfii 



UJfli 



h 



= n x 



fit X Et 



and since the phase factors in E^, exp i{(jjt — ho^-fk) are the same, we obtain 

ki 



Ufli 



n X {714 X 



Eq^ + h X (hr X Ei^ 



nx (htX E2j . 



209 



Using the relation A x {B x C) = {A - C)B - {A • B)C and the fact that 

h ' Eq = fi ^ El = h ^ E2 = we obtain 

n ' TiiEo + n • n^Ei ) = n • ntE2 



However 

fi ' fii = cos(7r — 9i) — — cosOi 

h ' hr = COS ^7. = COS 

n • = cos(7r — 9t) — — cos 9t 

and then, we obtain 



^"^ ^Eq cos 6*^ — ^1 cos 9i) — -^^2 cos 6^^ 



Since ^'s are all in the same direction, we might as well drop the vector 
signs. Thus 

Eo + Ei = E2 (46) 

Eq cos 6i — El cos 6i = E2 cos 6t . (47) 

kifi2 

Eliminating Ei using (46), Ei — E2 — Eq^ we get 

^2/^1 

£"0 cos 6i — (E2 — Eq) cos 6i — E2 cos 9t 

fclM2 



2Eo cos = I cos 9i + -^-^ cos | £"2 

fcl/i2 



2kifi2EQ cos = (A;i/i2 cos 6^^ + A;2/ii cos 9t) E2 



and then 



2kifi2COs9i 
/ci/i2 COS 6^^ + A;2/ii cos 9t 
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Since 



ki sin Oi = sin 6i 



we can eliminate 6t 



COS 6t = \Jk2 — ki sin^ 6^^ 



and obtain 



_ 2fci/i2COS^; _ , 

£^2 = ^ ^ = ^0 . (48) 

A;i/i2 cos + /ii — fc^ sin^ 6^^ 

Similarly, eliminating E2 = Eq + Ei from (47) above, we obtain 

fci/i2 cos Oi- fjiiJk^- kl sin^ (9^ 
/ci/i2 COS Oi + fjbiJ k2 — kl sin 6^^ 

Equations (48) and (49) are called Fresnel equations for the electric field 
amplitudes. 

The corresponding H fields are not parallel to each other, but their rela- 
tive magnitudes can be deducted from equations of the form 

^ k ^ kE 
H — — n X E i.e. H — — . 

16.2 in the plane of incidence 

In this case Hi is tangential to the boundary plane and then Hj. and Ht are 
tangential too. Thus 

h X (^Hi + Hj}j = hx Ht 

becomes 

Hq + Hi = H2 . 
The continuity of tangential E is given by 

n X (Eq + El) = n X ^2 . 
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It is simpler to work in terms of H in this case. 
We express E in terms of H 

E = — — h X B = — — h X H 
k k 

Thus 



n X 



-—rii X Ho + -—rir x Hi 



By a procedure similar to the case of Ei normal to plane of incidence, we 
obtain 



k^Uo cos 9i — Uikix kn — A;? sin^ 9i 
^2/^1 COS Oi + /X2^i Y ^2 ~ ^1 sin^ 6i 

H, = 2fc|^^^cos^__ ^ ^^^^ 

fcl/ii COS + /i2fci Y ^2 ~ ^1 si^^ ^« 

Equations (50) and (51) are called Fresnel equations for the magnetic 
field amplitudes. 

This time the E^s are not parallel, but their relative amplitudes can be 
deducted from E — {fiu;/k)H, 



16.3 Fresnel Equations for dielectric media 



In a dielectric: conductivities ai = cr2 = 0, /i-i = /i2 = /^05 ^ = 27r/A =real, and 

u; 1 



Vp = - = 



k y^ejl^ 

Since ~ 1/A ~ 1/vp ~ ^/e^ the Snell's law {ki sin 9i — ^2 sin^t) becomes 



^/el sin 9i = sin 9t 
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and then 

sin 9i _ fe^ _ vi 
sin 6t V ^1 ^2 



Consider two examples: 

(1) £^ normal to the plane of incidence. In general 

Eq kifi2 cos 9i + fc2/ii cos 9t 



Since 



we obtain 



Hence 



- ^2 sin 9i 
Ml = /^2 = /^o and — = - , 

ki sm Ut 



E2 2 cos 9i 2 cos 9i 

Eo ~ cos 9, + ^ cos 9t ~ cos 9, + ^ 

2 COS 6i sin 6t 

cos 6i sin + sin 6i cos 

£^9 2 cos 6i sin 



Eo sin(^i + ^i) • 



Similarly 



Eo sin {dt + di) ' 



(2) For E in the plane of incidence 



E2 2 cos di cos 



Eo sin {0t + ^i) cos - dt) ' 

Eo tan {Ot + 9^) ' 
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17 Applications of the Boundary Conditions 
and the Fresnel Equations 



Let us now examine some of the consequences of the Snell's law. There 
are two cases possible: n2 > rii and ni > 712. In the first case, an optical 
wave travels from an optically "rarer" to optically "denser" medium. In the 
second case, we have the inverse situation. We will consider these two cases 
separately for dielectrics and conductors. 

17.1 Applications in dielectrics 

17.1.1 Polarization by reflection 



\ E 










\ ^ £2 



For the case of E in the plane of incidence the ratio of reflected to incident 
amplitude is 

El ^ tan(g, - 9t) 
Eo ~ tan(^i + Ot) * 

If 6^^ + = I then tan(^i + 9t) = 00 and consequently Ei — 0. 

At the same time E\ will not be zero for the electric field component 
normal to the plane of incidence. 
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Thus if Ei has arbitrary polarization then Er will be plane polarized with Ej 

normal to the plane of incidence. 

If ^. + = I then = f - Thus 



sin 9i _ ^ _ _ sin Oj _ ^ 
sinOt \j Si sin(| - 9i) 

Hence, the angle of incidence for total linear polarization of the reflected 
wave is 



9i = arctan J— . 

V ^1 

It is known in the literature as the Brewster's angle. 

An alternative proof that if E is in the plane of incidence, then Ei = 0 
(no reflected field polarized in the plane of incidence). 

We can use the continuity conditions for the tangential components of E 

and from which we have 

Eo cos 9- El cos 9 = E2 cos 9t , (52) 

and 

Ho + Hi = H2, (53) 

Since 

H=—E, 

OJfl 

and = I — ^, we get 

Eo- El = E2t8in9 , 

Eo + Ei = ni2^2 , (54) 

where ni2 = ^2/^1. 

However, from the Snell's law we have that 

sin 6^ sin (9 sin6' 



sin 6^ 



t sin (1 - 9) cos 9 



= tan 9 = ni2 . 



Thus, Eqs. (54) will be satisfied simultaneously only if Ei = 0, i.e. there is 
no reflected field in the plane of incidence. 
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17.1.2 Total internal reflection 



Consider the SnelFs law 



sin^f = 



sin 9^ 

^21 



^2 



where n2i — \l — . 

^1 



In the case < 1 (62 < ^i)? goii^g from an optically more dense to optically 
less dense medium e.g. from water in to air, real angles 9t are obtained only 
for sin Oi < 7121 . 





Z COSBf 

/ X sin9, 

/ / X 








X -(x, z) 


z, 





For greater sinOt > 1, and then the angle of refraction 6t becomes 
imaginary. In this case, there is no real refracted wave, only a reflected wave. 
The Fresnel equations are still capable of giving a formal solution in this case 



cos 9t — V 1 — sin^ 9t — 



1 - 



sin^ 9a 



nil 



= I 



sin^ 9^ 



where j3 is real positive. Put, e.g. 



nil 



- 1 = , 



Sin Ut = = a . 

n2i 
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Now consider the propagation of the transmitted wave in the less optically 
dense medium, for which 

with 

fit ' r — z cos 6t + X sinOt = —i^z + ax , 

and then 

Here = kp and a' = fca, i.e. there is attenuation in the z direction but no 

phase propagation. Phase propagation occurs in the x direction along the 
boundary. 

This illustrates a general method of applying the Presnel equations. For 
only a limited range of circumstances will all the angles 6^^, 6^^, Of be real. We 
can however always apply a generalized SnelVs Law ^2 sin^t = ki smOi to 
find sin 9t and cos 6t and proceed as above. 

Note that the planes of constant phase are normal to the boundary (i.e. they 

have their normals tangential to the boundary). The phase of the trans- 
mitted wave below the boundary must match the incident wave above. The 

wavelength in the second medium is 

^ 27V 27V 27V Xon2i 
~~o^ ~ka~ ~ sin^i ' 

Ao n2i 

where Ao is the wavelength of freely propagating waves in this medium. 

The planes of constant amplitude in the transmitted medium are parallel to 
the boundary. 



17.2 Transmission and Reflection at a Conducting 
Surface 

Propagation of EM fields in conductors (metals) is more complicated phe- 
nomenon than in dielectrics. Consider a propagation of an EM wave in a 



217 




medium composed of a dielectric and a conductor, and assume that the in- 
cident wave originates in the dielectric. 

From previous work, we know that in the dielectric 

and in the conductor 



Moreover, from the Snell's law 

sin 9t ^ ^ sin 6i . 

k2 

Assume a good conductor. Since ^2 ^ 00 as ^ 1, we see that sin6'f 0 
independent of 9i, Thus, the direction of the transmitted wave is normal to 
the surface independent of the angle of incidence 9i. 

It follows that in the conductor the field vectors E and H lie tangential to the 

boundary and so the normal components of these vectors on the conductor 
side of the boundary are zero. 

It follows that: 

• Since the normal component of B (or H) is continuous across the 
boundary, the normal component of if or 5 is zero on the dielectric 
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side also. 

Thus the normal component of B of the reflected wave must be equal 
and opposite to that of the incident wave. 



In the conductor 



This happens because 



= K 



1 - 



SCO 



oo 



cr 

as — > 1 



This means that the electric field in the conductor (which is tangential 
to the surface) 0. 

Since E\\ is continuous across the boundary E\\ is zero also in the di- 
electric at the boundary. 

In summary, we have two useful special boundary conditions at the 

surface between a dielectric and a perfect conductor: 

1. The tangential component of E = 0, 



2. The normal component of 5 or = 0< 



17.2,1 Field vectors at normal incidence 

We now consider the special case of normal incidence at a boundary, i.e. 
when the wave propagation vector coincides with the normal to the boundary. 
Figure 9 shows how the field vectors must look in the plane of the surface 
between dielectric and conductor. 
Since for a good conductor 

a 

^ 00 , 

the field components 

Et — Ei — Er ^ 0 , 
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Hi A A 


E, 


E, 

H. 



Figure 9: The field vectors in the plane of the surface between dielectric and 
conductor. 

Ht — Hi-\- Hr 2Hi , 
In this case, the power reflection coefficient 



One could think that under the normal incidence, there is only the incident 

and transmitted wave with no reflected wave. Here, we prove the necessity 
of assuming the existence of the reflected wave. 

From the continuity of the tangential components, we have 

Ei-Er = Et, (55) 
Hi + Hr = Ht. (56) 

However 

H=,l^E = ,^,^E = nME. 
V /^o V ^0 V ^0 V 

Thus, Eq. (56) can be written as 

riiEi + niEr = n2Et , 
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and then we get two equations for the ampHtudes of the electric field 



Ei — Er 

Ei + Ej. 



Et 




Et. 



If Ej. is missing, we could not simultaneously satisfy both equations. Thus, 
there always is a reflected wave in the normal incidence. 

17.3 Wave Propagation in a Conducting Rectangular 
Pipe (The Rectangular Waveguide) 

In this lecture, we discuss the propagation of bounded EM waves by consid- 
ering the propagation of radiation through a waveguide where the radiation 
is fully confined in the transverse plane. We will consider the case where the 
bounding walls are planar and cross section is rectangular. 



From the nature of the problem we see that we are likely to get standing 
waves in the x direction due to reflections between the walls x — {) and rr = a, 
standing waves in the y direction due to reflections between the walls y = 0 
and y — h and waves propagating in the z direction. 

We solve Maxwell's equations subject to the good conductor boundary 
conditions being satisfied at a; = 0, x = a, ^/ = 0 and y — h. We write the z 
dependence of any field component in the form 




Thus 



d 



e 



dz 
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where 7 describes the propagation conditions, e.g. 7 purely imaginary de- 
scribes a wave propagating without loss. 

We describe the electromagnetic field by the vector pair and we use 

Maxwell's equations in the form 

V- = - = 0 

€ 

V'B = OorV-^ = 0 

_ dB OH 

^ dE - dE 
Vxif = J + e-=.E + e- 

Certain characteristic modes of propagation are found. In general a 
characteristic mode is one which propagates with constant polarization. 

In a rectangular waveguide, one may have TE (transverse electric) waves, 

or TM (transverse magnetic) wave, but not TEM wave. If both E and B 
fields are transverse, the wave would be going straight down the guide. How- 
ever, such a wave would not satisfy various boundary conditions. 



17.3.1 Transverse Electric (TE) Modes 



It is possible to propagate a wave with the electric field polarized in the xy 
plane by lifting the restriction that both the electric and magnetic fields are 

transverse. 

Look for solutions to Maxwell's equations with E^ = 0 that satisfy the 

good conductor boundary conditions. 

For the TE modes we use Maxwell's equations to derive a wave 
equation for H^^ the longitudinal component of H, This equation 
is used to satisfy the boundary conditions. 

Write out Maxwell's equations in Cartesian coordinates. 



Prom V X E + fi^ = 0 



dt 



A. 

dx 



d_ 

dy 



k 



•7 



+ iu>fj,H = 0 



E^ K. 0 



y 
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X component 
y component 

z component 



^Ey + iu^Hj, = 0 

- ^E^ + ioJiiHy = 0 
dE.. dE^ ^ . ^ 



'y 



dx 



dy 



Prom V • = 0 



dH^ dHy dH, 
dx dy dz 



= 0 



Hence: 



dx 



+ 



y 



dy 



- 7^. = 0 



From V X H - {aE + = 0 



_d_ 

dx 



J 

d_ 

dy 



k 



■1 



H., -ff. 



- {a+ju)e)E=- 0 



m 

dy 



+ ^Hy - (cr + iuje)E^ ^ 0 



—^Hx — — (cr + i(je)Ey = 0 



5a; 
5i7 



= 0 (E, = 0) 



From V ♦ ^ = p/e = 0 

And using the TE condition E^; = 0 



dEx dE 



dx 



rom Eqs. (57) and (58) 



E. 



X 



E: 



y 



y 



X 



1 



= Z 
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where Z is called the wave impedance. 

In Eq. (59), we can then express and Ey in terms of Hy and and show 
that it becomes identical to Eq. (60). 

Use Eq. (65) in Eqs. (61) and (62) 
In Eq. (61) substitute for Er^\ 

+ ^Hy — (cr + lijje) Ey 0 



% 7 



Hence 



-1 dE 

li-i, — 



" 7-(a + ^a;£)^ dy 

7 dH, 



7^ — i(jfi{a + icjs) dy 
In Eq. (62) substitute for Ey and similarly get 



(66) 



H,= , . -\. '-^ (67) 

Using Eqs. (65), (66), and (67), we find that Eqs. (63) and (64) are automat- 
ically satisfied. 

Putting Eqs. (66) and (67) into Eq. (60), we obtain 

-7 d^H, 



7^ — iojfi^a + ioje) dx^ 

7 d^H, 



7^ — iu^{a + iue) dy'' 
Dividing by —7 and multiplying by fc^ = 7^ — iu>fj,{a + ius), we obtain 

If we can solve Eq. (68) for then we can solve Eqs. (66) and (67) for Hj, 
and Hy, Then we can solve Eq. (65) for E^ and Ey^ and then we will know 
all the field components. 



224 



17,3.2 Boundary Conditions 



The boundary conditions at the surface of a good conductor are 

1. H normal to boundary (in xy plane) = 0. 

2. E tangential to boundary (in xy plane) = 0. 




Thus, we must have at the boundaries x = 0 and x = — 0 and 

Ey = 0. 

Looking at Eq. (62), we see that this means that dHzjdx must equal zero at 

X — ^ and X — a. 

Similarly, at the boundaries y = 0 and y — h we must have Hy — 0 and 
E^ = 0. 

Looking at Eq. (61), we see that this means that dH^/dy must equal zero at 
y = 0 and y — h. 

The solution of Eq. (68), which satisfies these boundary conditions is of 
the form 

Hz = i^o cos{kj.x) cos{kyy) 

with 

k^a — mir and L,b = nir 
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= Ho cos cos (^) e^-- . 

Here m, n = 0, 1, 2 An m, n combination represents a possible TE mode 

of propagation. The modes are designated in the form TE, 
Then Eq. (68) gives 



b 

For a non-trivial solution, ^ 0, and then 



(-^1 +1"^) =7 + 



2 /m7r\2 /n7r\2 ^ 
7 - (~^) + l"^] + ^/^^(^^ + ^^^) 

In general 7 is complex. We can write 7 = /5 + ia. 
Then we can write 

= ifo cos (^) cos (^) e-^^ e^^— ) 

Then a — 27r/Xg^ where Xg is the wavelength in the waveguide at frequency oj. 
The phase velocity in the waveguide is Vp = cj/a. 

The parameter /? is the attenuation coefficient describing losses in the 
waveguide. Energy loss may be due to: 

• Ohmic resistivity (cr of the medium finite). 

• Dielectric losses (described by £ having an imaginary component) 

• Magnetic losses (described by fi having an imaginary component) 

17.3.3 TE Modes in a Lossless Waveguide 

Assume a lossless propagation for which we have a = 0 and e and fi both 
purely real. In this case 

^ = -e^' + ( -) + 
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or introducing the phase velocity, we obtain 



2 uJ^ ( m-K^^ ( nix\ ^ 



7 = 2+ + 

\ a J 



b 



where vq is the phase velocity of propagation of waves in an infinite 
(unbounded) medium of the type filling the waveguide, 

(J _ 27r/ _ 27r 

where Aq is the infinite medium wavelength at that frequency. 



2 



The parameter 7^ can be negative or positive, and the nature of the 

propagation changes according as: 

1. Consider 7^ negative. Then 7 = ikg is purely imaginary and we have a 
propagating wave 



= H, cos (^) cos e^(-^^^) 



where kg is the guide propagation constant. 

The wave then propagates with a guide wavelength A^^ given by 



27r , lf27T\^ frmr\'^ (rni^^ 

— kn — 



Ag ^ VVAo^ \ a J \ b 



'9 



a) 



There is a maximum Aq (=Ac say) such that kg is real 



1 /m\^ f n^^ 

+ 



\l \2aJ \2b 

Thus, there is a minimum frequency fmn such that the TEmn mode will 
propagate down the waveguide 



fmn ^-V(2aj +V26. 
This can be derived from the cut-off condition fmn = ^m- 
Thus, the waveguide acts as a high- pass filter for any m, n mode. 
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2. Consider 7^ positive. Then 7 = ^ say is purely real and 



fm7rx\ 

Hz — 120 cos 

V a 



) 



n7ry\ -^^ 
cos I I e ^ e 



There is no phase propagation but amplitude attenuation. Note there 
is no energy loss mechanism available. This is an evanescent mode at 
frequencies less that fmn analogous to the case of total internal reflec- 
tion. 



Field components in the TE modes. 



Since 



m7rx\ 

^ ^ a J \ b 



- cos ( 1 cos ( ^ ] e^"*-^^ 



From Eq. (66) 



dHz\ / m7rx\ 

j.eos^— j 



sm —-- e ^ 



From Eq. (67) 



H 



X 



m 

dx 



. /m7rx\ 



~ sm 



COS I I e ^ 



And from Eq. (65) 



frmrx\ , f n7ry\ -^^ 
Ej, Hy COS [—^j sm I — ^ ) e ^ 



~ i?^ ~ sin 



. / m7rx\ 



'y ^^x 



cos I —7— I e 



Allowed values of m and n for TE modes 
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m = n = 0 is not allowed since E^^ Ey^ Hy all contain sine terms so 
all fields vanish for m = n = 0. All other TEmn modes are allowed. 

m = n = 0 mode is never possible in a transmission line consisting of 
a single closed conductor like a rectangular waveguide. It is possible in 2- 
conductor lines e.g. the coaxial line or the twin wire transmission line. 

TM (transverse magnetic) modes 

Put Hz — Q and go through the whole procedure again. E^ ^ 0 now. 

Equations analogous to Eqs. (66), (67), and (68) appear for components of 
the E vector this time. Consequently the previous discussion about modes 
and their cut-off frequencies for TE modes is also true for TM modes. The 
only difference is that more modes are not allowed. 

TM modes with either m = 0 OR n = 0 are NOT ALLOWED. 

Special properties of the TEiq mode 

If the transverse dimensions of the rectangular waveguide are different 
(a 7^ b) there is a finite range of frequencies over which the TEio mode is the 
only allowed mode. This means that a waveguide can be designed to allow 
propagation in one mode only. We have 




where we adopt the convention that a > b. In the frequency range /lo foi 

the TEio mode is the only mode allowed. 
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17-3.4 Phase and Group Velocities for TE and TM Modes 



The phase velocity 




Thus 

^0 




where Vq = cj/ko is the phase velocity in the unbounded medium and fco = 
27r/Ao. 

We see that the phase velocity of the wave inside the waveguide is greater 
than in an unbounded medium, and so may be greater than the speed of 
light in vacuum. 

Alternately we could write 

/m7r\2 /n7r\2 /27r\^ ,9 

(it) + (t) = (a:) = 
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where Ac is the infinite medium wavelength at the cut-off frequency for 
the m, n mode. 



LJ UJ Vq 



*=o - (fe)^ f^) 




provided f > fc i-e. in the pass-band for that mode. 

As / ^ the phase velocity Vp ^ oo. 

If a vacuum or air fills the waveguide then Vq = c and Vp > c. 



The group velocity 



Since Vp > c typically we see that Vp is not the velocity of propagation 
of energy or information down the waveguide. The waveguide is a dispersive 
medium. The group velocity Vg — is then different from the phase veloc- 
ity Vp = Lu/kg, 



The group velocity is the velocity of propagation of some modulation of 
the wave that carries information. A single frequency harmonic wave carries 
no information. It is just there. A finite bandwidth is required to carry in 
formation. Consider the sum of two frequencies cu and u; -\- duj 

cos{(jt — kz) + cos[(cc; + dijj)t — (fc + dk)z\ 



— 2 cos 



= 2 cos 





dk 






w- 


2 




dk 






< 2 


2 



COS 



200 + duj 2k + dk 



t- 



z 



cos(cc?t — kz) , 



where we have assumed that dto <C 2(jJ and dk <C 2k, 
The velocity of the amplitude modulation is 



duj 



f dk 
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In the rectangular waveguide 



k 



g 



2lT 



^0 "~ 



u 



m 



where ko = cj/vm and kc = oJc/'^m- 

Remember is the infinite medium phase velocity and uj^ is the cut-off 
(angular) frequency of the m, n mode. Then differentiating 



dkg 
du 



1 1 



(J 



Finally 



dcu 



^9 = 



dk 



9 



(J 



Thus 



7? 7? = 7? 



V 



m 



= V 



m 



In a vacuum-filled waveguide VgVp — c^. Thus, relativity is still all right. 



Questions: 

(1) Show that under the Brewster's angle of incidence there is no reflected 
electric field in the plane of the incidence. 

(2) Show that in a vacuum-filled rectangular waveguide Vp > c 
and VgVp — c^. 
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18 Relativistic Transformation of the Elec 
tromagnetic Field 



18.1 The Principle of Relativity 

1. The laws of physics are the same in all inertial reference frames. 

2. The speed of light in vacuum is independent of the uniform motion of 

the observer or source. 

The constancy of the velocity of light, independent of the motion of the 
source, gives rise to the relations between space and time coordinates in 
different inertial reference frames known as Lorentz transformations. 

Consider a stationary reference frame S and a inertial frame moving 
with a velocity u parallel to the x axis. 

The time and space coordinates in are related to those in S by the 
Lorentz transformations 

7(x - Pet) , 

y , 

where 7 = (1 — P^)~^'^ is the Lorentz factor, and P — u/c. 



X — 

y' = 

z' = 

ct' = 



The above transformation corresponds to a situation of parallel to the x 
axis. If the axis in S and remain parallel, but the velocity u of the frame 
is in an arbitrary direction, the generalization of the above transformations is 

r' = r+{^-l)^^-^Pct 
ct' = ^(ct-P-A , 

where (5 = u/c 
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Proof: 



Decompose the vector r into two components: 
parallel and normal to P 

r = rji + r j_ 

Then, using the one dimensional Lorentz transformations, we have 



rjj = 7 (^rji — Pctj 
fj_ = f± 

We can write the parallel and normal components as 

^ ^ (r-/3)/3 

f± = r — fii 



Hence 



Similarly 



However 



which gives 



as required. 



= 7 — ^ct^ + r — fji 
= r + (7-l)^^l^-7/3ct 



(i' = 7 (c^ - /5a;) = 7 (ct - /Jry) 



ct' = -^(ct-P-f) 
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The inverse Lorentz transformation is 



{r'_f3)f3 
ct = ^{ct' + i5-f') 



The principle of relativity indicates that the Maxwell equations and the 
continuity equation should be invariant under the Lorentz transformation. 



In the frame S: 



V-D = p 
VB = 0 

V X E = - 



Vxif = J + 



dt 
dD 



VJ = - 



dt 
dp 

'dt 



In the frame S': 



V-D' = p' 

VB' = 0 

^/ ^, dB' 

V xE' = - 



V xH' = J' + 



dt' 

dD' 



V-J' = - 



dt' 

d^ 
dt' 



where the prime variables are functions of the transformed variables 

[t' and r'). 

First, we will illustrate the transformation of an arbitrary scalar function 
^ and the divergence V • F. 



235 



Consider the transformation of a time derivative /d{ct) 

dct' dx' dy' 



d{ct) d{ct')d{ct) dx'd{ct) dy' d{ct) 
dm dz' dm Oct' dr' 



dz'd{ct) d{ct')d{ct) d{ct) 



However 



dct' dr' -> 

7 , TTTTV = -IP 



d{ct) ' ' d{ct) 



which gives 



f d 



= 7 



d{ct) ' \d{ct 
Consider now the divergence 



— B-V] ^ 

d{ct') ^ ^ )^ 



_ - dF^ dFy dF^ 
V ■ F = — - + — - H 

dx dy dz 



Since 



dF^ dF^ dct' dF^dx' 
dx d{cf ) dx dx^ dx 

dy^ dx dz^ dx 



and 



dct^ 
dx 



dy' dz' 
— = — = 0 
dx dx 



we obtain 



dx ^ d{ct^) ^ dx^ 
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where 



a, = 1 + (7 - 



Similarly, for Fy and Fz^ and finally we get 



V • F = a * (V' • F) - 7/3 



where a is a 3 x 3 diagonal matrix 



a — 



l + (7- 

0 
0 



1) 



/3? 



a- 



0 



l + (7 



0 



dF 

d{ct') 



1 '\ 



0 
0 



i + (7-i)i 



Using the above transformations, we can derive transformations for the 
current density J and the charge density p. 

In order to do it, we consider the continuity equation, that can be written as 



Hence 



Since 



V- J = - 



dcp 

d{ct) 



a*{V -J)- 7/3 



dJ 



d{ct') 



= -1 



d 



d{ct 

d 



V-{a*J)-^0-V{cp) = (^cp-P-j) 



P-V{cp)^V-{cpP) 



we obtain 



V' • (a* J — "ycpi 



d 

d{ct') L 



7 



[cp-^.j) 
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Thus, the continuity equation will be invariant under the Lorentz transfor- 
mation if 

= a^J-^cpP, (69) 

In order to understand the physical meaning of these equations, consider the 
following example. 



Example 

Assume that in the S frame there is a stationary volume charge of density 

p 7^ 0. Since p is stationary, there are no currents in the S frame ( J = 0). 
What are the charge and current densities as seen in the frame? 



In the S frame 

/=0, p^O 
According to Eq. (69), in the frame 

f = -^cpP , p' = 7P 

Thus, there is a current in the 5" frame. As seen from a given part of the 

charge is length contracted in the direction of motion so the charge density 
is correspondingly increased by the factor 7 > 1. The length contracted 
charge density appears from to move in the opposite direction. We can 
understand this result: The stationary charge in the S frame moves with 
velocity —u in the 5" frame. 



Less obvious and more interesting is the following situation. 

In the S frame 

J^O , p = 0 . 

Then, someone will see charge density 0 in the frame. 
This is a pure relativistic effect. 
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18.2 Transformation of Electric and Magnetic Field 
Components 



To find the transformation rules for electric and magnetic field components 
we will use the transformations of the time and space derivatives derived 
above. 

Consider two of the Maxwell equations that in the S frame are 



V-i) = p 

V xH = J + 



dD 
'dt 



These equations should be equivalent of two equations 



VD' = p' 

V xH' = J' + 



dP' 



in the S' frame. 



Using the transformations of the time and space derivatives, we have 



Q! * V' — 7/3 



d 



X H 



d{ct')_ 



a * V' — 7/3 



d{ct' 

d 

d{ct') 



= J 



cD = cp 



Substituting the transformations of J and p, we find that the D and H 

vectors transform as 



cD' = 7 (a-^ *cD + i3xH) 

H' = ^ xcD + a-^ * h) 
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where a ^ is the inverse of the matrix a 



a 



-1 



l + (^-l)|t 0 
0 l + 

0 



0 



0 
0 

l + (^- 



a- 



From the Maxwell equations 



0 



V X £ = - 



dB 



we find that the E and B vectors transform as 

E' = ^ (a-^ * E + P X cB) 

cB' = j(-PxE + a-^*cB) 



18.3 Transformation Rules in Terms of Parallel and 
Normal Components 

Suppose that the frame S' is moving with speed u in the direction paral- 
lel to the z axis. In this case, = Py = Pz = P ^ ^■iid then the 

transformations take the form 

— ^ /s /s ^ ^ — ^ 

cD = ^cD^i + ^cDyj + cD;,k + ^pk X H 

= -^pck X D + ^Hj + ^Hyj + H;,k 

= jEj + jEyj + E^k + ^cf3kxB 

cB^ = — 7/3fc X ^ + ^cB^i + ^cByj + cB^k 

It is useful to rephrase the transformation rules in terms of components 

parallel and normal to u. The parallel components are the z components and 
the normal components lie in the xy plane. For example 

E = E^ + E\\ = (eJ + Eyj) + E,k 
and the same for H and B, 
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18.3.1 Rules for Parallel Components 



Ei = E'k = E^k = En 

\\ Z Z II 

and similarly 

5|| = E|| , ^|[ = , 5|j = D\\ 
18.3.2 Rules for Normal Components 

/N /\ 

= ^cDJ + icDyj + ^(5Hxj - iPHyi 
= 7 {cD, - (3Hy) I + 7 {cDy + j 

Hi = H^l + Hfj 

= j{H^ + (3cDy) i + ^{Hy- (3cD^) j 

E[_ = Eii + E;^j 

= 7 - (3cBy) i + ^{Ey + (3cB^) j 

cB[ = cB'J^cB;) 

= 7 {cB^ + (3Ey) i + 7 {cBy - PE^) j 

In general 

E[ = ^[e±+0xcBa_) 
cB[ = 7 (cB± -Px E±) 
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Example 1 - purely electric field in S 



Suppose that in E ^ 0 but 5 = 0. 
Then from the transformation rules, in 5": 

En = E[ = ^E± 



b!\=0 b[=-^ux E± 



Thus 

f u X E[_ u X E 



B =5, =- 



since u x En = 0. 



Thus what appears to be purely an electric field to one observer is seen 
as both an electric and a magnetic field to a second observer moving with 
respect to the first. 

Example 2 - purely magnetic field in S 

Now suppose that in S", £^ = 0 while B ^ 0, 
Then using the transformation rules, in S': 

5|J = Bn B[ = jB± 



E|| =0 E^=^uxB± 

Thus 

E' = E'j_ = u X B'j_ = u X B' 

We see that what appears to be a purely magnetic field for one observer 
will appear to be both an electric and a magnetic field to a relatively moving 
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observer. 

This result could be used to calculate the emf in an electric dynamo from 
the point of view of an observer watching the conductor move in a magnetic 
field or from the point of view of an observer moving with the conductor. 

18.4 Transformation of the Components of a Plane EM 
Wave 

Example 1 

Suppose that a plane wave propagates in vacuum along the z axis. Then 



Hence from the transformation rules, in moving in the same direction: 



E 
B 



0 



E' 



7 LEo* - PcBoi 



7 {Eo - uBq) i 
7 {cBo - /3Eo) J 



cB' 



7 (cBoj - pEoj 



Since in vacuum 



cBo = Eo 



we obtain 




and 
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Thus, the ratio = |£^o|/|5o| is constant and independent of u. 



Example 2 

Suppose that the frame moves in the direction of the electric field, i.e. 
u — ui. 



X i 




u 



In this case 



and 



u 



cB ' = ^cBoj 



Thus, the magnetic field remains unchanged but the electric field turns to- 
wards the direction of propagation of the wave. 



18.5 Doppler Effect 



Consider a plane wave propagating in vacuum 

E{t) = Eoe*^"*-^-^^ 
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In the moving frame this wave will have a different frequency cj^ and 
the wave vector A;', but the phase of the wave will remain unchanged as in 
invariant under the transformation 

(j) — (jjt — k ' r — iJi! — y - r^ 

Using the Lorentz transformations 

f — a * r ' + 7/3ct' 
ct = 7(ci' + ^.r') 

Hence 

— y ' — ujt — k ' r 
= — (ct' + /3 • f') - ^ * • r' - ^k' Pet' 



— oj^t' — ^t'k ' u-\- [ — (3 — a * ) • r ' 



c 



— ^ [ijj — k ' u\i! — y ' T 




Thus 



{uj —k ' 



k' = 



a* k —u 
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If the wave propagates in vacuum, k = -k^ and then 



uj = 7u; (1 cos 01 



where (j) is the angle between the direction of propagation of the wave and u. 
For = 0 



c — u 
c + u 



18.6 Transformation of Energy of a Plane EM Wave 



Consider an EM wave propagating in the k direction and an observer moving 
in the z direction, as shown in the Figure 18.6. 




Let E = E i. Then, in the frame 5': 

= ^E + ^PxcB± 

— ^E + 7 cos (j) P X cB 
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Since B = Bj + B\\ k (the wave propagates in the plane yz)^ and cB = 
we obtain 

= 7jB — 7 cos /5c5 i 
= 7(1 -/3cos(/))^ 

Consider now energy of an electric field of a plane EM wave confined in 
a volume AV, In the S frame 

1 



In the frame 



However 



AV = AxAyAz 



Assume that S' frame moves in the direction of the z axis, and the wave 
propagates in the direction k •u = u cos <j). Then 



where 



AV' = AxAyAz' 



Az' = 



7 (1 - ^COS0) 



Hence 



AV 



— - cos 

c 



It is interesting to compare the transformation of energy with the transfor- 
mation of frequency 



= iye7(^l-^C0S(/)^ 

oj = cc;7 (1 COS0J 
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We see that the energy and frequency transform in the similar way, indicating 
that We ~ uj. This proportionahty was predicted in quantum physics as We — 
hu and forms backgrounds of the so called quantum electrodynamics. 



Questions: 

(1) Find the condition under which the continuity equation for p and J 
is invariant under the Lorentz transformation. 
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Appendix A 

PHYS3050 Facts and Formulae 



Gauss' Divergence Theorem : 

s V 

Stokes's Theorem : ^F - di= Jv x F-fidS 

i s 

Numerical values in SI units: 
So = 8.85 X 10-12 
/io = 47r X 10"'^ 

c = 3 X 10^ [ms-i] 

For the electron: e = 1.6 x 10"^^ [C], m = 9.11 x 10"^^ [kg] 



The Lorentz force law : F — q{E -\-v x B) 



Coulomb's Law : F — — ^^^f 

—*■ jjjQ I d^ X v 

Biot — Savart Law : dB — — 

47r 

Gauss' Law : (f E - hdS = — 

J So 

Ampere's Circuital Law : ^ ^ ' dl— fiol 



Maxwell's Equations in vacuum: 



dB ^ 1 dE 
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Maxwell's Equations in material bodies: 

VD = p VH = 0 



V X E = - 



dB 

'at 



dP 

'at 



Poynting vector: 



Poynting's Theorem: 



N = Sqc^E X B 



isoc^{E xB)-hdS=- J E- JdV [ ( IsoE^ + ^— I dV 



s 



V 



dtJ \2 

V 



2 /io 



Or in polarizable materials where it is convenient to define D = SqE + P and 



H = ^-M: 

MO 



({ExH'hdS=- J H'^dV-jE~dV-J E ^ JdV 



s 



V 



V 



V 



A theorem on the calculation of the mean Poynting vector from complex 
fields: 

1 



7V= - Re 
2 



E^xH 



The rate of doing work in magnetization 



dW dB 
= H- 



dt 



dt 



and B = fio{H + M). 



An arbitrary vector F can be written as: 
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Fields and potentials: 



E = -V$- 



M 
'dt 



B = V X A 



Differential equation for the vector potential: 



V^A - — ^ = -/xoJ + V V • A + — 



C2 



c2 



In the Lorentz gauge: 



V • A = 

c2 



the difTerential equations for the electromagnetic potentials are; 



1 d^<^ 



V $ 

V A 

and these have solutions of the form: 



^ _P_ 

So 



$ = 



1 



A = 



1 



-dV 



r 



The field of a Hertzian dipole: 



Er 



2 2 

+ ^ 



J.{ijjt—kr) 



Eq = 



/oA£sin^ 

AtTSqC 



1 1 iA: 

+ ^ + — 



J,{ujt—kr) 



Ba — 



loAisine 



ik 



1 



- + -2 



J,{Ljt—kr) 
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The mean energy flux from the Hertzian dipole: 



SEqC 



Ai 

T 



2 



A series solution in 2 dimensions to Laplace's equation in Cartesian co- 
ordinates: 

^(x^z) = ^[A/c sm{ax) + Bk cos{ax)][Ck sinli(a^) + Dk cosh(az)] 
k 

A series solution in 3 dimensions to Laplace's equation in spherical polar 
coordinates: 



"^[aim cos{m(f)) + bim sin(m0)] (cos 6*) 



m 



Useful properties of trigonometrical functions: 



sin(a ± P) 
cos(a ± P) 



sin^ a 



cos a = 



sin a cos ± sin cos a 

cos a cos P =F sin a sin 

-(1 — cos 2a) 

1/ 

-(1 + cos 2a) 



Jo 



sin^ 9 d9 = - 



4 
3 



/•27r 

/ 8m(m(f)) sin(n<^) (i<^ = 

/•27r 

/ cos(m0) cos(n<j!)) = 
Jo 



0 for m ^ n 

TT for m = n 

0 for m ^ n 

TT for m = n 



/'27r 

/ sin(m^) cos(n^) d(j) = 0 for all m and n 

^0 
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Properties of Legendry polynomials: 



1 



unless m — n and / = k 



-1 



1 



2 (Z + m)! 



-1 



21 + 1 {I- m)\ 



Po = l 



Pi° = cos e Pl = sin 6 P2=\ (3 cos(2^) + 1) 

Pl = lsin{2e) P| = ^(l-cos(2^)) 
P;(l) = l for all I 



A theorem on the electrostatic potential due to a distribution of electric 

dipoles of moment per unit volume P: 



47r6:o J ^ 47r6:o ^ ^ 

s V 



dV 



A theorem on the vector potential due to a distribution of magnetic 
dipoles of moment per unit volume M: 



1 



V 



V X M 1 

dV + 



M xh 



s 



dS 



A dispersion equation: 



k = u 



e fjL — fi 



II a 
e +- 

00. 



II a 
e +- 

00. 



I II 



+ e fj, 
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The skin depth in a good conductor: 



General boundary conditions: 

— * 

• The normal component of B is continuous across a boundary. 

• The normal component of D is continuous across a boundary. 

• The tangential component of E is continuous across a boundary. 

• The tangential component of H is continuous across a boundary. 

Special boundary conditions at the surface between a dielect 
and a perfect conductor: 

• The tangential component of ^ = 0. 

• The normal component oi B or H — 0, 
The characteristic impedance of free space: 



Zo = ^ = J— 377 - 1207r ohms 

H y So 



The Fresnel equations: 
Case 1: E normal to the plane of incidence. 

Reflection: 



£/l — / -OjQ 



Transmission: 





cos 6i — 111 




— k\ sin^ 9i 




cos 9i + /ii 




— k\ sin^ 6i 




2kifi2 


cost 





E2 = / Eq 

kifi2 cos 0i + fJ^iy k2 — kl sin^ 6i 
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Case 2: E in the plane of incidence. 

Reflection : 

fxikl cos 9i - iJi2ki\Jkl - kl svd? Oi 
til = / == 

/ii^l cos 9i + /X2fci Y ~ ^1 ^2 

Transmission: 

2/xifc|cos(9i 

/xi^l cos 9i + /^2^i v ^2 ~ ^1 

In dielectric media the Fresnel equations become: 
E normal to the plane of incidence: 





El 






Eq 


sin(^i + Ot) 




E2 


2 cos 6i sin^t 




Eq 


sin(^i + ^t) 


E in the plane of incidence: 








El 


_ tan(^i - dt) 






~ tan(^i + Ot) 


E2 




2 cos 6i cos 6t 


Eq 


sin 


{6i + Ot) cos(^i 



Rectangular waveguides: 

If the fields vary as e~^^ : 
For TE modes the longitudinal component of H satisfies: 

where 

A;^ = 7^ — iu;fi{a + iojs) 

Then satisfying the boundary conditions (assuming the walls are perfect con- 
ductors) requires: 

7 = I — ) + I ] +^/^^(^ + ^^^) 
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For the lossless waveguide: 
Cut-off frequency for the mn mode: 



Phase velocity: 



Group velocity: 




VECTOR FORMULAS 



V($ + 


= V^ + V* 




= V-A + V-B 


V X {A + B) 


= V X A + V X B 


V($^') 


= + 


V • ($1) 


= .4-V$ + $V-A 




= B ■ {W X A) - A- {W X B) 


V X {^A) 


= V$ X i'+ $V X A 


V x{AxB) 


= AV ■ B - BV ■ A+ {B ■V)A- {A-V)B 


V- V$ 


= 


V • (V X A) 


= 0 


V X V$ 


= 0 


V X (V X 1) 


= V(V-l)-V^i 


1- {B X C) 


= B-(CxA) = C-iAxB) 


Ax{BxC) 


= B{A-C)-C{A-B) 
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FORMS OF VECTOR OPERATIONS IN CYLINDRICAL 

COORDINATES 



V$ = 



dp p d(j) dz 



VA = - 



p dp P d<f> dz 



V X A = p - 



1 OA, OA^ 



p d4> dz 



+ 0 



fdAo OA 



\ dz 



dp 



,, lfd{pA^) dA, 



P\ dp 



d<j) 



= - 



1 5 / 5$ \ 1 

P— I + — — + 



p dp \ dp I p^ d^ dz^ 



VECTOR AND DIFFERENTIAL OPERATIONS IN 

SPHERICAL COORDINATES 



A = {Ax sin 6 cos ^-\- Ay sin 6 sin + Az cos 6) r 

+ {Ax cos 6 cos <!) + Ay cos 6^ sin ^ — yl;^ sin 6) 9 

+ {—Ax sin 0 + Aj, cos 0) 4> 

= A^ f + Ae ^ + A^ 0 
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dr r dO r sin 9 d(j) 



dr 



r sin^ 



de 



rsin^ d(j) 



VxA = 



+ 



r sin 6' 



dismOA^) _ dAe 



e 

+ - 

r 



0 



dr 



de 



1 d{rA^) 



SlYld d<j) 



dr 



= — 



1 d 



d^ 

dr 



+ 



1 



sin 6 dO 



d ( . d^ 



+ 



1 



sin 6* c^^^ 
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